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Feshbach resonances are the essential tool to control the interaction between atoms in ultracold 
quantum gases. They have found numerous experimental applications, opening up the way to 
important breakthroughs. This Review broadly covers the phenomenon of Feshbach resonances in 
ultracold gases and their main applications. This includes the theoretical background and models 
for the description of Feshbach resonances, the experimental methods to find and characterize 
the resonances, a discussion of the main properties of resonances in various atomic species and 
mixed atomic species systems, and an overview of key experiments with atomic Bose-Einstein 
condensates, degenerate Fermi gases, and ultracold molecules. 
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I. INTRODUCTION 

A. Ultracold gases and Feshbach resonances: Scope of the 
Review 

The great impact of ultracold atomic and molecular 
quantum gases on present-day physics is linked to the 
extraordinary degree of control that such systems of- 
fer to investigate the fundamental behavior of quantum 
matter under various conditions. The interest goes be- 
yond atomic and molecular physics, reaching far into 
other fields, like condensed matter, few- and many-body 
physics. In all these applications, Feshbach resonances 
represent the essential tool to control the interaction be- 
tween the atoms, which has been the key to many great 
breakthroughs . 

U ltracold gases are generally produc ed by laser cool- 



ing ( Metcalf and van der Stratenl.ll999l) and subsequent 
evaporative cooling ( Ketterle and van Drutenl . ll997t ). At 



temperatures in the nanokelvin range and typical number 
densities somewhere between 10^^ cm^'^ and 10^^cm~'^, 
quantum-degenerate states of matter are formed when 
the atomic de-Broglie wavelength exceeds the typical in- 
terparticle distance and quantum statistics governs the 
behavior of the system. The attainment of Bose- Einstein 
condensation (BEC) in dilute ultracold gases marked the 
start i ng point of a new er a in physics ( Anderson et all . 
119951: iBradlev et d\ . Il995t iDavis et all llQQSl ). and"de- 
generate at omic Fermi gas e s entered the stage a few 
years later foeMarco et all Il999l : ISchreck et all . l200ll: 
iTruscott et all l200l[ ). The developments of the tech- 
niques to cool and trap atoms by laser ligh t were rec- 
ognized with the 1 997 Nobel priz e in ph ysics ( Chul . [l998t 
ICohen-Tannoudiil , [l998: PhiUig^, Ulii) . Only four years 
later, the achievement of BEC in dilute gases of al- 
kali atoms and early fundamental studies of the prop- 
erties of the condensates led to the 2001 Nobel prize 
(ICornell and Wiemanl . [200l iKetterld . I2002D \ 

In this Review, we give a broad coverage of Feshbach 
resonances in view of the manifold applications they have 
found in ultracold gases. Regarding theory, we focus on 



^ For overviews on laser cooling and trapping, BEC, and ultracold 
Fermi gases see the proc eedings of the Va r enna summer schools 
in 1991, 19 98, and 2006 l lArimondo et ali Il992l: llnguscio eroH . 
l2008l. Il999l 'l. For reviews on the theory of degenerate quan- 
tum gases of bosons and fermions see llDalfovo et al\ , Il999 



the underlying two-body physics and on models to de- 
scribe Feshbach resonances. In the experimental part we 
include applications to few- and many-body physics; we 
discuss typical or representative results, instead of the 
impossible attempt to exhaustively review all develop- 
ments in this rapidly growing field. Several aspects of 
Feshbach resonances and related topics have already been 
reviewed elsewhere. An e arly review on Feshbach r eso- 
nance theory was given by ( Timmermans et all . 
another theoretical review, ( Duine and Stooi 



2004) fo- 



1999 ). In 



cusse d on atom-molecule c ohere nce. ( Kohler et all . \206di ) 



and llGi orgini et al\. I200S ), res pecti vely, and the textbooks b y 
llStringari and Pitaevskiil . boOSl) and llPethick and Smithl. [200^) . 



and l|Hutson and Soldanl . l2006l) reviewed the formation 
of ultracold molecules near Feshbach resonances. The 
closely related topi c of photoassociation was reviewed by 
(| Jones et flZ.I . l2006t ). 

In Sec. ini we start with a presentation of the theoreti- 
cal background. Then, in Sec. lIIII we present the various 
experimental methods to identify and characterize Fesh- 
bach resonances. There we also discuss the specific inter- 
action properties of different atomic species, which can 
exhibit vastly different behavior. In Sec. IIVI we present 
important applications of interaction control in exper- 
iments on atomic Bose and Fermi gases. In Sec. El we 
discuss properties and applications of ultracold molecules 
created via Feshbach association. Finally, in Sec. IVII we 
discuss some related topics, like optical Feshbach reso- 
nances, interaction control in optical lattices, few-body 
physics, and the relation to molecular scattering reso- 
nances and cold chemistry. 



B. Basic physics of a Feshbach resonance 

The physical origin and the elementary properties of a 
Feshbach resonance can be understood from a simple pic- 
ture. Here we outline the basic ideas, referring the reader 
to Sec. inifor a more detailed theoretical discussion. 

We consider two molecular potential curves Vbg(i?) and 
Vc{R), as illustrated in Fig. [TJ For large internuclear 
distances R, the background potential Vbg{R) asymptot- 
ically connects to two free atoms in the ultracold gas. 
For a collision process, having the very small energy E, 
this potential represents the energetically open channel, 
in the following referred to as the entrance channel. The 
other potential, representing the closed channel, 

is important as it can support bound molecular states 
near the threshold of the open channel. 

A Feshbach resonance occurs when the bound molec- 
ular state in the closed channel energetically approaches 
the scattering state in the open channel. Then even weak 
coupling can lead to strong mixing between the two chan- 
nels. The energy difference can be controlled via a mag- 
netic field when the corresponding magnetic moments are 
different. This leads to a magnetically tuned Feshbach 
resonance. The magnetic tuning method is the com- 
mon way to achieve resonant coupling and it has found 
numerous applications, as will be extensively discussed 
in this Review. Alternatively, resonant coupling can be 
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FIG. 1 Basic two-channel model for a Feshbach resonance. 
The phenomenon occurs when two atoms colliding at energy 
E in the entrance channel resonantly couple to a molecular 
bound state with energy Ec supported by the closed channel 
potential. In the ultracold domain, collisions take place near 
zero-energy, i5 ^ 0. Resonant coupling is then conveniently 
realized by magnetically tuning E^ near 0, if the magnetic 
moments of the closed and open channel differ. 
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FIG. 2 Scattering length a (Panel (a)) and molecular state en- 
ergy E (Panel (b)) near a magnetically tuned Feshbach reso- 
nance. The binding energy is defined to be positive, £(, = —E. 
The inset shows the universal regime near the point of reso- 
nance where a is very large and positive. 



achieved by optical methods, leading to optical Feshbach 
resonances with many conceptual similarities to the mag- 
netically tuned case; see Sec. IVI.Al Such resonances 
are promising for cases where magnetically tunable reso- 
nances are absent. 

A magnetically tuned Feshbach resonance can be 
describe d by a sim ple expression^, introduced by 
(Moerdiil Tet fflZ.I . Il995l ). for the s-wave scattering length 
a as a function of the magnetic field B, 



a{B) = Obg 1 



A 



B- Br 



(1) 



Figure [UJa) illustrates this resonance expression. The 
background scattering length Obg, which is the scatter- 
ing length associated with Vbg(i?), represents the off- 
resonant value. It is directly related to the energy of the 
last-bound vibrational level of Vhg{R)- The parameter 
denotes the resonance position, where the scattering 
length diverges (a — > ±cx)), and the parameter A is the 
resonance width. Note that both abg and A can be posi- 
tive or negative. An important point is the zero crossing 
of the scattering length associated with a Feshbach res- 
onance; it occurs at a magnetic field B = Bq -\- A. Note 
also that we will use G as the magnetic field unit in this 
Review, because of its near-universal usage among groups 
working in this field; 1 G = 10^^ T. 

The energy of the weakly bound molecular state near 
the resonance position Bq is shown in Fig.l^^b), relative 



^ This simple expression applies to resonances without inelastic 
two-body channels. Some Feshbach resonances, especially the 
optical ones, feature two-body decay. A more general discussion 
including inelastic decay is given in Sec. III. A. 31 



to the threshold of two free atoms with zero kinetic en- 
ergy. The energy approaches threshold at E = Q on the 
side of the resonance where a is large and positive. Away 
from resonance, the energy varies linearly with B with a 
slope given by Sji, the difference in magnetic moments of 
the open and closed channels. Near resonance the cou- 
pling between the two channels mixes in entrance-channel 
contributions and strongly bends the molecular state. 

In the vicinity of the resonance position at Bq, where 
the two channels are strongly coupled, the scattering 
length is very large. For large positive values of a, a 
"dressed" molecular state exists with a binding energy 
given by 



Eh 



(2) 



where fi is the reduced mass of the atom pair. In this limit 
Eh depends quadratically on the magnetic detuning B — 
Bq and results in the bend seen in the inset to Fig. (5] This 
region is of particular interest because of its universal 
properties; here the state can be described in terms of 
a single effective molecular potential having scattering 
length a. In this case, the wavefunction for the relative 
atomic motion is a quantum halo state which extends to 
a very large size on the order of a; the molecule is then 
called a halo dimer, see Sec. IV.B.2I 

A very useful distinction can be made between reso- 
nances that exist in various systems; see Sec. III. B. 21 For 
narrow resonances with a width A typically well below 
1 G (see Appendix) the universal range persist only for 
a very small fraction of the width. In contrast, broad 
resonances with a width typically much larger than 1 G 
tend to have a large universal range extending over a 
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FIG. 3 Observation of a magnetically tuned Feshbach reso- 
nance in an optically trapped BEG of Na atoms. The upper 
panel shows a strong loss of atoms near the resonance, which is 
due to enhanced three-body recombination. The lower panel 
shows the dispersive shape of the scattering length a near the 
resonance, as determined from measurements of the mean- 
field interaction by expansion of the condensate after release 
from the trap; here a is normalized to the background value 
abg. The magnetic field is given in G, where IG = 10"* T. 
Repr inted by perrn i ssion from Macmillan Publishers Ltd: Na- 
ture (|lnouve et all\l99^ ]. copyright 1998. 



considerable fraetion of the width. The first class of res- 
onances is referred to as closed-channel dominated res- 
onances, whereas the second class is called entrance- 
channel dominated resonances. For the distinction be- 
tween both classes, the width A is not the only relevant 
parameter. Also, the background scattering length Obg 
and the differential magnetic moment dfi need to be taken 
into account. Section Hi. B .21 discusses this important dis- 
tinction in terms of a dimensionless resonance strength. 

Figurc[3]shows the ob servation of a Feshba ch resonance 
as reported in 1998 by ( Inouve et all llQQSl ) for an opti- 
cally trapped BEC of Na atoms. This early example 
highlights the two most striking features of a Feshbach 
resonance, the tunability of the scattering length accord- 
ing to Eq. [1] and the fast loss of atoms in the resonance 
region. The latter can be attributed to strongly enhanced 
three-body recombination and molecule formation near a 
Feshbach resonance; see Sec. IIII.A.21 

A Feshbach resonance in an ultracold atomic gas can 
serve as a gateway into the molecular world and is thus 
strongly connected to the field of ultracold molecules; 



see Sec. [V] Various techniques have been developed to 
associate molecules near Feshbach resonances. Ultracold 
molecules produced in this way are commonly referred to 
as Feshbach molecules. The meaning of this term is not 
precisely defined, as Feshbach molecules can be trans- 
ferred to many other states near threshold or to much 
more deeply bound states, thus being converted to more 
conventional molecules. We will use the term Feshbach 
molecule for any molecule that exists near the thresh- 
old in an energy range set by the quantum of energy for 
near-threshold vibrations. The universal halo state is a 
special, very weakly bound case of a Feshbach molecule. 



C. Historical remarks 

Early investigations on phenomena arising from the 
couplin g of a bound state to the continuum go back to the 
1930s. (|Ricel . ll93,'^ consid ered how a bound state predis - 
sociates into a continuum, ( Fanol . [l~935l ; lFano et al. L I2005D 
described asymmetric line shapes occurring in s uch a sit- 
uatio n as a result of quantum interference, and ( Beutleil . 
Il935h reported on the observation of highly asymmetric 
lineshapes in rare gas photoionization spectra. Nuclear 
physicists considered basically the same situation, having 
nuclear s cattering experiments in m ind instead of atomic 
physics. ( Breit and Wigneil . Il936t ) considered the situa- 
tion in the limit when the bound state plays a dominant 
role and the asymmetry disappears. Later, interference 
and line-shape asymmetry were taken into account by 
several authors ( Blatt and Weisskopi 1195^ ). 

Herman Feshbach (1917-2000) and Ugo Fano (1912- 
2001) developed their thorough treatments of the res- 
onance phenomena that arise from the coupling of a 
discrete state to the continuum. Their work was car- 
ried out independently, using different theoretical ap- 
proaches. While Feshbac h's work o r iginat e d in t he con- 
text of nuclear physics ( Feshbachl . Il958l Il962[ ). Fano 
approac hed the prob lem on the background of atomic 
physics ( Fand 1961), reformulating and e xtending his 
earlier work (iFanol . 119351 ; iFano et all HqoI). Nowadays, 
the term "Feshbach resonance" is most widely used in the 
literature for the resonance phenomenon itself, but some- 
times also the term "Fano-Feshbach resonance" appears. 
As a curiosity Feshbach himself considered his name be- 
ing attached to a well-know n resonance phenomenon as a 
mere atomic physics jargon ( Kleppneil . 120041 ; [Raul . |2005[ ) . 
Fano's name is usually associated with the asymmet- 
ric lineshape of such a resonance, well known in atomic 
physics as a "Fano profile" . 

A prominent example for the observation of a Fesh- 
bach resonance in a tomic physics is the experiment of 
( Bryant et all Il977( ) on photodetachment by the nega- 
tive ion of hydrogen. Near a photon energy of 11 eV 
two prominent resonances were seen, one of them be- 
ing a Feshbach resonance and the other one a "shape 
resonance"; see Sec. III. A. 31 Many more situations 
where Feshbach resonances play an important role can 
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FIG. 4 Number of publications per year (from 1985 to 2008) 
with "Feshbach resonances" appearing in the title. Data from 
ISI Web of Science. 



up many new applications in the field. This is reflected 
in the steep increase of the publication rate in the period 
from 2002 to 2004. 



II. THEORETICAL BACKGROUND 

This review primarily concentrates on magnetically 
tunable resonances, described in detail in the next Sec- 
tions, while Section IVI.AI discusses optical changes in 
scattering lengths. Here we describe the two-body 
physics of collision resonances, not the few-body or many- 
body aspects. Properties of a number of magnetic Fesh- 
bach resonances are tabulated in the Appendix; see Ta- 
ble EVl 



A. Basic collision physics 



be f ound in atomic, molecula r , and chemical physics ; 
see (Gauvacg and Herzenberg ^ 19821: MacArthur et ai . 



1985 



1975 



.Nieh and Valentiiiil . Il99u iSpence and Noguchi , 
IWeber et all 19991 ) for a few examples. In such 



experiments, the resonances occur when the scattering 
energy is varied. This is in contrast to the experiments 
on ultracold gases, where scattering takes place in the 
zcro-cncrgy limit and the resonances occur when an ex- 
ternal field tunes bound states near threshold. 

In the context of qua ntum gases, Fes hbach resonance 
were first considered by ( Stwallevl . Il976t ). who suggested 
the existence of magnetically induced Feshbach reso- 
nances in the scattering of spin-polarized hydrogen and 
deuterium atoms (H-l-D and D-l-D). He pointed to en- 
hanced inelastic decay near these resonances and sug- 
gested that they should be avoided to maintain sta- 
ble spin-polarized hydrogen gases. A related loss res- 



onance in hydrogen was observed by (jRevnolds et al\ . 
Il986l ). The pos itive aspect of s uch r esonances was first 
pointed out by ( Tiesinga et aLl . [19931 ). who showed that 
they can be used to change the sign and strength of the 
interaction between ultracold atoms. In 1998, the possi- 
bility of interacti on tuning via Feshb ach resonances was 
demonstrated by ( Inouve et aZ.l . ll998l ) for a ^^Na BEG, as 
alread y discussed in the prec eding Section. In the same 
year, ( Courteille et all Il998l ) demonstrated a Feshbach 
resonance in a trapped sample of *^Rb atoms through 
the enhancement of photoassociative loss induced by a 
probe laser. 

The very important role of Feshbach resonances in 
present-day quantum gas experiments can be highlighted 
by looking at the number of publications per year with 
"Feshbach resonances" in the title; see Fig. U) Before 
1998, one finds just a few publications with the great ma- 
jority not related to ultracold atoms. Then, after 1998, 
a substantial increase is observed as a result of the first 
successful experiments with Feshbach resonances in ul- 
tracold gases. It then took a few years until Feshbach res- 
onances had become a fully established tool and opened 



The theory for describing 2-body collisions is de- 
scribed in a number of textbooks ( Messiahl . Il966t 
iMott and Massevl . Il965l : iTavlod . Il972t) . Let us first con- 
sider the collision of two structureless atoms, labeled 1 
and 2 with masses mi and 7712 interacting under the in- 
fluence of the potential V^(R), where R is the vector be- 
tween the positions of the two atoms with magnitude R. 
The separated atoms are prepared in a plane wave with 
relative kinetic energy E = h^k"^ /{2^) and relative mo- 
mentum ?ik, where ^ = mim2/(mi + 7712) is the reduced 
mass of the pair. The plane wave in turn is expanded in a 
standard sum over spherical harmonic functions Yim, (R), 
where £ is the relative angular momentum, mi is its pro- 
jection along a space-fixed z axis, and R = R/J? is the 
direction vector on the unit sphere ( Messiahl . Il966l ). This 
expansion is called the partial wave expansion, and the 
various partial waves ^ = 0,1,2,... are designated s-, p-, 
d-, ... waves. 

If the potential T^(R) is isotropic, depending only on 
the magnitude of R, there is no coupling among partial 
waves, each of which is described by the solution ipe{R) = 
(pi{R)/R to the Schrodinger equation 



h"^ <f(t>i{R) 

2JI dR^ 



+ Vi{R)MR) = EMR), (3) 



where Ve{R) = V{R) + h^e{£ + l)/{2fiR^) includes the 
centrifugal potential, which is repulsive for £ > and van- 
ishes for the s-wave. We assume V{R) — > as i? — > 00, so 
that E represents the energy of the separated particles. 
This equation has a spectrum of Ng bound state solu- 
tions at discrete energies Eni for E < and a continuous 
spectrum of scattering states with E' > 0. While bound 
states are normally labeled by vibrational quantum num- 
ber w = 0, . . . , iVf — 1 counting up from the bottom of the 
potential, we prefer to label threshold bound states by 
quantum number n = —1,-2,... counting down from 
the top of the potential for the last, next to last, etc., 
bound states. The bound state solutions \n£) are normal- 
ized to unity, \{n£\n£)\'^ = 1, and 0„£(i?) = {R\n£) 
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as i? oo. The scattering solutions, representing the 
incident plane wave plus a scattered wave, approach 



,{R,E) 



sm{kR-Tr£/2 + r]e{E)) 



(4) 



as i? ^ oo, where rji{E) is the scattering phase shift 

and c = \J 2/x/(7r?i^) is a constant that ensures the wave 
function \Ei) is normalized per unit energy, {E£\E'i) = 
(l>*i{R,E)(j)e{R,E')dR = 5{E - E'). The scattering 
phase shift is the key parameter that incorporates the 
eff ect of the whole potentia l on the collision event. 

(jSadeghpour et al\ . [2003) reviews the special proper- 
ties of scattering phase shift near a collision threshold 
when fc ^ 0. If V{R) varies as l/R" at large i?, then 
ian-qi oc fc^^+^ if 2^ + 1 < s — 2 and tanyy^ oc fc'*~^ if 
2£ + 1 > s — 2. While Levinson's theorem shows that 
rji —^ NiTT as fc 0, we need not consider the NgTr part 
of the phase shift in this review. For van der Waals po- 
tentials with s = 6, the threshold tanrii varies as k and 
for s- and p-waves, and as fc^ for all other partial waves. 
The properties of s-wave collisions are of primary interest 
for cold neutral atom collisions, where near threshold, a 
more precise statement of the variation of tan rjo with k 
is given by the effective range expansion. 



kcotrio{E) = - - + ^rofc^ . 

a 2 



(5) 



where a is called the s-wave scattering length and tq 
the effective range. For practical purposes, it often suf- 
fices to retain only the scattering length term and use 
tan77o(i?) = —ka. Depending on the potential, the scat- 
tering length can have any value, — oo < a < -|-cx3. 

When the scattering length is positive and sufficiently 
large, that is, large compared to the characteristic length 
scale of the molecular potential (see Section riI.B.l[) . the 
last s-wave bound state of the potential, labeled by index 
n ~ —1 and £ = 0, is just below threshold with a binding 
energy Ei, = — i?_i,o given by Eq. ^ in the Introduc- 
tion. The domain of universality, where scattering and 
bound state properties are solely characterized by the 
scattering len gth and mass, are discussed a t length in re- 
cent reviews ( Braaten and Hammeil [200I lK6hlcr et all . 
1200611 . The universal bound state wave function takes 
on the form 0_i.o(i?) = •\/2/aexp (—R/a) at large R. 
Such a state exists almost entirely at long range beyond 
the outer classical turning point of the potential. Such a 
bound state is known a s a "h alo state," also studied in nu- 
clear physics f|Riisagertll994l ) and discussed in Sec. lV.B.2l 



1. Collision channels 

The atoms used in cold collision experiments generally 
have spin structure. For each atom i = 1 or 2 in a colli- 
sion the electronic orbital angular momentum L.^ is cou- 
pled to the total electronic spin angular momentum 



to give a resultant j^, which in turn is coupled to the nu- 
clear spin li to give the total angular momentum f^. The 
eigenstates of each atom are designated by the compos- 
ite labels Qi. At zero magnetic field these labels are fiirii, 
where rrii is the projection of f^. For example, alkali- 
metal atoms that are commonly used in Feshbach reso- 
nance experiments, have ^Si/2 electronic ground states 
with quantum numbers Li ~ and Si ~ 1/2, for which 
there are only two values of /i = — 1/2 and li + 1/2 
when li ^ 0. Whether fi is an integer or half an odd in- 
teger determines whether the atom is a composite boson 
or fermion. 

A magnetic field B splits these levels into a manifold 
of Zeeman sublevels. Only the projection along the 
field remains a good quantum number, and B = levels 
with the same rrii but different fi can be mixed by the 
field. Even at high field, where the individual fi values 
no longer represent good quantum numbers, the fi value 
still can be retained as a label, indicating the value at 
B ^ with which the level adiabatically correlates. 

Figure [5] indicates the Zeeman energy levels versus B 
for the ^Li ator n, a fermion, according to the classic Breit- 
Rabi formula ( Breit and Rabil . Il93lh . The two fi levels 
are split at _B = by the hypcrfinc energy, E'hf /h = 228 
MHz. At large fields the lower group of three levels are 
associated with the quantum numbers ms = —1/2, while 
the upper group has 7715 = -1-1/2. The figure also shows 
our standard notation for atomic Zeeman levels for any 
species and any field strength. We label states by lower 
case Roman letters a, 6, c, . . . in order of increasing en- 
ergy. Some authors prefer to label the levels in order 
numerically as 1, 2, 3, . . . The notation qi can symboli- 
cally refer to the firrii, alphabetical, or numerical choice 
of labeling 

The collision event between two atoms is defined by 
preparing the atoms in states qi and (72 while they are 
separated by a large distance R, then allowing them to 
come together, interact, and afterwards separate to two 
atoms in states q[ and If the two final states are the 
same as the initial ones, 91,(72 = 91,921 the collision is 
said to be elastic, and the atoms have the same relative 
kinetic energy E before and after the collision. If one 
of the final states is different from an initial state, the 
collision is said to be inelastic. This often results in an 
energy release that causes a loss of cold atoms when the 
energetic atoms escape from the shallow trapping poten- 
tial. We will concentrate primarily on collisions where 
the two-body inelastic collision rate is zero or else very 
small in comparison to the elastic rate, since this corre- 
sponds in practice to most cases of practical experimental 
interest. This condition is necessary for efficient evapo- 
rative cooling or to prevent rapid decay of the cold gas. 
Section IIII.A.2I discusses how atom loss due to 3-body 
collisions can be used to detect the presence of 2-body 
resonances. 

In setting up the theory for the collision of two atoms, 
the scattering channels are defined by the internal states 
of the two atoms 1 and 2 and the partial wave, \a) = 
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FIG. 5 Atomic energy levels of the ®Li atom, which has S — 
1/2, I ^ 1, and / = 1/2 and 3/2. The figure shows both 
the projection m of f and the alphabetic shorthand notation 
qi = a, b, c, d, e, and / used to label the levels in order of 
increasing energy. 



\qiq2)\£mi), where {R\£mi) = YimiCR)- Since for col- 
lisions in a magnetic field the quantum number M = 
mi + 7712 + ni^i is strictly conserved, a scattering channel 
can be conveniently labeled by specifying the set of quan- 
tum numbers {qiq2iM}. For s-waves, where £ = = 
and M ~ mi + m2, it is only necessary to specify the 
quantum numbers {9192} to label a channel. 

When the two atoms are of the same isotopic species, 
the wave function must be symmetric (antisymmetric) 
with respect to exchange of identical bosons (fermions). 
We assume such symmetrized and normalized functions, 
as described by ( Stoof et~cil\ . Il988l ) . Exchange symmetry 
ensures that identical atoms in identical spin states can 
only collide in s-, d-, . . . waves for the case of bosons and 
in p-, /- ... waves in the case of fermions; in all other 
cases, collisions in all partial waves are allowed. 

The channel energy Ea ~ E{qi) + E{q2) is the internal 
energy of the separated atoms. Let us assume that the 
atoms are prepared in channel a with relative kinetic 
energy E so that the total energy is Etot = Ea + E. Any 
channel (3 with Ep < Etot is called an open channel, and 
any channel with Ep > Etot is called a closed channel. 
A collision can produce atoms in an open channel after 
the collision, but not in a closed channel, since the atoms 
do not have enough energy to separate to the product 
atoms. 



2. Collision rates 

The partial collision cross section for starting in open 
channel a with relative kinetic energy E and ending in 
open channel (3 can be expressed in terms of the Sa.piE) 



element of the multichannel unitary scattering matrix S. 
The cross section for elastic scattering at energy E in 
channel a is 



(6) 



whereas the unitarity property of S allows us to express 
the cross section for loss of atoms from channel a as 



Crioss,aiE) = 5a 77 (1 - \Sa,aiE)\" 



(7) 



The corresponding partial elastic and inelastic rate coef- 
ficients Kci,a{E) and Kioss,a{E) are found by multiply- 
ing these partial cross sections by the relative collision 
velocity v = hk/^. The factor = 1 except for cer- 
tain special cases involving identical particles. The fac- 
tor Qa = 2 for describing thermalization or inelastic col- 
lisions in a normal Maxwellian gas of two atoms of the 
same species in identical spin states. Inelas tic decay of a 
pure Bose- Einstein cond ensate has go, = 1 ( Kagan et al\ . 
119851 : IStoof et aZ.l . ll989D . 

If only one open channel a is present, collisions are 
purely elastic and Sa.a{E) = eyip{2ir]a{E)) . For s-waves 
the real- valued tanrjaiE) rhOj(y clS fc — > and 

is the scattering length for channel a. When other 
open channels are present, the amplitude \Sa,OL{E)\ is 
no longer unity, and for s-wave we can represent the 



-kcLry for k 



in terms of 



complex phase rja{E) - 
a complex scattering len gth ( Balakrishnan ef"all . Il997t 
iBohn and Juliennd . Il996[) 

Oa = - iba , (8) 

where a and h are real, and 1 — \Sa.a{E)\'^ Akba > 
as k 0. The threshold behavior is 

cr^laiE) = 4TTga{ai + bl) , (9) 

for the s-wave elastic collision cross section and 

^1oss,q(£') = gaba ■ (10) 

M 

for inelastic collisions that remove atoms from channel a. 
Both (Tci,a and /irioss,a approach constant values when E 
is sufficiently small. 

The unitarity property of the ^-matrix also sets an up- 
per bound on the cross sections. Since there is a rigorous 
upper bound of \Sa^a{E)\ < 1, we find that the elastic 
scattering cross section is maximum 



CTc1,q(£^) 



47r 



(11) 



for any channel a (and thus any partial wave £) when 
Sa,a{E) = —1. Furthermore, (Tioss.q(-E'), if nonvanish- 
ing, has a maximum value of (7ioss.q(-E') = daT^/k'^ when 
Sa,a{E) = 0. These limits are called the unitarity limits 
of the cross sections. For s-wave collisions this limit is 
approached at quite low energy approximately equal to 
E w ?i^/(2/ia„), where fca^ w 1. 
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In order to compare with experimental data the par- 
tial rate coefficients must be summed over partial waves 
and thermally averaged over the distribution of relative 
collision velocities at temperature T. This defines the 
total rate coefficients Kc\^qiq2{T) and K\oss,qiq2 

(T) when 

the atoms are prepared in states qi and (72, respectively. 
Often the temperatures arc sufficiently small that only 
the s-wave entrance channel contributes. 



3. Resonance scattering 

The idea of resonance scattering in atomic and molec- 
ular systems has been around since the earliest days of 
quantum physics, as described in the introduction. A 
conventional "resonance" occurs when the phase shift 
changes rapidly by w tt over a relatively narrow range 
of energy, due to the presence of a quasibound level of 
the system that is coupled to the scatttering state of 
the colliding atoms. Such a resonance may be due to 
a quasibound level trapped behind a repulsive barrier 
of a single potential, or may be due to some approxi- 
mate bound state which has a different symmetry and 
potential from that of the colliding atoms. The former 
is commonly known as a "shape resonance" , whereas 
the latter is often called a "Fcshbach resonance" , in 
honor of Herman Fcshbach, who developed a theory 
and a classification scheme for resonance scatt ering phe- 
nomena in the context of nuclear physics ( Feshbachl . 
I1958L 1X9621) . We will follow here Fano's configura- 
tion interaction treatment of resonant scattering ijFanol . 
I1961D . which is common in atomic physics. A va- 
riety of treatments of the two-body physics of reso- 
nances in the context of ultracold Bose gases have been 

I n I — ) f 1 

given, for example. (IDuine and Stoof. 2004; Gora l et al. 



2004: Marcelis et al. 



Timmermans et al. 



.2004 
19991) . 



Raoult and Mica 12004 : 



We first consider the standard scattering picture away 
from any collision threshold defined by a two-channel 
Hamiltonian H. Assume that we can describe our system 
to a good approximation by two uncoupled "bare" chan- 
nels, as schematically shown in Fig. [TJ One is the open 
background scattering channel \bg) with scattering states 
\E) = (j)hg{R, E)\bg) labeled by their collision energy E. 
The other is the closed channel |c) supporting a bound 
state |C) = (j)c{R)\c) with eigenenergy Ec- The fimc- 
tions 4>c{R) and (phgiR, E) are the solutions to Eq. ([3]) for 
the background potential Vhg {R) and the closed channel 
potential Vc{R) respectively. Here (pdR) is normalized 
to unity. The scattering in the open channel is charac- 
terized by a background phase shift rihg{E). When the 
Hamiltonian coupling W{R) between the two channels is 
taken into account, then the two states become mixed, 
or "dressed" , by the interaction, and the scattering phase 
picks up a resonant part due to the bound state embed- 
ded in the scattering continuum. 



fl{E)=r]^g{E)+r],,,{E), 



(12) 



where r]rcs{E) takes on the standard Breit-Wigner form 
(Mott and Massey, 1965; Taylor, 1972): 



?/rcs(i^) 



tan 



E-Ec- 6E{Ec 



(13) 



The interaction W{R), which vanishes at large R, de- 
termines two key features of the resonance, namely, its 
width 

riE) = 27:\{C\WiR)\E)\\ (14) 
and its shift SE to a new position at Ec + SE{E), 

\{C\WiR)\E')\-' 



r 

SE{E) 

J —t 



E-E' 



-dE' 



(15) 



where V implies a principal part integral, which includes 
a sum over the contribution from any discrete bound 
states in the spectrum of the background channel. When 
the resonance energy is not near the channel threshold, 
it is normally an excellent approximation to take the 
width and shift as energy-independent constants, r(iJc) 
and SE{Ec), evaluated at the resonance energy Ec, as in 
Eq. (|13p. The resonance phase changes by w tt when E 
varies over a range on the order of F from below to above 
resonance. 

The essential difference between conventional and 
threshold resonance scattering is that if Ec is close to 
the open channel threshold a,t E — 0, the explicit 
ener gy-dependence of the wid t h and shift becorne cru- 
cial (iBohn and Juliennel . Il999l ; lJulienne and Gad . 120061 ; 
iMarcehs et al\ . \2004 ): 



■qrcsiE) = - tan ^ 



hnE) 



E- Ec- 5E{E) 



(16) 



The threshold laws for the s-wavc width and shift as k 
are 



1 



-T{E) ^ (fcabg)ro. 



Ec + 5E{E) -> Ei 



• 



(17) 
(18) 



where Fq and Eq are ^-independent constants. Since 
T(E) is positive definite, Fq has the same sign as Obg. 
Combining these limits with the background phase prop- 
erty, rftig{E) — fcobgj and for the sake of generality, 
adding a decay rate j/h ior the decay of the bound state 
into all available loss channels, gives in the limit of k 0, 



a ^ a — ib = ahe 



abgFo 



-£^0 + i(7/2) ■ 



(19) 



The unique role of scattering resonances in the ultra- 
cold domain comes from the ability to tune the threshold 
resonance position Eq through zero by varying either an 
external magnetic field with strength B or optical field 
with frequency v. 

Both magnetically and optically tunable resonances 
are treated by the same theoretical formalism given 



above, although the physical mechanisms determining 
the coupling and tuning are quite different. In the 
case of a magnetically tunable resonance, the channel 
can often be chosen so that 7 is zero or small enough 
to be ignored, whereas optical resonances are always ac- 
companied by decay processes 7 due to decay of the 
excited state. The resonance strength Tq is fixed for 
magnetic resonances, but ro(/) for optical resonances 
can be turned off and on by varying the laser inten- 
sity /. It may also be possible to gain some control 
over F n by using a combina tion of electric and magnetic 
fields (jiviarcelis adl2008[) . 

In the case of a magnetically tunable resonance, there 
is a difference dfi = /iatoms — fJ-c between the magnetic 
moment //atoms of the separated atoms and the magnetic 
moment fic of the bare bound state |C). Thus, the 
energy Ec of the state \C) relative to the channel energy 
of the separated atoms. 



= Sn{B - Be) 



(20) 



can be tuned by varying the magnetic field, and Ec is 
zero at a magnetic field equal to Be- Then, given that 
7 = 0, the scattering length takes on the simple form 
given in Eq. 



a{B) = abg 



A 



B-Bo' 



where 



£0 



and Bo = Be + SB 



(21) 



(22) 



are the width and the position of the singularity in 
the scattering length, shifted due to the interaction be- 
tween the closed and open channels by an amount SB = 
—SE/S^. Note that A has the same sign as (5/i/abg. Fig- 
ure m schematically illustrates the scattering length near 
the point of resonance Bq. 

The complex scattering length of an optically tunable 
resonance at laser frequency 1/ includes the coll i sional 
loss due to exci t ed stat e decay ( Bohn and Julienn3 . Il999t 
iFedichev er^ . Il996af) . 



d{v, I) = Obg 



abgro(/) 



h{iy-i^c-5iy{I)) + i{-f/2) 



(23) 



where the optically induced width Fo(/) and shift 5v{I) 
are linear in /, and Vc represents the frequency of the 
unshifted optical transition between the excited bound 
state and the collisional state of the two atoms &i E = Q. 

Whenever bound state decay is present, whether for 
magnetically or optically tunable resonances, Eq. (fT9|) 
shows that resonant control of the scattering length. 



a = abg - Arc 



7^0 



is accompanied by collisional loss given by 



b = -arc 



7 



(24) 



(25) 



The resonant length parameter 



To 



(26) 



i s useful for defining the st r ength of an op t ical re sonance 



(|Bohn and Juliennel . Il997t ICiurvlo et all, l2005^ or any 



other resonance with strong decay ( Hutsoni r 20071 ). Fig- 



ure [6] gives an example of such a resonance. The scat- 
tering length has its maximum variation of abg ± arcs at 
-Bo = ±7/2, where h = arcs- Resonances with arcs < l^-bgl 
only allow relatively small changes in scattering length, 
yet 6 remains large enough that they are typically ac- 
companied by large inelastic rate coefficients. On the 
other hand, if arcs S> |abg|j losses can be overcome by us- 
ing large detuning, since the change in scattering length 
is a — abg = —a-cesin / Eq) when |Bo| ^ 7, whereas 
6/|a-abg| = \\l/Eo\ < 1. 

The resonance length formalism is quite powerful. By 
introdu cing the idea of an energy-dependent scatter ing 
length (jBlume and Greend . I2OO2I : iBolda et alY l2002t ) it 



can be extended to Feshbach resonances in reduced di- 
mensional syste ms such as pancake or cigar shaped opti- 
cal lattice cells ( Naidon and Juliennd . l2006l ). 

While this discussion has concentrated on resonant 
scattering properties for B > 0, the near-threshold reso- 
nant properties of bound Feshbach molecules for energy 
B < are very i mportant aspects o f Feshbach physics; 
see Fig. H and (|Kohler et al\ . l2006f ) . In particular, as 
the bound state becomes more deeply bound, the closed 
channel character of the bound state increases and the 
binding energy Bb is no longer described by the univer- 
sal expression in Eq. The "dressed" or true molecular 
bound state of the system with energy — B;, is a mixture 
of closed and background channel components. 



\MR)) - VZ0,(i?)|c) + Xbg(i?)|&.9) , (27) 

where < Z < 1 represents the fraction of the 
eigenstate lil^biR)) in the closed channel component 
(|Duine and Stooi l2003l ). Unit normalization of \ipii(R)) 
ensures that J IxbgiR)]"^ dR = 1 — Since the variation 
of the energy — B^ with a parameter x of the Hamiltonian 
satisfies the Hellman-Feynmann theorem d{^Ei,)/dx ~ 
{ijb\dH/dx\i)b), it follows from Eq. ^ that 



Z = 



d{-Eb) SfXb 



dEc 



(28) 



Here Sfib = dEb/dB = A^atoms — fJ-b is the difference 
between the magnetic moment of the separated atoms 
and the magnetic moment fib of the "dressed" molecu- 
lar eigenstate. Since SfXb vanishes in the limit B ^ Bq, 
where Eb according to the universality condition in 
Eq. then Z vanishes in this limit also. Section llI.C.51 
develops more specific properties and conditions for Bf, 
and Z in this limit. 
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FIG. 6 Scattering length for an optically tunable Feshbach 
resonance as a function of laser tuning i/ — i/q. The lengths 
a and 6 are defined in Eqs. (|24p and (|25p . Here Obg — 5.29 
nm, Vo/h = 21 MHz at J = 500 W/cm^, a^os = 5.47 nm, 
and 7//t = 20 MHz. Numerical values for the strength and 
spontaneous linewidth of the resonance a re typ ical for *^Rb 
and are taken from Fig. 1 of (jTheis et all\2004 ). 



B. Basic molecular physics 



Most atoms that can be trapped at ultracold tem- 
peratures have ground S-states with zero electronic or- 
bital angular momentum (L ~ 0), as for alkali- metal or 
alkaline-earth-metal atoms. The collision between two 
atoms is controlled by the electronic Born-Oppenheimer 
interaction potential(s) between them. All potentials are 
isotropic for the interaction of two S'-state atoms. We 
restrict our discussion of molecular physics to such cases. 
Figure [7| shows as an example the and '^S+ poten- 
tials for two ground state Li atoms, which are anal- 
ogous to the similar potentials for the H2 molecule or 
other alkali-metal atoms. The superscripts 1 and 3 refer 
to singlet and triplet coupling of the spins of the un- 
paired electrons from each atom, i.e., the total electron 
spin S = Si -I-S2 has quantum numbers S = and 1. The 
I] refers to zero projection of electronic angular momen- 
tum on the interatomic axis for the 5-state atoms, and 
g/u refers to gerade/ungerade electronic inversion sym- 
metry with respect to the center of mass of the molecule. 
The g/u symmetry is absent when the two atoms are not 
of the same species. 

The Born-Oppenhcinier potentials are often available 
from ab initio or semi-empirical sources. When R is 
sufficiently small, typically less than i?ox ~ 20 ap ~ 1 
nm for alkali-metal atoms, electron exchange and chem- 
ical bonding effects determine the shape of the poten- 
tials. For R ^ i?ex, the potentials are determined by the 
long-range dispersion interaction represented by a sum of 
second-order multipolar interaction terms. 
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FIG. 7 Molecular potentials V{R)/h versus R of the two 
electronic states of Li2 that correlate with two separated '^S 
atoms. The inset shows an expanded view of the long-range 
s-wave potentials of ^Li at B = 0, indicating the five hyperfine 
states of the separated atoms (see Fig. [S} for which the to- 
tal angular momentum has projection M = 0. The inset also 
shows the last two nearly degenerate bound states (unresolved 
on the figure) of the ^Li2 molecule from a coupled channels 
calculation. It is a good approximations to label these nearly 
degenerate levels as the / = and 2 components of the total 
nuclear spin I = Ii + 12 of the last v = 38 vibrational level of 
the ^E^ potential. 



1. Van der Waals bound states and scattering 

Many aspects of ultracold neutral atom interactions 
and of Feshbach resonances, in particular, can be un- 
derstood qualitatively and even quantitatively from 
the scattering and bound state properties of the long 
range van der Waals potential. The properties of this 
potential relevant for ultracold ph otoassociation spec - 



troscopy have been reviewed in (iJones et al. 



Its a i ialytic properties are discussed bv (IGao 



2OOI. 



1998b . 



20001: iGribakin and Flambauml . Il993t iMott and Massev 
19651) 



In the case of S'-state atoms, the lead term in the long- 
range part of all Born-Oppenheimer potentials for a given 
atom pair has the same van der Waals potential charac- 
terized by a single Cg coefficient for the pair. Conse- 
quently, all qiq2 spin combinations have the long-range 
potential 



ViiR) 



9i 

R^ 



2^ i?2 



(29) 



A straightforward consideration of the units in Eq. ([29|) 
suggests that it is useful to define length and energy scales 



R 



_ 1 ^ 2/iC6 

vdw-2l f2 



1/4 



and E'v 



1 



(30) 
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( Gribakin and Flambauml . Il993f) defined an alternative 
van der Waals length scale which they called the mean 
scattering length: 



a = 47r/r(l/4)2 R^^„ = 0.955978 ... R 



vdw ; 



(31) 



where T{x) is the Gamma function. A corresponding en- 
ergy scale is E ^ ^1^(2^^) = 1.09422 ... i^vdw The 
parameter a occurs frequently in formulas based on the 
van der Waals potential. Table U gives the values of 
i?vdw and Evdw for several c ases. Values of Ce for 
other systems are tabulated in (iDerevianko et ~ai. 1. 119991: 
iPorsev and Dereviankol . 120061: iTang et al\ . \l%T^ 

The van der Waals energy and length scales permit a 
simple physical interpretation ( Juhenne and Miea . ll989t ). 
A key property for ultracold collisions is that Ce / be- 
comes large compared to the collision energy E when 
R < Rvdw Thus, the wave function for any partial wave 
oscillates rapidly with R when R < Rydw , since the lo- 
cal momentum hk{R) = ^J2^{E — V{R)) becomes large 
compared to the asymptotic hk. On the other hand, 
when R > i?vdw , the wavefunction approaches its asymp- 
totic form with oscillations on the scale determined by 
the long de Broglie wavelength of the ultracold collision. 
The energy scale i?vdw determines the nature of the con- 
nection between the long- and short-range forms of the 
wavefunction. The de Broglie wavelength A = 27r(i?vdw) 
for E = i?vdw When E < E^dw so that A > i?vdw, a 
WKB connection can not be made near i?vdw between the 
asymptoti c s-wave and the sh ort-range wavefunction (see 
Fig. 15 of ( Jones et aZ.I . l2006l )). Consequently, the quan- 
tum properties of the collision are manifest for E < E^dw ■ 

The van der Waals length also characterizes the ex- 
tent of vibrational motion for near-threshold bound state. 
The outer turning point for classical motion for all low £ 
bound states is on the order of Rvdw The wave function 
for £ = oscillates rapidly for R < i?vdw and decays ex- 
ponentially as e^^*"^ for R S> Rvdw, where h^kl/{2^) is 
the binding energy. The only case where the wave func- 
tion extends far beyond i?vdw is that of the last s wave 
bound state for the case of the universal halo molecule, 
where a ^ Rvdw', sec Sees. III.Al and lV.B.21 

The van der Waals potential determines the interaction 
over a wide zone between i?vdw and the much smaller i?ox 
where chemical forces become important. Thus, near- 
threshold bound and scattering state properties are de- 
termined to a large extent by the long-range van der 
Waals potential. The effect of short-range is then con- 
tained within the phase of the wave function, or equiv- 
alently, the log de rivative ( Moerdiik and Verhaail . Il994l : 
IVoeels aLl . l200Qll . More precisely, for any Rz satisfying 
Rex < i?z <C i?vdw so that k{Rz) ^ k, the wave function 
phase is nearly independent of E and almost the same 
for all near-threshold bound or scattering states. In fact, 
the phase is nearly independent of partial wave i as well, 
since the centrifugal potential is typically small compared 
to the van der Waals potential for such an R^ . Using this 
phase as a boundary condition for propagating the wave- 
function to large R in the asymptotic domain determines 



the energy-dependent scattering phase r]i{E) and bound 
state energies. In fact the phase of the wave function in 
the zone i?ox < Rz ^ Rvdw i s uniqu ely related to the s- 
wave scattering length ( Gaol . [l998a[ ). Consequently, to a 
good approximation the near-threshold bound states and 
scattering properties for all low partial waves are deter- 
mined by the s-wave s cattering le ngth, the Cg constant, 
and the reduced mass ((Gaol . 1200 Ih . 

(iGaol . I2000h has worked out the energies Enj of the 
bound states of all partial waves for a van der Waals 
potential as a function of the s-wave scattering length, 
where n = —1,-2,... is the vibrational quantum num- 
ber and £ is the rotational quantum number of the bound 
state. He showed that the energies of weakly bound 
states have a A£ = 4 periodicity. Figure |H| shows bound 
state energies as function of £ for two values of a. In 
the left panel a = ±oo so there is a s-wave bound state 
with E = 0. The figure shows that for £ = 4 the r e also 
is bound state with E/Evdw — 0. In fact ( Gaol . |200G| ) 
showed for £ = 8,12,... there will be a bound state at 
zero energy as well. The right panel in Fig. |5| shows that 
when a = a there is a bound state at zero energy for 
£ ^ 2. There also will be a bound state at zero energy 
for £ = 6, 10, . . .. 

Figure |5| can also be used to define the concept of "en- 
ergy bins" in which, regardless of the value of a, there 
must be a bound state. Bins arc most easily defined by 
starting from a case with a bound state at zero bind- 
ing energy. By changing the short-range log derivative 
its binding energy can be increased, or its energy low- 
ered, and at some point the binding energy is so large 
that a new bound state appears at zero binding energy. 
This is exactly the situation depicted in Fig. [5^ for s 
and g waves. In other words, for s waves there must be 
a. n — —1 bound state between — 39.5£'vdw and Oi^vdw, 
while for g waves there must be a n = — 1 bound state be- 
tween — 191 -Evdw and Evdw- The n = —2 s-wave bound 
state appears between — 272iJvdw and — 39.5i?vdw Fig- 
ure IHJd can similarly be used to define the bins for other 
waves. 

When the scattering length is large compared to a and 
positive, a simple expression for the van der Waals correc- 
tion to the binding energy o f the last s-wave bound state 
can be worked out (jGribakin and Flambauml . 1 1993[ ): 



E^ 



1.0 



2/-t(a — ay 



(32) 



The universal formula in Eq. ^ only applies in the limit 
that a 3> a and l-E-i^ol ^ -Evdw CGao, 2 00^ has worked 
out higher order corrections to the binding energy due to 
the van der Waals potential, which can be recast as 



9ia 



520 



2^{a — a)^ 



(a - ay 



(33) 



Here gi = r(l/4)V67r2-2 = 0.9179..., = (5/4)gJ-2 
—0.9468... are constants. 
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TABLE I Characteristic van der Waals scales i?vdw and -Bvdw 
for several atomic species. (1 amu — 1/12 mass of a ^'^C atom, 
1 au= 1 EhaQ where Eh is a hartree and 1 ao= 0.0529177... 
nm) 



Species 


mass 


C6 


-^vdw 


-tivdw/ Kb 






(amu) 


(an) 


(ao) 


(mK) 


(MHz) 


«Li 


6.0151223 


1393.39" 


31.26 


29.47 


614.1 


2^Na 


22.9897680 


1556'' 


44.93 


3.732 


77.77 




39.9639987 


3897'' 


64.90 


1.029 


21.44 


«Ca 


39.962591 


2221= 


56.39 


1.363 


28.40 




86.909187 


4698'' 


82.58 


0.2922 


6.089 


»«Sr 


87.905616 


3170'= 


75.06 


0.3497 


7.287 


Cs 


132.905429 


6860"= 


101.0 


0.1279 


2.666 



a. (lYan ef ami996f) b. (iDerevianko et a/.Ul999D 
c. (iPorsev and DereviankoL l2002ll 
d. (|van Kempen et all 120021 ') e. (|Chin et all [2004bl ) 




FIG. 8 Bound-state energies of the last vibrational levels of 
two atoms interacting via a van der Waals potential as a func- 
tion of partial wave £. The zero of energy is at two free atoms 
with zero collision energy ("at threshold"). The lowest par- 
tial waves are shown. Panel a) shows the bound state struc- 
ture up to ^ = 8 when the scattering length of the collid- 
ing atoms is infinite or equivalently that there is an s-wave 
bound state with zero binding energy. Panel b) shows the 
bound state structure up to = 6 when the scattering length 
is a = 0.956 . . . -Rvdw = a, or equivalently there is a d-wave 
bound state with zero binding energy. The length i^yd w and 
energ y -Evdw are defined in the text. Adapted from (|Gacil . 
l2000l i 



Similarly, the effective range of the potential in Eq. ([5|) 
is also determined from the van der Waals potential , 
given the s-wave scattering length ( Flambaum et all . 
ll999l : [Gaol . ri998aD : 



r(i/4)4 

^■o = — a 

DTT 



1-2^ 
a 



(34) 



where r(l/4)'*/(67r) w 2.9179. When a > a, this simpli- 
fies to To = 2.9179a. Note that tq diverges as a ~> 0. 
The energy levels of the van der Waals potential are 



not exact due to the slight influence from the actual 
short-range potential and extremely long-range retarda- 
tion corrections. They are nevertheless relatively accu- 
rate guides to the expected energy spectrum for real 
molecules. For example when the scattering length is 
slightly larger than a, which corresponds to Fig. [5)3 
with all bound states shifted to slightly more posi- 
tive energies, the d-wave bound state becomes a shape 
resonance, that is, a decaying quasibound state with 
E > trapped behind the d-wave centrifugal barrier. 
For ^^Na and *^Rb the experimentally observed scat- 
tering length is 10%-20% larger than a and, indeed, 
in both cases a d-wave shape res onance has been ob- 



served under various circumstances (IBoesten et al.. 1997 



Buggle et ai\ . |2004 ISamuelis et'd\ . l2000HThomas al 



20041 ) . Similarly, a p-wave shape resona nce occurs when 
a is s lightly larger t han 2d, as for "^ ^ K (iDeMarco et all . 
I1999D and ^^^Yb (iKitaeawa aZ.! . l2008^ . In addi- 
tion ( Kitagawa et aZ.I .~ 20081 ) show how the scattering 
length and binding energies of the last few bound states 
for the single potential of the Yb -I- Yb interaction are 
related as the reduced mass is changed by using different 
isotopic combinations of Yb atoms. The scattering length 
and binding energies can be "tuned" over a wide range by 
choosing different pairs of atoms among the seven stable 
isotopes of Yb. 

The 6£ = A characteristic of van der Waals poten- 
tials also has practical consequences for ultra-cold scat- 
tering. For ^^Rb the scattering length has been found 
to be large compared to R ydw and a q-wave shap e res- 
onance has been observed ( Boesten et all Il996a|) . For 
^^■^Cs the scattering length is large compared to i?vdw, 
and numerous (7- wave bound states with bin ding energies 
much smaller than -Evdw were observed by ( Chin et all 
l2004bD at low magnetic field. In fact, some of these bound 
states appear as magnetic Feshbach resonances in the col- 
lision of two Cs atoms. Recently, a weakly- bound I = 8 
or /-w ave state has been observed as well ( Mark et~al\ . 
l2007al) . 



2. Entrance- and closed-channel dominated resonances: 
Resonance strength 

The van der Waals theory is very useful for characteriz- 
ing and classifying the basic properties of the resonances 
discussed in Section ITl. A. 3 1 by expressing lengths in units 
of a and energies in units of E] see Eq. pip. The numer- 
ator of the resonant term in Eq. p9p defines a resonance 
strength parameter to be abgFo, where Fq — SfiA; see 
Eq. (P^ . It is helpful to define a dimensionless resonance 
strength parameter s^es to be 



Fo _ flbg SfiA 



E 



E 



(35) 



where ?'bg ~ flbg/d is the dimensionless background scat- 
tering length. The sign of s^cs is always positive. The 
resonance phase in Eq. p6p is determined by the tunable 
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resonance position and the resonance width and shift. In 
the hmit E ^ 0, both the width 



(36) 



and the shift (jGoral et aZ.l . l2004lJuhenne and Gaoll2006D 



SE = 



1-rbg 



1 + (1 - rbg)- 



-{ESrcs) , 



(37) 



are proportional to Es^cs- Section ril.B.SI describes widths 
and shifts for some typical resonances. Sections III.C.3| 
III. C. 41 and III. C. 51 give additional analytic properties of 
threshold scattering and bound states associated with 
Feshbach resonances and show how Eq. ((37|) can be de- 
rived. 

The strength parameter Sros allows us to classify Fesh- 
bach resonances into two lim iting cases. ( Kohler et all . 
120061: IStoU and Kohled . [200l used r] = 1/sres to do this. 
When Sres ^ li the resonance is called an entrance chan- 
nel dominated resonance. Here, the near-threshold scat- 
tering and bound states have the spin character of the 
entrance channel for detuning Eq over a large fraction of 
the width Tq and thus for B — Bq over a large fraction of 
the resonance width A. In this regime, the resonance can 
be well modeled by the B-dependent scattering length of 
Eq. ([1]). The bound state is universal with Z <C 1 (see 
Eq. (j27p ) over this large detuning range and with a bind- 
ing energy well-approximated by Eq. Resonances of 
this type have the largest resonance width A and are 
conventionally called "broad resonances" . 

Resonances with Sjos 1 arc called closed channel 
dominated resonances. Here, the near-threshold scatter- 
ing and bound states have the spin character of the en- 
trance channel only over a small fraction of the width 
Tq near Eq = and thus over a small fraction of the 
resonance width A near B ~ Bq. A universal bound 
state only exists over this small detuning range. Thus, 
the closed channel fraction Z is only small near B ~ Bq 
and is near unity over a wide detuning range away from 
B = Bq. Such resonances need to be modeled by a cou- 
pled channels description. Resonances of this type often 
have a small width A and are conventionally called "nar- 
row resonances." 

It should be emphasized that the conventional use of 
"broad" or "narrow" resonances referring to those that 
can or cannot be modelled by single channel model is not 
rigorously defined. Exceptions exist where resonances 
with apparently broad widths are actually closed channel 
dominated. The terms introduced here, entrance and 
closed channel dominance, better reflect the nature of the 
near threshold states over a detuning range on the order 
of the width A and can be unambiguously assigned to a 
resonance by evaluating Sros- 

Section III. B. 51 illustrates the differences between en- 
trance and closed channel dominated resonances by giv- 
ing specific examples of such resonances. Section III. C. 51 
develops a simple model for the bound states for any 



type of resonance, and shows that the norm Z of the 
closed channel part vanishes in the limit that Eq 
near the point of resona nce Bq, even for c losed channel 
dominated resonances. (jSzvmanska et all [20051 ) discuss 
in detail the implication of the distinction between open 
and closed channel dominance for the modeling of many- 
body systems, a topic that is beyond the scope of this 
review. 

Figure [49] in the Appendix illustrates the wide range of 
resonance strengths Sics and widths A observed for var- 
ious alkali atom resonances. Broad resonances with A 
larger than 1 G tend to have Sres > 1 and thus be en- 
trance channel dominated ones. Narrow resonances with 
A smaller than ~ 1 G tend to have s^cs < 1 and thus be 
closed channel dominated ones. A notable exceptions is 
the ^Li 737 G resonance with Sros < 1 that is very broad 
yet tends towards being closed channel dominated; see 

Section ULRII 

Equation ([55)1 allows us to address the question 
whether a sharp resonance feature appears at small but 
finite collision energy above threshold. A condition for 
having a sharp resonance is that the width lT{E) should 
be smaller than the collision energy E. It is convenient 
to rewrite Eq. p6)) as ^T{E) = {s,:cB/{ka))E. For an en- 
trance channel dominated resonance with s^cs S> 1 and 
ka < 1 OT E < E, it follows that ^r{E) > E . Thus, 
there can be no sharp resonance features evident in the 
above-threshold phase i]{E, B) of an entrance channel 
dominated resonance when E < E. A sharp resonance 
featu re can only appear when E ^ E. |Nvgaard et all . 
1200611 illustrate this case for a resonance involving 
atoms. On the other hand, for a closed channel domi- 
nated resonance with Si-es <C 1 a sharp resonance feature 
in T]{E,B) with ^r{E) < E can appear immediately 
above threshold. 



3. Coupled channels picture of molecular interactions 

While many insights can be gained from the proper- 
ties of the long range van der Waals potential, actual 
calculations require taking into account the full molec- 
ular Hamiltonian, including not only the full range of 
the Born-Oppenheimer potentials but also the various 
spin-dependent couplings among them. In general, the 
potential should be viewed as a spin-dependent potential 
matrix, the elements of which account for the interac- 
tion among the various spin states of the atoms. The 
wave function for atoms prepared in channel a can be 
written as a coupled channels expansion in th e separated 
atom spin b asis describ ed in Section II. A. II (iGao . 1996 : 
Hutson et al. 200 8: Mies et all. l2000l : 
1965l:IStoof 6^011 . 119881) 1 



Mott and Massev 



|Va(i?, E)) = \P)MR, E)/R . (38) 


The allowed states in this expansion are those that 
have the same projection quantum number M = mi + 



14 



7712 +1^11 (see Section fll. A. 1[) . In addition parity conser- 
vation implies that all partial waves in the expansion 
are even if £a is even and odd if £a is odd. 

Subsituting Eq. ([38| in the Schrodinger equation gives 
the coupled channels equations for the vector 0(i?, E): 



(fcj){R,E) 
2^ dR2 



-V{R)4){E,R) ^ E(f>{E,R). (39) 



The solutions to these coupled equations give the bound 
and scattering states of the interacting atoms. The po- 
tential matrix V(i?) gives the matrix elements of the 
Hamiltonian between the channel basis sets. It takes on 
the following form in the asymptotic spin basis: 



Er, 



flh^ji^ + 1) 



Safl + VintMR) (40) 



The interaction matrix Vint(^) contains the electronic 
Born-Oppcnheimer potentials discussed in Section III.BI 
and the relativistic electron spin-dependent interactions: 



Vint(i?) = Vel(i?)+V,,(i?) 



(41) 



For collisions of S state atoms the term Vci(i?) represents 
the strong isotropic electronic interaction that is diagonal 
in £ and mi but off-diagonal in the atomic channel quan- 
tum numbers (71(72- The diagonal elements Vd^aa vary at 
long range as the van dcr Waals potential (see Fig. [7]). It 
has normally been unnecessary to include small retarda- 
tion or nonadiabatic corrections to long range molecular 
potentials in order to fit experimental data on ground 
state c ollisions within their ex perimental error; see for ex- 
ample (jKitagawa et "^■ 120081 ). The off-diagonal elements 
KjI.q/S, where /3 7^ a, decrease exponentially at large R 
as the exchange potential and become small compared 
to the atomic hyperfine splitting for R > i?ex- The Vd 
coupling is responsible for elastic scattering and inelas- 
tic spin-exchange collisions and gives rise to the largest 
resonance strengths. 

The term Vss(7') in Eq. (|41]) represents weak rel- 
ativistic spin-dependent interact ions. These include 
the spin-sp i n dip olc interaction ( Moerdiik et all Il995l : 



IStoof ei al-l. 119881) and the second -order spin-orbit inter- 
action ( Kotochigova et all . l2000[ l. important for heavy 



atoms (jLeo et aLl . 120001 ). The two contributions are both 
anisotropic and are off-diagonal in both qiq2 and i. Thus, 
Vss(7') couples different partial waves. At long range 
Vss(?') is proportional to /R^, where a = 1/137.0426 
is the fine structure constant. This anisotropic potential 
only contributes diagonal terms for partial waves £ > 1, 
and does not contribute to the potential Vint,QQ(^) when 
a represents an s-wave channel. The 'Vsa{i') coupling 
is responsible for weak inelastic relaxation and normally 
gives rise to small resonance strengths. 

The Born-Oppenheimer potentials are normally never 
known with sufficient accuracy to permit accurate calcu- 
lations of threshold scattering properties. Consequently, 
it is usually necessary to vary the short range potentials 



over some range of i? < i?ox to calibrate theoretical mod- 
els so they reproduce measured threshold bound state or 
scattering data. In some cases the van der Waals co- 
efficients are accurately known, whereas in other cases 
they need to be varied to fit the data as well. Once this 
is done, coupled channels theoretical models typically 
arc robust and predictive of near-threshold collision and 
bound state properties. Some examples of high quality 
theoretical models b ased on fitting Feshbach resonance 



data are given 
for 23Na, ilChin et al. 



van Abeelen and Verhaai]. 1999al ) 



"^Cs , (iMarte et all |2002|) for ^ ^ Rb. (iBartenstein et ai\ . 



2OOOI) and (iLeo et all I2OOOD for 



I2OO5D for ^Li (IWerner gLl. I2005D f or '^^Cr. and 



(iFerlaino et oil . 120061: iPashov al\ . l2007h for ^OK^^Rb. 



Information on other models can be found in the refer- 
ences listed in Sections IIII.BI and IIII.Cl 



4. Classification and molecular physics of Feshbach resonances 

The previous Sections have laid the groundwork for 
classifying and understanding the properties of Feshbach 
resonance states in ultracold collisions of ground 5-state 
atoms. This classification can be made according to the 
quantum numbers {qiq2£M\qc£c\ , where {qiq2£M} char- 
acterize the entrance channel (see Section III. A. II) and 
{qc£c} characterize the "bare" closed channel bound state 
that gives rise to the resonance. Such a bound state 
has the same M as the entrance channel. Some possible 
choices for quantum numbers comprising the composite 
qc are given below. 

It is important to note that £c need not be the same as 
the entrance channel partial wave £. Parity conservation 
ensures that |^ — is even. In the case of two L = 
atoms the Vd term is isotropic and only gives rise to 
nonzero matrix elements when £c = £. On the other 
hand, £c can be different from £ for the anisotropic Vgs 
term. We are primarily concerned with entrance channel 
s-waves, although some resonances with p- wave (e.g., ^Li, 
^"K, ^^"^Cs) or c^-wave (^^Cr) entrance channels are known 
in the fiK domain. 

We find it is convenient to designate resonances ac- 
cording to the value of the closed channel bound state 
quantum number £c, as shown in Table [III If £c is even 
(odd), we assume an s- (p-) wave entrance channel un- 
less otherwise stated. The strongest resonances with the 
largest widths A are s-wave resonances with £ = £c = 0, 
and are due to the Vd term in the Hamiltonian. A num- 
ber of weak resonances with small A are known where 
the s-wave entrance channel is coupled through the Vgg 
term to bound states with even £c such as 2 or 4. Follow- 
ing Table im the latter are designated as d-wave or g-wave 
resonances, respectiv ely. For example, d -wave resonances 
are known for ^ ^Rb (iMarte et d. I I2OO2D and (7- wave res- 



onances for Cs ( Chin et all . l2004bl ) . For s-wave entrance 



channels the g-wave resonances are only possible due to 
second-order coupling in Vgs . Entrance channel p- waves 
can be coupled to resonant bound states of odd = 1 or 
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3, although the latter would tend to be quite weak and 
rarely observed. 



TABLE II Classification of magnetic Feshbach resonances in 
collisions of ultracold atoms. The type of the resonance is 
labeled by the partial wave 1^ of the closed channel bound 
state rather than the entrance channel partial wave I. Almost 
all cases known experimentally have = or 1. Note that 
identical bosons (fermions) in identical spin states can only 
interact with even (odd) partial waves. All other cases permit 
both even and odd partial waves. 



Type 


I 




s-wave resonance 


0, 2 ... 





p-wave resonance 


1, 3 ... 


1 


d-wave resonance 


0, 2 ... 


2 


/-wave resonance 


1, 3 ... 


3 


g-wavo resonance 


0, 2 ... 


4 



The long-range potential Vci.apijf) is diagonal for the 
interaction of two ground state alkali metal atoms for 
all combinations q\q2 of Zeeman sublevels, and all chan- 
nels have the same van der Waals coefficient Cg. Conse- 
quently each channel will have a spectrum of vibrational 
and rotational levels for a van der Waals potential as de- 
scribed in Fig. [HI The n = — 1, — 2, . . . levels associated 
with closed spin channels a' can become scattering res- 
onances for entrance channel a if they exist near energy 
Ea. The value of n can be one of the values compris- 
ing the set of approximate quantum numbers Qc- Some 
examples of the use of n in r esonance classification are 
given bv (iMarte et all . |2002[ ) for 87Rb (see Fig. ^ or 
(jKohler et all . |2006[) for *^Rb. The vibrational quantum 
number can be either n, counting down from the top, or 
V, counting up from the bottom of the well. 

The approximate spin quantum numbers in qc are de- 
termined by whatever set of quantum numbers blocks the 
Hamiltonian matrix into nearly diagonal parts. This will 
depend on the nature of the coupling among the various 
angular momenta of the problem so that no unique gen- 
eral scheme can be given. For alkali-metal dimers the 
spacing between vibrational levels, which is on the order 
of tens of Eydw as seen from Fig. [51 must be compared 
to the spacing between the channel energies Ea. For ex- 
ample, in a light molecule like ^Li2 or ^^Na2 they are 
large compared to the atomic hyperfine splitting Em = 



\E 



1+1/2 - 



E 



1-1/2 



and Zeeman interactions. In this case. 



the vibrational levels arc to a good approximation clas- 
sified according to the electronic spin coupling, S* = or 
1 of the respective ^S]+ and Born-Oppenheimer po- 
tentials, with additional class ification according to their 
nuclear spin substructu re . ( van Abeelen and Verhaail 
Il999al : iLaue et all . l2002t ) give an example of suc h a clas- 
sification for ^^Na2, and ( Simonucci et aZ.l . l2005l ) give an 
example for ^Li2. 

In contrast to light species, heavy species like Rb2 or 
Cs2 have vibrational spacings that are smaller than Em, 
so that near-threshold bound states of of the and 



potentials are strongly mixed by the hyperfine in- 
teraction. The near-threshold molecular states do not 
correspond to either 5* = or 1, but often can be charac- 
terized by the approximate quantum number /c, where 
f = fi -|- f2 . As with the /i or /2 atomic quantum num- 
bers, fc is not a good quantum number at large B but can 
be used as a label according to the low field sta t e with 
which it adiabatically correlates. (Marte et Ml, I2002D 
give examples of such r esona nce classification for ^^Rb, 
and (|Chin et all l2004b[ ) and (iKohler et all. I2006D d o so 
for "3Cs and ^'^Rb respectively. (|Hutson et aLl . l2008[ ) de- 
scribe improved computational methods for calculating 
the coupling between bound state levels and characterize 
a number of experi mentally observed avoided crossings 
([Mark et all l2007a[ ) between Cs2 levels having different 
approximate quantum numbers. 



TABLE III Separated atom channel labels for the five s-wave 
M — channels of ®Li. The (/1/2) quantum numbers are only 
exact at B = 0. 



Q 


ih h) 




ah 


{--) 


' 2 ' 2 


ad 


(--) 


' 2 ' 2 


he 


\2 2I 


2 ' ' 2 


cf 


\2 2! 


3 ,3 
2 ' 2 


de 


\2 2l 


^2- 2 



5. Some examples of resonance properties 

We will use the fermionic species ^Li to illustrate some 
basic features of Feshbach resonances. Figure[5|shows the 
atomic Zeeman levels. The inset to Fig. [7| shows the 5 
channels and the potentials Vaa {R) at long range needed 
to describe the s-wave collision of an gi = a atom with a 
q2 = b atom. These 5 channels summarized in Table IIIll 
have the same van der Waals Cg coefficient and the same 
projection M = 0. Due to the light mass, the last bound 
state of the van der Waals potential must lie in a "bin" 
that is 39.5Svdw//i = 24.3 GHz deep. 

Figure [7| shows that the last two M = coupled chan- 
nels s-wave bound states for B ~ have the character 
of the 71 = — 1 or w = 38 level of the potential. 
They have a binding energy of «1.38 GHz relative to the 
separated atom energy Eat, associated with the positive 
scattering length of a = 45 . 17 ag of the S = singlet 
potential ( Bartenstein et al] . l2005l ). The two levels have 
total nuclear spin / = or 2 and projection mj ~ 0, 
where I = Ii -I- I2. The next bound states below thresh- 
old are three M ~ spin components of the v = 9 level 
of the ^5]+ potential, far below threshold with binding 
energies «24 GHz near the bottom of the n = — 1 "bin" . 
These deeply bound levels are associated with the large 
negativ e scattering length of --2140 ap for the S = 1 p o- 
tential ([Abraham ad [l997t [Bartenstein et all [2005[ ). 
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FIG. 9 Last coupled channels bound states of the ®Li2 dimer 
with M = 0. The arrows indicate the locations of the 543 
G and 834 G Feshbach resonances, where the binding energy 
of a threshold bound state equals 0. While the low B field 
7 = 2 = 38) level retains its spin character as it crosses 

threshold near 543 G, the 1 = level mixes with the entrance 
channel and switches near « 550 G to a bound level with 
ah spin character, eventually disappearing as a bound state 
when it crosses threshold at 834 G. 



Figure [5] shows how the channel energies and the ener- 
gies of the last two M = s-wave bound states of the ^Li 
mole cule vary with magnetic field 



(Simonucci et al 



I2005D give a detailed description of the molecular physics 
of these multichannel bound states. At high B field. Eat 
varies linearly with B with a slope of nearly dEat / dB = 
—2fiB, where /xb is a Bohr magneton. Since both bound 
states have 5 = character near B = 0, their magnetic 
moment vanishes, i.e, dEc/dB = near 5 = 0. The 
1 = 2 state crosses Eat near B = 543 G, where it inter- 
acts weakly with the ab entrance channel and makes a 
very narrow Feshbach resonance, shown in Fig. 1111 On 
the other hand, the energy of the strongly interacting 
1 = bound state changes dramatically above about 540 
G and becomes nearly parallel to the energy of the ab 
entrance channel. This state switches to 5 = 1 character 
near the 540 G crossing region, and transforms into the 
V = 10 level of the '^Sj potential at higher B. This level 
becomes a very weakly bound "universal" halo state of 
dominantly entrance channel character above around 650 
G, and does not disappear until it reache s the Eab thresh- 
old near 834 G (|Bartenstein et aZ.l . l2005f) . where it makes 
a very broad Feshbach resonance, shown in Fig. [TOl 

Figure [TOlshows the near-threshold bound and scatter- 
ing state properties of the ®Li ab channel between 500 G 
and 1200 G, while Fig. [TT] shows an expanded view of the 
narrow resonance near 543 G. The energy range is typ- 
ical of the ultracold domain. The two figures illustrate 
two extremes of resonance behavior. The 834 G reso- 
nance is strongly open channel dominated with Sj-cs — 59 
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FIG. 10 Near-threshold bound and scattering state proper- 
ties of the ^Li ab channel. The upper panel shows the coupled 
channels scattering length versus magnetic field B using the 
model of (|Bartenstein et all |2005| ) . The double-headed arrow 
indicates the point of singularity -Bo for the broad resonance 
near 834 G, which is an entrance channel dominated resonance 
with width A = 300 G and Srcs = 59 (see Sections III. A. 31 and 
III.B.2|I . There also is a narrow resonance with a singularity 
near 543 G. The lower panel shows for E < the energy 
of the bound state (solid line) that merges with the contin- 
uum at Bq. The zero of energy at each B is the ab channel 
energy Eab{B). An energy of E/h = 0.4 MHz is equivalent 
to E/ks = 19 /iK. The universal bound state energy from 
Eq. ^ is indistinguishable on the scale of this graph from 
the coupled channels bound state energy. The nearly verti- 
cal dotted line shows the energy of the "bare" bound state 
Ec{B) of ^T.+ {v = 38,1 = 0) character that crosses thresh- 
old at Be near 540 G. The shaded contour plot for E > 
shows sin^ ri{E, B). The broad light-colored region near the 
point of resonance indicates the region where sin^ 77 « 1 and 
the cross section is near its maximum value limited by the 
unitarity property of the S'-matrix. Since Sres 1 the width 
r(£') in Eq. (|36|l is larger than E in the near-threshold region 
so that there is no above-threshold "resonance" feature in the 
collision cross section versus E. 



(see Eq. ([55]) ) and is well represented by a universal halo 
bound state of entrance channel character over a large 
fraction of its width A. The 543.2 G resonance is strongly 
closed channel dominated with Sjos = 0.001. It exhibits 
open channel character and universal behavior only over 
a negligible detuning range spanning at most a few /iG 
when B is tuned near Bq. 

It is instructive to examine the wave functions for the 
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FIG. 11 Expanded view near 543 G of Fig. \T0\ The up- 
per panel shows the coupled channels scattering length versus 
magnetic field strength B, where the double-headed arrow in- 
dicates the calculated point of singularity Bo for the narrow 
resonance at 543.18 G with a width of A = 0.10 G in ex- 
cel lent agreement with t he measured resonance at 543.26(10) 
G (jStrecker et aZ.I . l2003t ). This is a closed channel dominated 
resonance with Srcs = 0.001. The lower panel shows for E < 
the energy of the coupled channels bound state that merges 
with the continuum at Bq. The dashed line shows the uni- 
versal bound state energy from Eq. ([2|. Universality does 
not apply for detunings over most of the width of the res- 
onance but is only applicable in an extremely narrow range 
very close to Bo. The shaded contour plot of sin'^ vi^y B) for 
-E > shows a very narrow and sharp resonance emerging 
above threshold with a width T{E) <^ E, with a linear vari- 
ation of position with B, and with a very small domain of 
unitarity of the S-matrix. 



coupled channels bound states with the same binding 
energy near each resonance. For example, the binding 
energies in Figs. [10] and [TT] are « 200 kHz near 700 G 
and 543.1 G respectively. We calculate that the projec- 
tion Z on the closed channel components, Z = 1 — Ua 
are 0.002 and 0.98 for these respective cases, where 
Ua = /q°° \(f>ab.ab{R)\^dR IS thc norm of the entrance 
channel component (j)ab.ab{R) of the bound state from 
the coupled channels expansion in Eq. (|38p . The small 
projection Z at 70 G is in good agreem ent with the 
value measured by ( Partridge et al\ . |2005[ ): see Fig. [SO] 
These projections for levels with the same near-threshold 
binding energy illustrate the very different character of 
entrance and closed channel dominated resonances. 



The width A itself does not determine whether a res- 
onance is entrance or closed channel dominated. Rather, 
it is necessary to apply thc criterion in Eq. ([35]) . A 
good example of this is provided by the bosonic ''Li sys- 
tem, which has a very broad resonan ce in the aa channel 
near 737 G with a widt h of 192 G (I Junker et al\ . 120081 : 
iKhavkovich et al\ . 120021 : iPoUack et all l2009l ). where a 
represents the state which correlates with the |/ = 1, to = 
1) state at i? = 0. Because the background scattering 
length is nearly two orders of magnitude smaller for this 
^Li case than for the 843 G ^Li resonance, the Sres pa- 
rameter for the ^Li aa resonance is only 0.80 instead of 
59. Consequently, this broad ^Li resonance is tending 
towards closed channel dominance according to our clas- 
sification scheme and only has a small region of universal- 
ity spanning a few G when B is tuned near Bq. Figure [T^ 
shows the pronounced differences between the ^Li ab and 
^Li aa resonances, in spite of the similar magnitudes of 
their widths; see also Section [ILC3] 

In order to illustrate the difference between the 
"bare" and "dressed" resonance states introduced in Sec- 
tion III.A.3[ Fig. [13] shows the coupled channel bound 
state energies and scattering phases in the near-threshold 
region for thc "'"K ah channel. The a and h states 
correlate at i? = with the |/ ~ |,to 

1/ = |,m = — |) atomic states^ 

This resonance wa s observed by ( Loftus et al\ . 120021 : 
iRegal e^"^]^ 2003a ) and has additi onally been character- 
ized by ( Szvmanska et al\ . 120051 ) and ( Nygaard et al\ . 
1200611 . The actual eigenstates of the "dressed" system 
(solid lines) result from the avoided crossing of the ramp- 
ing closed channel "bare" state energy Ec = Sii{B — Be) 
and the last "bare" bound state at E-i of the background 
potential (dashed lines). The shift in the location of the 
singularity in a{B) at Bq from the threshold crossing of 
the "bare" state at Be is given from Eq. (|57)) : 



|,TO ~ — |) and 
atomic states of fermionic ^"^K. 



Bn — Br 



rhgjl - r-bg) 
'l + (l-rbg)' 



(42) 



This same formula predicts the large difference between 
Bq and Be evident in Fig. \W\ for the broad ^Li ab reso- 
nance. 

Finally, we give an example of resonances for a heavy 
species *^Rb where classification of near-threshold bound 
states using the electronic spin S = and 1 quantum 
numbers is not possible. The bin size « 0.240 GHz for the 
last bound state is much less than the ^^Rb ground state 
hyperfine splitting E^f/h = 6.835 GHz so that the last 
few bound states of the ®^Rb2 molecule are mixed by the 
hyperfine interaction. Fig. [TJ] shows the c oupled chan- 
nels s-wave bound states calculated by ( Marte et al\ . 
120021 ) for channels with M = mi -\- TO2 = 2, 1, and 0. 
Unlike the ^Li case in Fig. [S] there are a number of 
bound states within a few GHz of threshold. The lev- 
els are labeled at i? = by the spin quantum numbers 
(/1/2) of the separated atoms and the vibrational quan- 
tum number n counting down from the separated atom 
dissociation limit. Thc figure shows the last 3 vibrational 
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FIG. 12 Scattering length and entrance channel fraction for 
the ^Li 834 G open channel dominated resonance (solid lines) 
and the ^Li 737 G resonance (dashed lines), which is tend- 
ing towards closed channel dominance. The narrow ^Li res- 
onance is not shown. The upper panel shows the scattering 
length a{B) versus magnetic field strength B, whereas the 
lower panel shows the norm na{B) of the entrance channel 
spin component of the coupled channels wave function. The 
vertical lines indicate the location Bo of each resonance, and 
the horizontal double arrows indicate the width A of each. 
The dotted lines on the lower panel indicates the slope of 
na{B) = 1 — Z{B) predicted near the resonance position 
-Bo by Eq. (fSTj) of Section III. C. 51 While both resonances 
have quite large A, the ^Li 834 G one is clearly "open chan- 
nel dominated" with 1 — Z{B) remaining near unity when 
\B — Bo\ ranges over half of its width. On the other hand, the 
737 G ''Li resonance has Srcs = 0.80 and is tending towards 
being "closed channel dominated," since 1 — Z{B) drops off 
rapidly from unity as \B — Bo\ increases from resonance, with 
71^(5) > 0.5 only for |(B-_Bo)/A| < 0.11. Furthermore, this 
resonance has a universal bound state (not shown) only over 
a relatively small fraction of its width, with the calculated 
binding energy departing from Eq. ([2} by 10 percent when 
|(B-Bo)/A| =0.06. 
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FIG. 13 Bound states and scattering phase near the ah 
202 G resonance. The energy is between E/h = ±40 MHz 
{E/ks = ±1.9 mK) where the zero of energy is taken to be 
the separated atom energy of an a and a h atom. The horizon- 
tal solid line shows E^^^/h — 21 MHz. The horizontal dashed 
line shows the last "bare" bound state energy E-i of the back- 
ground potential and the sloping dashed line shows the "bare" 
resonance level with energy (5^(_B — Be), whe re = 2.35 

MHz/G. The resonance width A = 7.8 G (jGreiner et all 
I2OO3I ). The solid lines for E < indicate the coupled chan- 
nels "dressed" energies of the molecule. Away from reso- 
nance these approach the "bare" energies. The strong avoided 
crossing between the two "bare" states leads to the shift in 
the point of singularity Bo from the "bare" crossing at Be- 
For positive energies the interference between the background 
and resonant phases is evident. Since this is a open channel 
dominated resonance with Srcs ~ 2, no sharp resonant feature 
appears in sin^ '?(£', B) versus E for < Jf? < iJvdw. A sharp 
resonance feature only emerges when E ^ Bvdw. 



levels of the lowest (/1/2) = (11) separated atom limit. 
The B = energy of the (12) separated atom limit is 
Ehf = 6.835 GHz and only the ~ —4 vibrational level 
appears in this range. Similarly only the n = —5 level 
appears for the (22) separated atom limit. The closed 
channel dominated resonance (srcs = 0.17) near 1007 
G in the (11) M = 2 cha nnel has be e n used to make 
molec ules in atomic gase s ( Diirr et all l2004a[ ) and lat- 
tices (jThalhammer et all\200^ . 



C. Simplified models of resonance scattering 

While coupled channels models are very valuable for 
understanding the near threshold molecular physics of 
scattering resonances and for highly quantitative predic- 



tive calculations for a range of B field and multiple spin 
channels, they can be quite complicated to set up and 
use. Consequently, it is highly desirable to have simpli- 
fied models that are accessible to experimental and theo- 
retical researchers. Fortunately, a variety of high quality 
models are available, each valid over a limited domain of 
energy. 

The key to practical approximations for the near 
threshold bound and scattering states for ultracold neu- 
tral atom interactions is the separation of the length 
and energy scales associated with the separated atoms 
on the one hand and the molecular interactions on the 
other. The molecular interactions are characterized by 
various energy scales associated with the van der Waals 
potential, the potential at i?oxj or the minimum of the 
potential. This scale should be compared with the hy- 
perfinc, Zeeman, and kinetic energies of the ultracold 
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FIG. 14 Coupled channels *^Rb s-wave bound states. The 
dotted lines show the channel energies for four different en- 
trance channels with /i = /2 = 1, labeled by the projection 
quantum numbers (mi, 7712). The bound-state energies (solid 
lines) for these channels are shown as a function of magnetic 
field B with quantum numbers (/i/2)w assigned at i? = 0, 
where n is the vibrational quantum number in the /1/2 chan- 
nel counting down from the separated atom limit in that chan- 
nel. A Feshbach resonance occurs when a bound state with 
quantum number M — mi -\- m2 crosses a dissoci ation thresh- 
old h aving the same Af (solid dots). From (jMarte et all 
120021 '). 



atoms. For ranges of internuclear separation R where the 
molecular energy scale is much larger than the atomic 
one, the phases and amplitudes of the coupled chan- 
nels wave function components are nearly indepen- 
dent of energy and partial wave over energy range on 
the order of the atomic scale. In effect, the short range 
wave function provides an energy-independent boundary 
condition for connecting to the near threshold asymp- 
totic bound or scattering states, which arc strongly 
energy dependent. While this separation of scales 
can be made explicit in met hods based on long range 



1999al: iTsai et ai. 1997) or multichannel Quantum de- 


fect methods (Burke et ai. 


19981: Julienne and Gao. 


20061: Julienne and Miesl. 1989 


; Raoult and Miesl. 2004 


Voeels et ai. 19981). it remains imolicit as the basis for 



many other approximation schemes. 

A variety of simplified treatments show that it is suf- 
ficient to use the basic framework in Section III. A. 31 to 
parameterize threshold resonances in ultracold atoms. 
Thus, resonances are characterized by a reduced mass /i, 
a background scattering length Obg, a tunable position Eq 
selected by an external field, and an energy- and tuning- 
insensitive width Tq. Resonances that decay, whether by 
emission of light or by relaxation to lower energy open 
channels, can readily be treated by introducing the de- 
cay width 7 (jBohn and Julienn i, fl999l: iFedichev et 'dl . 



ll996bl:lHutsonl. l2007l : [K6hler et aZ.I . [20051) . The review by 
(jKohler et aiX\200^ describes how two-channel models 
can be especially effective when the resonance parame- 
ters are already known. 



1. Contact potential model 

The simplest approximation for resonance scattering is 
to use the Fermi pscudopotential ( Huang and Led . Il957l ) 



F(R) 



a{B)SiK)—R-, 

H on 



(43) 



with a strength proportional to the scattering length 
a{B). This zero-range delta-function pscudopotential is 
an excellent approximation for the full molecular inter- 
action when fc|a(_B)| <C 1 and fca ^ 1 and becomes exact 
in the limit E ^ 0. It can be used for positive or nega- 
tive a, and its phase shift is ta,nr]{E, B) = —ka{B). For 
a > it has a bound state given by the universal energy 
of Eq. @. 

If the resonance parameters Cbg, A, and Bq are known, 
the effect of tuning near a resonance can then be fully in- 
corporated using a{B) from Eq. ([T]). For an entrance 
channel dominated resonance with s^-cs ^ 1, so that 
the universal binding energy in Eq. ^ applies (see 
Eq. (I32[l ). the scattering length is the only parameter 
nee ded to treat near-thr eshold bound states and scatter- 
ing ( Kohler et al However, more robust approxi- 
mations are needed, since universality will only apply for 
detunings that at most span a range on the order of the 
width A and can be much less, depending on Sj-cs- 



2. Other approximations 

Although the underlying molecular physics often in- 
volves a number of coupled channels, many resonances 
are isolated in energy and magnetic field. Then the prop- 
erties of the "bare" resonance level are determined by 
energy scales large compared to the the small kinetic 
energies of the ultracold domain, and the level can be 
accurately approximated as coming from a single bound 
state channel, as in the Fano treatment summarized in 
Section [lI.A.3l A number of groups have developed a va- 
riety of simplified methods for characterizing the prop- 
erties of ultracold scattering resonances, but we can not 
review this work exhau st ively or in detail. F o r exa mple. 



Moerdiik efai]. Il995l). (iTimmermans et all Il999f ) and 



Kokkclmans et ali . 120021 ) introduce the standard Fesh- 
bach formalism of separating the system into bound and 
scattering subspaces, Q and P, to characterize magnet- 
ically tu nable resonances f or ground state alkali metal 
atoms. ( Goral et all |2004[ ) use a Green's function for- 
malism and introduce a separable potential that is cho- 
sen to accurately represent the two-body scattering and 
bound states of the background channel. ( Marcelis et all 
|2004[ ) are especially interested in representing the case of 
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a large neRa t ive qu f, which is relevant to the ^^Rb system. 
(|Mies et all \200d ) use the resonances of two ^'^Na atoms 



to show how to reduce a coupled 5-channel problem to an 
effective 2-channel problem using a Lennard- Jones pseu- 
dopotential with the correct van der Waals coefficient. 
(|Nvgaard et a^.l ■ l2006^ illustrates this method for the '^^K 
system. 

One model that shows great promise for practical 
and accurate fitting of resonance data is the asymp- 
totic bound sta t e (AB M) model based on the work of 
(|Moerdiik et all Il995l ). Rather than solving for the 
bound states of a set of coupled equations in order to lo- 
cate resonance positions, it uses an expansion in the last 
bound states of the Born-Oppenheimer potentials. The 
model is far less computational ly demanding tha n full 
coupled channels calculations. ( Stan et~al\ . |2004[ ) used 
a simplified version of this model to characterize mea- 
sured resonances due to the triplet m olecular state in the 
^Li+^^Na system. Recently (jWille et al. , 200§i) used this 
model to quantitatively characterize a number of reso- 
nances of ^Li-I-^'^K that involved strong mixing of the 
singlet and triplet molecular states. 



3. van der Waals resonance model 
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FIG. 15 Two-channel square well model of a magnetic Fesh- 
bach resonance. The potential for the "bare" entrance back- 
ground channel (blue line) has a well depth — Vbg- The spatial 
width is chosen to be a so as to simulate the length scale of a 
van der Waals potential. The potential for the "bare" closed 
channel (red line) has a well depth —14 relative to the sep- 
arated atoms and is infinite for R > a. The value of Vhg is 
chosen so that the background scattering length is abg. The 
value of Vc is chosen so that there is a bound state with en- 
ergy Ec near E = 0. We assume that this energy varies with 
magnetic field as Ec — 5^{B — Be), where 5n is a relative 
magnetic moment and Ec = at B = Be- 



By introducing the van der Waals Ce coefficient as an 
additional model parameter, the properties of bound and 
scattering states can be extended away from their very 
near-threshold domain into the domain where ka ^ 1 
and the binding energy is much larger than E'vdw or E. 
The reason is that there is a large range of i?, namely 
Rex < R < Rvdw, where the potential is accurately rep- 
resented as —Gq/R^ and is much larger in magnitude 
than i?vdw • The properties of the van der Waals poten- 
tial have been discussed in Section Iff. B. II 

Feshbach resonances are characterized by a width T{E) 
and shift 5E{E). These are given in the _E — > limit by 
Eqs. (|36|) and (|37|) . which depend on the dimensionless 
resonance strength Sros and rbg. The E result can 
be generalized to finite energy by introducing two stan- 
dard functions Cbg(£')~^ and tanAb g(i^) of multichan- 
nel quantum defect theory (MQDT) (lJulienne and Miesl . 



nei q uantum aeiect tne ory (Muui ) luunenne ana Mi 
llQSa iMies and Raoultl . l2000t iRaoult and Miesl . |2004| ) 



'ME) 



■C^^iE)- 



SE{E) = - tan Abg (i;). 



(44) 
(45) 



whe re for the van der W aals background poten- 
tial (jJulienne and Gaol . |2006|) 



r , 1 

2 1 + rbg) 



l + (l-rbg)^ 



(46) 



is proportional to s^^s a-nd is independent of energy. The 
MQDT functions have the following limiting form as — > 
0: ChgiE)-^ = ka{l + {1 - r^g)^) and tanAbg(£^) = 1 - 



Tbg. When E-> E, Cbg(_E)-2 ^ i and tanAbg(-E) ^ 0. 
Consequently, T{E) = T and SE{E) vanishes when E 
becomes large compared to E. If |rbg| » 1, the Chg{E)~^ 
function has a maximum and | tan Abg (E) \ has decreased 
to half its = value at £; w ?i^/(2^ (abg - a)^). 

The functions Chg{E) and tan Abg (i?), as well as 
Vhsj E), depend on on l y thr ee parameters, Cg, /i, and 
abg (| Julienne and Gaol . l2006l ) . The near-threshold phase 
r]{E) in Eq. (fT6|) can be evaluated over a wide range of 
energy on the order of E and larger from a knowledge 
of these three parameters plus Sj-cs, the magnetic mo- 
ment difference Sfi, and the resonant position Bq. The 
sin^ r]{E, B) function evaluated using Eqs. (|44|l and (j45|) 
are virtually indistinguishable from the coupled channels 
results shown in Figs. UmiT^ 



4. Analytic 2-channel square well model 

A very simple square well model, because it is analyt- 
ically solvable, can capture much of the p hysics of near - 
threshold bound and scattering states. (iBethd . I1935D 
used such a model to successfully explain the threshold 
scattering of cold neutrons from atomic nuclei, where the 
neutron de Broglie wavelength was very large compared 
to the size of the n ucleus. (jKokkelmans et I2002D 



and (|Duine and Stoot . .2004) have introduced 2-channel 
square well models to represent Feshbach resonances in 
ultracold atom scattering. 

Figure [15] shows the "bare" background and closed 
channel potentials for a square well model where for con- 
venience of analysis we take the width of the wells to 
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be the van der Waals length a. The background en- 
trance channel and the closed channel are designated 
by \bg) and |c) respectively. Using a two-state cou- 
pled channels expansion as in Eq. \tl.'{R,E)) = 
\c) (f>ciR, E) / R + \bg)4}hgiR, E)/R, the potential matrix 
in Eq. gO]) is 



~Vc W 
W -14g 



for R < a 



(47) 




for R> a, 



The off-diagonal matrix clement W describes the weak 
coupling between the two channels. 

In order to simulate a magnetically tuned Fcshbach 
resonance, the model parameters need to be chosen so as 
to give the correct parameters for that resonance. The 
well depth Vbg is chosen so that the background channel 
scattering length is Obg. The well depth Vc is chosen so 
that the well has a "bare" bound state at Ec- The tuning 
of the bound state as Ec = Sfj,{B — Be) can be simulated 
by varying Vc linearly with the external magnetic field 
B. Finally, weak couphng requires \W\ ^ |Vbg — 141. 
The coupling parameter W can then be chosen to give 
the right resonance width T(E) = 2kahgSfJ.A at low en- 
ergies (see Eq. (|17p ). using the known resonance width 
A. Analytically calculating the matrix element defining 
r(£:) in Eq. ^ relates to A as follows. 



''bg 



[V^g - Vc^ (1 - r^J' 



(48) 



With the chosen parameters, the square well model yields 
analytic form for the scattering phase shift as in Eq. (jl6[) 
and the scattering length as in Eq. ([1]). 

The square well model also permits an analytic evalu- 
ation of the weakly bound state below the continuum. 
Assuming an eigenstate lipi,) exists at energy — E't, = 
—fi^kll {2^) < and |abg| ^ a, we get 



1 



rsq/2 



Ubg - a a{Eb + Ec) 



(49) 



where rsq/2 (5^Arbg(l-rbg)~^ (jMarcelis et a/.l . l2004l ) 
derived a similar result for a contact potential. Note 
that when the coupling term ^ so that Egq 

0, the solutions Et = —Ec and Et = h? /[2^{abg — a)^] 
correspond to the bare states of the square well in the 
closed and open channel (for abg > a), as expected. Since 
the resonant singularity in the scattering length occurs 
when Eh 0, taking this limit of Eq. allows us to 
calculate the resonance energy shift 5E ~ Sii{Bq — Be) 
as 



6E =^il-r,,). 



(50) 



Both Eq. ([501) ^■'^d Eq. can also be derived from the 



('^bg — 1)~^]- Note that Tsg and f are nearly the same for 
{rbgl ^ 1. The modified version of Eq. (|50p is equivalent 
to Eq. ([37]) and Eq. ([42]) , derived from the van der Waals 

m odel. 

( Lange efai] . l2009l) extends the above model to pre- 
cisely determine the scattering length and the resonance 
parameters in the magnetic field regime where multiple 
Feshbach resonances overlap. 



5. Properties of Feshbach molecules 

A variety of properties of Feshbach molecules can be 
calculated by solving Eq. for the binding energy Ei, . 
For example, the closed channel fraction Z of the eigen- 
state can be found by differentiating Eh with respect to 
Ec; see Eq. ((25)) . In the limit B ^ Bo where Ei, vanishes 
and a — > -t-oo, we have 



B-Bn 



1 



where the dimensionlcss proportionality constant 



Fbgl 



'bg 



E 



(51) 



(52) 



determines the rate at which the Feshbach molecular 
state deviates from the entrance channel dominated 
regime or, cquivalcntly, the halo molecule regime, when 
B is tuned away from Bq. Eq. (|5ip shows that having 
a small closed channel fraction Z <^ 1 requires the mag- 
netic field to be close to resonance, \B — Bo\ <C Cl^l- 
Figure [T^ of Section ITl.B. 51 compares 1 — Z from Eq. ((5T|) 
for the respective open and closed channel dominated ^Li 
834 G and ^Li 737 G resonances. 

For open channel dominated resonances, where s^cs S> 
1, it is usually true that |rbg| > 1 and C 3> 1, and Z re- 
mains small over a large fraction of the resonance width 
A. The bound state wave function takes on primarily en- 
trance channel character over this range. See the exam- 
ples of the f^Li 834 G or ^^K 202 G resonances in Fig. [TH 
Closed channel dominated resonances have Srcs ^ 1 and 
small C 1- Gonsequently, Z remains small only over 
a small range of the resonance; see the example of the 
s^Rb 1007 G or ^Li 543 G resonances in Fig. [TH A smaU 
Cm = 0.09 for the Rb 1007 G resonance implies that 
the molecular state is entrance channel dominated only 
within ^ 9% of the resonance width. Table IIVI in the 
Appendix lists C for several other resonances. 

Expanding Eq. (|49)) at small binding energies, the 
molecular binding energy has the following form in the 
threshold limit: 



Eh 



2fi{a - a + R*y^ 



(53) 



van der Waals model with Tsq replaced by F 



^sg[l 



where R* ~ a/s,-^s- This expression applies in the limit 
that a ^ a and a 3> 4i?*. The binding energy Ei, shows 
two corrections to the universal l/a^ threshold law in 
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FIG. 16 (Color) Closed channel fraction Z (panel a) and en- 
ergy —Eb (panel b) of weakly-bound molecular states as a 
function of magnetic field for selected Feshbach resonances. 
The open circles come from coupled channels calculations as 
described in Section III. B. 31 The solid curves are calculated 
based on Eq. (|49)) . Eq. ([28}, and the parameters in Table HVl 
The slope of Z at B = Bo defines the C parameter; an ex- 
ample for the ^''Rb 1007 G resonance is shown as the dashed 
line with Cflt = 0.09. This ^'^Rb resonance and the ''Li 543 G 
one are closed channel dominated and have Sred,C < 1- The 
® Li 834 G and the '^°K 202 G resonances are open channel 
dominated resonances with Sres,C > 1- 



Eq. (21). One is the finite range correcti on a shown i n 
Eq. (1211); the other one, R*, introduced by (|Petrovl . [200l . 
is unique to Feshbach resonances. The correction R* is 
negligibly small for open channel dominated resonances 
with Sios ^ 1- Closed closed channel dominated reso- 
nances can have large R* 3> a. Such resonances only 
have a regime with a bound state of predominant en- 
trance channel character near Bq where a ^ 4i?*. This 
condition is consistent with the one given in the previous 
paragraphs, namely, \B — Bo\ <^ (^\A\. 



III. FINDING AND CHARACTERIZING FESHBACH 
RESONANCES 

Magnetically tunable Feshbach resonances have been 
experimentally observed in essentially all alkali-metal 
species, in some mixtures of different alkali-metal atoms, 
as well as for Cr atoms. Experimental data on the reso- 
nance positions have in many cases enabled the construc- 
tion of accurate models to describe the near-threshold 
behavior, including scattering properties and molecular 
states. We will first review the various experimental 
approaches to identify and characterize Feshbach reso- 
nances in Sec. IIII.AI and then discuss observations of 
Feshbach resonances for alkali-metal atoms in Sec. IIII.BI 
and various other systems in Sec. IIII.Cl 



A. Experimental methods 

Experimental approaches to detect magnetic Fesh- 
bach resonances can be classified into several types. 
After some general considerations in Sec. IIII.A.li 
we will discuss detection by inelastic collisional trap 
loss in Sec. IIII.A.21 by elastic collision properties in 
Sec. IIII.A.3[ and loss in the presence of optical radiation 
in Sec. IIII.A.4I Finally, Section IIII.A.5I discusses preci- 
sion radio- frequency spectroscopy of Feshbach molecules. 



1. General considerations 

a. What is the magnetic field range to be explored? The 
typical spacing between two Feshbach resonances can be 
estimated from the ratio of the energy splitting between 
closed channel molecular levels and the relative magnetic 
moment 6fi between the entrance channel and the closed 
channel. The vibrational energy splitting between near- 
threshold bound states is determined by the long-range 
van der Waals potential to be on the order of 100i?vdw; 
see Sec. IILB.II For alkali- metal atoms is on the order 
of two Bohr magneton 2^b = 2.8MHz/G and 6 times 
larger for ^^Cr. For atoms with a small hyperfine split- 
ting compared to lOOSvdw, Feshbach resonances are in- 
duced by the last bound states. This leads to a typical 
s-wave Feshbach resonance separation of ~ 10, 000 G for 
^Li and ^ 100 G for ^^Cr. For atoms with hyperfine 
splittings much larger than 100i?vdwj resonances can be 
induced by much deeper bound states in the closed chan- 
nel (see Fig.[T3]for ^^Rb), and the expected spacings can 
be estimated accordingly. 

The density of Feshbach resonances increases when 
higher partial-wave scattering and multiple closed hyper- 
fine channels, defined in Sec. III. B. 31 are included. The 
relevant number of channels is determined by the angu- 
lar momentum dependence of the molecular potentials 
and identical particle statistics. For alkali-metal atoms 
there are on the order of 10 closed channels for the low- 
est partial wave ^ = 0. The number of channels increases 
rapidly for higher partial waves. The ultralow tempera- 
ture usually limits scattering to s- and sometimes p-wave 
entrance channels, but coupling to molecular states with 
up to = 4 has been observed; see Tables [Tl] and IIVI 

For species without unpaired electrons, e.g., Sr and 
Yb, one expects no or very limited magnetic tunability 
because there is no electron contribution to the magnetic 
moment and the nuclear contribution is very small. In 
these systems, Feshbach resonances can possibly be op- 
tically induced; see Sec. IVI.AI 

b. What is the required magnetic field resolution? The 
width of the resonance generally determines the mag- 
netic field resolution required for detection. Many s-wave 
Feshbach resonances have widths larger than 1 G; see 
Table IV in the Appendix. High partial wave Feshbach 
resonances are typically much narrower because of the 
weaker Feshbach coupling strength. Usually, a resolution 
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in the milligauss range is required to detect d- or g-wave 
Feslibach resonances. 

c. How to trap atoms fo r collision studies? Optical 
dipole traps, reviewed by ( Grimm et all . [2OOOI ). are the 
main tool to confine cold atoms for collision studies re- 
lated to Feshbach resonances. Optical potentials trap 
atoms in any sub-level of the electronic ground state and 
permit investigation of collisions in any corresponding 
spin channel. For many experimental applications, the 
lowest atomic state is of particular interest, which is a 
high-field seeking state and can therefore not be trapped 
magnetically. Optical dipole traps allow for the appli- 
cation of arbitrary homogenous magnetic fields without 
affecting the trapping potential. In contrast, magnetic 
traps can only confine atoms in low-field seeking states, 
and the application of a magnetic bias field for Fesh- 
bach tuning can strongly influence the trap parameters. 
This limits the application of Feshbach tuning in mag- 
netic traps to very few situations. 

d. How low a temperature is needed to observe the res- 
onances? In most experiments, a temperature of a few 
fiK is sufficiently low to observe a clear resonant structure 
in inelastic coUisional loss. Collision studies can be per- 
formed with thermal samples, BECs or degenerate Fermi 
gases. Elastic collision measurements are more complex. 
Enhancement of elastic collision rates near Feshbach res- 
onances is more prominent at lower temperatures < l^K. 
On the other hand, suppression of elastic collision rates 
due to a zero crossing of the scattering length can be 
easily seen well above l^K; see Sec. IIII.A3] 



2. Inelastic loss spectroscopy 

Resonant losses are the most frequently observed signa- 
tures of Feshbach resonances in cold-atom experiments. 
These losses can be induced by two-body or three-body 
processes. Loss occurs because of the release of internal 
energy into the motion when colliding atoms end up in a 
lower internal state or when a molecule is formed. The 
gain in kinetic energy is on the order of the Zeeman en- 
ergy, the hyperfinc energy, or the molecular vibrational 
energy, depending on the inelastic channel, and is gener- 
ally so large that all atoms involved in the collisions are 
lost. Near a Feshbach resonance, inelastic loss is strongly 
enhanced because the Feshbach bound states have strong 
couplings to inelastic outgoing channels. 

Two-body and three-body collision loss can be quan- 
tified based on the evolution of the atom number N{t), 
which for a single species satisfies 



N(^t) = - J [L2n\r, t) + L^n^r, t)] d\, 



(54) 



where r is the one-body lifetime, typically determined 
by background gas collisions, n(r, t) is the position- and 
time-dependent atomic density distribution and L2 (L^) 
is the thermally averaged two-body (three-body) loss co- 
efficient. 
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FIG. 17 Time evolution of the peak atomic density in a cloud 
of ^^^Cs atoms in the |/ = 4,m = -4) state at T = 5.3 ^K. 
The solid circles show the off-resonant evolution at B = 140 G 
whereas the open circles show the on-resonant evolution at 
B — 205 G (open circles), where a Feshbach resonance is lo- 
cated. The fit is based on Eq. (|55p wi th L3 — 0, and the peak 
density corresponds to 2\/2w. From (jChin et oLL l2004bl l. 
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FIG. 18 Two-body inelastic loss coefficient of cesium atoms in 
the \ f — 3,m = —3) state as a function of magnetic field. The 
loss coefficient is extracted from the atomic density evolution, 
as shown in Fig |17l Three resonanc es are ident ified here. The 
solid line is a Lorentzian fit. From (|Chin et a/.i . i2004by 



The loss equation can be further simplified under the 
assumption that thermalization is much faster than in- 
elastic loss. For example, for a thermal cloud with tem- 
perature T in a 3D harmonic trap, one finds 



nit) 



n{t) 



-L^nitf -{AlZf/^L^n{t)\ (55) 



where n = Nu}^{4:nkBT/m)^^^'^ is the mean density, 
TO the atomic mass, and uj the geometric mean of the 
three trap vibrational frequencies. Examples of density- 
dependent loss curves are shown in Fig. [T71 

The trap loss coefficient L2 is related to the inelastic 
loss coefficient in Sec. III.AT2] bv L2 ~ K\oss{T), where we 
have assumed both atoms are lost in one collision event. 
Near a Feshbach resonance, L2 is enhanced and has a 



24 



Lorentzian profile at low temperatures; see Sec. IILA.3I 
for more details on inelastic scattering resonances. Two- 
body collision loss has been observed in many cold atom 
system; see Fig. [T5] for an example. 

Three-body loss, as described by the loss coefHcient 
L3 in Eq. ((54)) . is also strongly enhanced near Feshbach 
resonances. In many experiments cold atoms are po- 
larized in the lowest ground state and two-body inelas- 
tic collisions (in the aa channel) do not occur so that 
three-body recombination loss is the dominant trap loss 
process. Three-body recombination occurs when three 
atoms interact and form a diatomic molecule and a free 
atom. In this process, the molecular binding energy is re- 
leased into the kinetic energies of the outgoing molecule 
and the third atom, which except for very small molecu- 
lar binding energies leads to immediate trap loss. 

In the fir st experimental rep ort on atomic Feshbach 
resonances, (jlnouve et aLl . ll998[) observed very fast trap 
loss of a sodium BEC near a Feshbach resonance. In this 
experiment, three-body recombination is the leading trap 
loss process. Recombination losses induced by Feshbach 
resonances have been observed and studied in numerous 



later experiments, for example, those by (IRoberts et al. 
2OOOI) on 85Rb, by (iMarte et all l2002l: ISmirne et al. 



20071) on 87Rb and by (IWeber et al\ . l2003aD on ^^^Cs. 

For bosonic atoms with large scattering length 
a ^ a and low temperatu r es, L 3 scales ge n eralh 
as a* 



2 , ana low temperatu r es, scales ge n erally 
(iBraaten and Hammed. l2006l: lEsrv et ali [Tooa 



iFedichev et all Il996bl : iNielsen and Macekl . Il999l ). but 
with additional quantum features (resonance and in- 
terference effects) as discussed in Sec. I VI. CI on Efimov 
physics. For fermionic atoms the situation is more 
complicated because of Pauli suppression effects; see 
Sec. IIV.B| but generally a loss feature accompanies a 
Feshbach resonance. 



3. Elastic collisions 

Elastic collisions refer to scattering processes in which 
the colliding atoms only change their motional state, but 
not the internal state. The cross section CTei for elastic 
s-wave collisions follows from Eq. Neglecting the 

effective range correction (rg = in Eq. (O), one obtains 
the simple expression 



(56) 



where a depends on magnetic field B and g is the sym- 
metry factor introduced in Sec. III.AT2] Near a Feshbach 
resonance, the scattering length becomes very large. The 
elastic s-wave cross section is g Atto^ at very low energy 
as /c — > and approaches its upper bound gAn/k'^ in the 
unitarity limit at finite k where /ca > 1. The latter can 
be reached in the regime if a becomes very large. 

A strong enhancement of the thermally averaged elas- 
tic collision rate n{aeiv) can indicate the occurrence of 
a Feshbach resonance. One experimental approach is to 
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FIG. 19 Fraction of trapped atoms remaining after evapora- 
tion of ^Li atoms in an optical dipole trap after 1 s (open cir- 
cles) and 3 s (closed circles) hold times. The peak near 530 G 
indicates a minimum of evaporative loss. This minimum re- 
sults from the zero crossing of the scattering l ength induced 
by th e broad 834-G Feshbach resonance. From (|Jochim et all 
|2002| ). 



measure the thermalization rate, which is proportional 
to the elastic collision rate as Kn{aciv). Here k is nu- 
meric ally calculated as 2.7 in the low temperature A: — > 
limit ([Monroe et al. I, [1991), 10.5 in the unitarit y limit 
([Arndt et ali Il997t) . Similarly, ([PeMarco et all Il999l ) 



found K =4.1 for p-wave collisions in Fermi gases. Find- 
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Near a resonance the thermalization rate can be lim- 
ited under hydrodynamic conditions, which are reached 
when the cross section is so large th at the collision rate in 
the trap exceeds the trap frequency ( Vuletic et aLl . [T999i) . 
The maximum collision rate is also bounded by the uni- 
tarity limit when a is large. In both cases resonance 
structure that is due to clastic scattering can become 
less evident near Bq. 

Another efficient method to identify Feshbach reso- 
nances based on elastic collisions is to locate the zero 
crossing of the scattering length, that is, the magnetic 
field for which the scattering length vanishes near res- 
onance where B = Bq + A; see Eq. ([T]). This can be 
monitored by measuring atom loss resulting from elas- 
tic collisions during the process of evaporation. Ther- 
malization and evaporation loss are suppressed at the 
zero crossing. Schemes t o locate the zero cro ssing have 
been apphed to ^^ Rb by (iRoberts ad. 1 19981) 



Chin et aLl.[2000l). ^°K by (iLoftus et al 



Jochim et aZ.! . 20021: O'Hara et al 



®^Rb mixture by ( Zaccanti et all 2006f ). An example is 



20021) . 



Cs by 
^Li by 



2002bD and to a ^'^K 



shown in Fig. While the zero crossing is evident, no 
resonance feature is seen near Bq = 834 G. 

Finally, resonant changes of elastic scattering can also 
be revealed through the d etection of co l lision shifts in 
atomic clock experiments ( Marion et all |2004[ ) and by 
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measurements of the mean-field interaction in Bose- 
Einstein condensates (ICornish et "all. l2000t[hrouve et all 
119981 : iRegal et all l2003bD; see discussion in Sec. IIV.A.2I 



4. Radiative Feshbach spectroscopy 

Radiative Feshbach spectroscopy makes use of red- or 
blue-detuned light to detect the variation of the coUi- 
sional wave function near a Feshbach resonance. The 
amplitudes of both the open and closed channel com- 
ponents of the wave function are strongly modified for 
distances on the order of or less than a when B is tuned 
near resonance. This is evident for the closed channel 
component since the outer turning point will be on the 
order of a; see Sec. lII.BTTl The optical transition induces 
loss of atoms by excitation of atom pairs at such separa- 
tionS; 

(|Courteille et all Il998[ ) adopted the idea of radia- 
tive spectroscopy to identify a Feshbach resonance in 
^^Rb. In this exper iment, a photoassociation laser beam 
( Jones et all l2006f ) is held at a fixed frequency to the 
red of the strong 5 — s- P atomic transition and serves as 
a sensitive probe to measure the resonance position Bq. 
The light excites the colliding pair of atoms to an excited 
molecular level in a state with an attractive potential. 
The excited level decays by spontaneous emission, giving 
rise to atom loss. The experiment monitors the atom loss 
as B is varied n ear Bq, thu s locat ing the resonance. 

In contrast, ( Chin et oZI . |2003| ) applied a laser with 
far blue detuning to detect Feshbach resonances in ce- 
sium samples. The blue-detuned light excites a molec- 
ular state with a repulsive potential so that the atoms 
are repelled, accelerated from one another, and lost from 
the trap. This method requires less detailed knowledge 
of molecular structure than the previous method. In this 
experiment, multiple narrow d- and g-wave resonances 
were identified in two different collision channels of ce- 
sium atoms; see Fig. [501 



5. Binding energy measurements 
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FIG. 20 Radiatively induced atom loss versus magnetic field 
due to a far blue-detuned laser beam applied to a sample of 
^^^Cs atoms confined in an optical dipole trap at a temper- 
ature of 3 nK. The laser wavelength is tuned 5 nm above 
the free atomic transition to excite the colliding atoms to a 
repulsive molecular state. Once excited, the atoms quickly 
dissociate and are lost from the trap. Multiple narrow d- 
and g-wave Feshbach resonances were observed. The inset 
shows an expanded view with more resonances resolved. From 
(|Chin et ali\200^ . 



bound" transition, which converts them into a different 
molecular s tate. Based on the lines hape functions cal- 
culated by ( Chin and Julienn"^ . |2005[ ). binding energy of 
Feshbach molecules can be measured to 1 kHz. This kind 
of precise data can be combined with theoretical model- 
ing to determine the position and the width of the reso- 
nance. 



Other efforts to spectroscopically probe weakly bound 
states inclu de oscillating magne t ic field spectroscopy em- 
ployed by ( PapD and Wiemanl . l2006l : iThompson et al 



2005a ) , and rf and microwa ve techniques by ([Mark et al. 



2007aHZirbel et aZ.I, [2008bl] . Using the t heoretical models 
described in Sec. III.C.4[ jhanse et all l2009t ) show how 
Feshbach resonance parameters can be extracted from 
molecular binding energy measurements. 



The detection methods discussed in the previous sub- 
sections provide good ways to determine the existence 
of resonances and their positions. Measurements of 
the magnetic-field dependent binding energies of near- 
threshold Feshbach molecules can yield information to 
precisely determine the scattering properties near a spe- 
ci fic resonance; see F ig . [2] an d Sec. fII.C.5l 

( Bartenstein et all l2005l : iRegal et all l2003a[) em- 
ployed radio-frequency (rf) spectroscopy on Feshbach 
molecules to measure very small molecular binding en- 
ergies. An example of the rf spectroscopy is shown in 
Fig. 1211 In this experiment, weakly bound molecules are 
first prepared near the Feshbach resonance; see Sec. IV.Al 
An rf field is then applied to drive either a "bound-free" 
transition, which dissociates the molecules, or a "bound- 



B. Homonuclear alkali-metal systems 

Feshbach resonances have been found and character- 
ized in essentially all single-species alkali-metal systems. 
The scattering properties show vast differences between 
the various species and also between different isotopes of 
the same species. Each system is unique and has partic- 
ular properties. Here, we give a brief account for each 
species or isotope, ranging from first observations to the 
best current knowledge, and we discuss the characteristic 
properties. A table of important resonances can be found 
in the Appendix (Table [IV|. 
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FIG. 21 Investigation of Feshbach resonances in the low- 
est two scattering channels in ^Li with radio-frequency (rf) 
spectroscopy on ultracold molecules. In this experiment, 10^ 
molecules, initially prepared in the ab channel and confined 
in a dipole trap, are radiatively excited to the ac channel. 
The rf excitation lineshape shows a markedly different be- 
havior below and above 690 G. Panel (a) shows the molec- 
ular bound state (solid) and channel energies (dotted) for 
the two channels. Below 690 G, a weakly bound state ex- 
ists in the ac channel, and a bound-bound transition is pos- 
sible (blue arrow). Above 690 G, no weakly bound state 
exists in the ac channel, and only a bound-free transition 
that dissociates the molecules is possible (blue arrow). Panel 
(b) shows the observed narrow and symmetric bound-bound 
excitation lineshape at 661.4 G. Panel (c) shows the ob- 
served broad and asymmetric bound- free lineshape at 720.1 
G. The disappearance of a bound state in ac indicates that 
a Feshbach resonance oc curs in this channel. Adapted from 
(jBartenstein et a/.l[2005h . 



1. Lithium 

Lithium has two stable isotopes, one fermionic (^Li) 
and one bosonic (^Li). 

'^Li - In 1998 (|Houbiers et all Il998[ ) predicted Fesh- 
bach resonances in cold collisions of fermionic ^Li. Ex- 
perimental evidence for a prominent re sonance was estab- 
lished by monitoring inelastic deca y (iDieckmann et al 



2002) and elastic collision properties (jJochim et aLl . l2002 ; 



O'Hara et all [20021J ). These experiments showed large 
variations in loss and thermalization rates as a function 
of the magnetic field strength for equal mixtures of the 
lowest two hyperfine ground states a and b prepared in an 
optical dipole trap. Because of the fermionic nature of ^Li 
even partial-wave scattering and, in particular, s-wave 
scattering only occurs between atoms in unlike hyperfine 
states. Consequently, if the gas is sufficiently cold such 
that the £ > centrifugal barriers are higher than the 
temperature of the gas, s-wave collisions in the ab chan- 



nel represent t he es s ential therma l ization mechanism. 
([Jochim et all l2002t lO'Hara et oil . l2002b[ ) reported a 
strongly suppressed thermalization rate near 530 G, in- 
dicating the zero crossing o f the scattering l ength in the 
ab channel. Most recently, ( Du et ai . 20081) located the 
zero crossmg to 527.5 ± 0.2 G. (jPieckmann et all I2002D 
observed enhanced inelastic loss near 680 G, about 150 G 
below the actual resonance location {Bq ~ 834 G). Note 
that this particular resonance is extrao rdinarily broad 
(A — 300G). Fu r ther experiments by ( Bourdel et all 
120031 : iGupta et all |2003[ ) provided confirmation of this 
Feshbach resonance. 

In order to furt h er pin point the resonance position, 
( Bartenstein et Ml, I2005D conducted radio-frequency 
spectroscopy on Feshbach molecules in the ab channel 
as described in Sec. IIII.A.41 The resonance could be lo- 
cated with an uncertainty of 1.5 G, more than two orders 
of magnitude smaller than its width. Broad resonances 
in two other entrance channels (ac and be) were reported 
as well. 

( Strecker et all I2003D identified a narrow s-wave reso- 
nance a^_543;25_G in thBo channel, and ( Schunck et all 
120051: IZhang et all l2004^ observed three p-wave reso- 
nances near 200 G, one in each of the channels aa, ab, 
and bb. All of these observed resonances are induced by 
s- and p-wave rotational levels of the u 38 bound state 
of the singlet X^I]+ potential, as shown in Fig. [5] for s- 
wave resonances. 

The two s-wave resonances in the ab channel illus- 
trate the concept of resonance strength, as described in 
Sec. III. B. 21 The broad resonance at 834 G is strongly 
entrance-channel dominated, while the narrow resonance 
at 543 G is an extreme case of a closed-channel dominated 
resonance. The former has an extraordinarily large mag- 
netic field range of universal behavior, while that for the 
latter is vanishingly small (detunings of a few fiG or less) ; 
see Sec. III. C. 51 and Fig. [1^1 This has important conse- 
quences for molecules formed near these resonances. 

The very broad resonance at 834 G in the ab channel is 
used to create molecular Bose-Einstein condensates and 
fermionic superfluids. This will be extensively discussed 
in Sec. 

""Li - (|Moerdiik and Verhaarl . Il994[ ) predicted Fesh- 
bach resonances in collisions of bosonic ^Li at o ms in 
1/ = l,m = 1) {aa channel). (|Strecker et all |2002| ) 
identified a zero crossing of the scattering length in- 
duced by a broad Feshbach resonance. This resonance 
was used to study the formation of bri ght so litons by 
(|Khavkovich et all l2002l : IStrecker et oil . 12002^ at small 
and negative scattering lengths near the zero crossing; 
see also Sec. IIV.A31 In these early experiments, the res- 
onance po sition was est i mated to 7 20 G. Later measure - 
ments by (jjunker et all |2008|) and (|Pollack et all l2009l ) 
accurately pinpointed the resonance position and the zero 
cr ossing to 736.8(2) G and 543.6( 1) G, respectively. 

( Gross and Khaykovichl |2008|) observed two reso- 
nances in the state |/ = l,m = 0) {bb channel). A 
narrower one was found at 831(4) G with a width of 7 G, 
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while a broader one (width 34 G) was located at ~ 884 G 
. In between these two resonances, a zero crossing was 
found at 836(4) G. 



2. Sodium 



{ab channel). One year later, t he same group repor ted on 
a p-wave resonance at 199 G ( Regal et all l2003bl ); they 



( Inouve et all Il998l ) pioneered experimental research 
in locating Feshbach resonances. In an optically trapped 
BEC of ^^Na (the only stable isotope) with all atoms in 
the lowest hypcrfinc state \f = l,m — 1), they identi- 
fied resonances at 853 G and 907 G. A third resonance 
i n a different c hanne l at 1195 G was later reported by 
(jStenger et all Il999t ). All three resonances are narrow 
and s wave in nature. The 1998 experiment showed both 
strongly enhanced trap loss and the dispersive tuning 
of the scattering length near the 907-G resonance; see 
Fig. [31 

The experimental determination of the resonances has 
enabled detailed models of the interaction potentials 
betwee n ultracold Na atom ( van Abeelen and Verhaail 
ll999aD. These models were further refined by 



( Samuelis et Idl. l200nt ) based on conventional molecular 
spectroscopy in combination with photoassociation data. 

The Feshbach resonan ce at 907 G wa s used to create 
ultracold Na2 molecules ( Xu et aJ} . l2003l) and to demon- 
strat e coherent molecule optics ( Abo-Shaeer et all . 
I2005D : see also Sec. lV.A.ll 



measured the resonantly enhanced elastic collision rate 
of atoms in the second-lowest hyperfine state (bb chan- 
nel) at a temperature of 3 //K. At an even lower tem- 
perat ure and by monitorin g the collision-induced heating 
rate, (jTicknor et al. I. I2004D found that this resonance is 
actually a doublet. This doublet structure in the p-wave 
resonance is due to a small energy splitting between the 
\£ = 1,771^ = 0) and \£ = l,me = ±1) molecular states. 
The anisotropic nature of the p-wave resonances has 
found interest i ng ap plications in low-dimensional traps 
(iGiinter et all\2005) . 

In spin mixture s with s- w ave i nteractions, 
(|Regal and Jinl . l2003t iRegal et all l2003b[ ) identified 



10 G wide s-wave Feshbach resonances in the ab and ac 
channels at 201.6(6) G and 224.21(5) G, respectively, 
by monitoring the thcrmalization rate and mean-field 
shifts. These resonances provide a convenient tool to 
study strongly interacting Fermi gases and fermionic 
condensates; see Secs lIV.Bl 

"^ ^K - A Feshbach reso nance was recently observed 
by ( Kishimoto et aZ.I . |2009[ ) at 51.4 G in the cc channel. 
The observation confi rmed a theoretical prediction by 
( D'Errico et al I I2007D . which was based on experimen- 
tal data available for the other potassium isotopes. 



4. Rubidium 



3. Potassium 

Potassium has three stable isotopes, two of them are 
bosonic ('^^K and ^^K) and one i s fermi onic (^°K). 

- In 1996 (|Boesten et all Il996bh predicted Fesh- 
bach resonance locations in collisions between ^^K atoms. 
Their results were based on spectroscopic data of bind- 
ing energies of ro-vibrational states of the X^S+ and '^S+ 
potentials. The data, however, were not sufficiently com- 
plete to give quantitative resonance locations, but it did 
show that the likelihood of resonances was large. An- 
oth er early prediction of resonance locations was made 
by (|Bohn fflUll999l) . 

Feshbach resonances in ^^K were observed at 402 G 



and analyzed by (jP'Errico et all 12007), and applied t o 



p iiea t o 

create a tunable BEC of this species (|Roati et adl2007t) . 
The zero crossing of the scattering length near the broad 
402-G resonance has found intriguing applications for 
atom interfer o metry with non-interacting condensates 
(iFattori et aLl . l2008a[) : sec also Sec. IIV.A.4I This reso- 
nance has an intermediate character between that of an 
entrance channel and a closed channel dominated reso- 
nance. 

""^K - Early predictions on Fesh bach resonances in 
were made by (|Bohn et all Il999l). The first experimen- 



tal observation was reported by i Loftus et aLl . [200^ . who 
demonstrated resonant control of elastic collisions via the 
202-G resonance in a mixture of the lowest two spin states 



^^Rb - (|Courteille et all Il998l) reported a Fesh- 
bach resonance in a magnetically trapped thermal sam- 
ple by observing enhanced photoassociative loss; see 
Sec. IIILA.41 This resu lt confirmed the p r edicti on by 



(|Vogels et all Il997[ ). (|van Abeelen et all Il998[ ) sug- 



gested using photoassociation as a probe to identify 
Feshbach resonances . By monitoring inelastic loss, 
( Roberts et aZ.I . [l998l ) determined the position of this 10 - 
G wide resonance to be at 155 G. ( Claussen et aLl . [20031 ) 
then used a BEC to perform a high-precision spectro- 
scopic measurement of the molecular binding energy and 
determined the resonance parameters Bq and A within 
20 mG. 

Attainment of BEC in '^^Rb crucially depended on 
the existence of the 155-G resonance; see Sec. IIV.A.II 
Only in a 10-G window nea r the resonan c e the scat- 
tering length is positive and ( Cornish et all [2OO0I ) were 
able to Bose condense ^^Rb. As the first available 
BEC with widely tunable interactions the system has re- 
ceived considerable attention; see Sec. IIV.A.3I Coherent 



atom- molecule coupling (jPonlev et a/.l . l2002D and the for- 
mation of ultracold *^Rb2 molecules ( Thompson et all 
l2005a[ ) wer e reported based on this Feshbach resonance. 

**^Rb - (jMarte et all |2002[ ) conducted a systematic 
search for Feshbach resonances in bosonic ®^Rb. For 
atoms polarized in various combinations of magnetic sub- 
levels in the / = 1 hyperfine manifold, more than 40 
resonances were observed between 300 G and 1200 G by 
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monitoring atom loss in an optical dipole trap. These res- 
onances are induced by s- and d-wave bound states and 
are all very narrow. For the s-wave states, the underlying 
molecular structure is shown in Fig. 1141 The widest and 
most often used resonance is located at 10 07 G and has a 
width of 0.2 G. In differe nt experiments, ( Erhard et all 



120041 : IWidera et aLl . l2004h observed a low-field resonance 
near 9 G in the ae channel. 

Several resonances in ^^Rb have been used to form ul- 
tracold Feshbach molecules. They are form ed by ramp- 
ing the magnetic field through the resonance ( Diirr et all . 
l2004a[ ): see Sec. IV.A.ll The great potential of com- 
bining Feshbach resonances with optical lattices has 
been demonstrated i n a s e ries of experiment s with 
^^Rb (ISvassen et all l2007t iThalhammer et aZi j2006t 
IVolz et aLl . 120061 : IWinkler et aZ.l . l2006D : see Sec. lVLBl 
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5. Cesium 



FIG. 22 Scattering length and bound state energies for ce- 
sium atoms in the lowest internal sta te as a function of mag- 
netic field. From (jChin et aZ.I . l2004bl ). 



Cesium, for which the isotope '^■^Cs is the only 
stable one, was proposed as the first alkali-metal 
species in which Fesh bach resonances could be observed 
(|Tiesinga et all\l99^ . With the limited state of knowl- 
edge of the interaction potentials no quantitative predic- 
tions could be made. The first observation of Feshbach 
reson ances in cesium colli sions was published seven years 
lat er (IVuletic aZ.I ITooj). 

( Chin et am2000ll2004bl : IVuletic et aLl . ll999f) reported 
on more than 60 Feshbach resonances of various types us- 
ing several detection schemes. These include s-, p-, d-, j- 
and (7- wave resonanc es in 10 different scattering channels. 
(|Vuletic et aZ.l . ll999l ) used cross-axis thermalization rates 
to identify resonances in the lowest aa chan nel and used 
trap loss measurements in the gg channel. ( Chin et all 
l2000t ) reported many more resonances by preparing the 
atoms in other internal states and by monitoring the 
evaporation rates to more efficiently m easure the elastic 
cross section. The resonances provided ( Leo et all . l2000f ) 
with the essential information to precisely determine the 
interaction potentials of ultracold cesium. 

Narrow r esonances were observed by ( Chin et all 
l200l l2004d) using radiative Feshbach spectroscopy; see 
Sec, nil. A. 41 In these experiments ^^'^Cs atoms were illu- 
minated by a far blue-detuned laser beam, whose wave- 
length was optimized to only remove atoms near a res- 
onance. These resonances are induced by g-wave bound 
states and are only strong enough to be observed due 
to the large second-order spin-orbit coupling of cesium 
atoms. The observed resonances are at low magnetic 
fields, which can be understood as a consequence of the 
large background scattering length in the aa collision 
channel; see Sec. III. B. II In additi on, /-wave Fe s hbach 
molecules have be en produced by ( Knoop et all l2008t 
iMark et aZll2007al) : their coupling to the s-wave contin- 
uum is too weak to lead to o bservable resonan ces in col- 
lision experiments. Finally, ( Lee et aTl . l2007t ) predicted 
a broad s-wave resonance at ~800 G in the aa channel. 



This is a magnetic field regime that has not been exper- 
imentally explored yet. 

Figure [21] shows the scattering len gth in the aa chan- 
nel. It has a zero crossing at 17.1 G (IChin et all l2004bl: 
iGustavsson et al ], [200i) and multiple narrow resonances 
below 50 G. The gradual change in a from —2500 ao to 
500 ao across 30 G in the Figure is actually the tail of 
a broad resonance with = —12 G; see Appendix and 
Table[IVl The negative Bq follows from fitting to Eq. 
Both the narrow and broad resonances have provided fa- 
vorable conditions for many exciting experiments. This 
included the attainment of BEG with Cs atoms, a Bose 
condensate ("Weber et all l2003cJ) the formation of Cso 
(jChiii _et_aL, 2003.; Herbig et all l2003l ). the observation 
of re sonances between ultracold molecules ( Chin et 'all . 
l2005h . and studies on Efimov physics; see Sec. IVI.GI 



C. Heteronuclear and other systems 

Most of the experimental and theoretical attention so 
far has been focused on locating and using magnetic 
Feshbach resonances in single-species alkali-metal atom 
gases. Over the last five years, however, considerable 
progress has also been made in locating Feshbach reso- 
nance in other atomic species and in mixtures of alkali- 
metal atoms. These systems are investigated for various 
reasons. Of particular interest is th e promise of more 
exotic quantum man y-body behavior (iBloch et all 120081 : 
iMenotti et d\ . 120081: iMicheh et all l2006l ). 

Heteronuclear systems provide the path to prepare 
mixtures of bosonic and fermionic quantum degener- 
ate gases. Intriguing applications include the cre- 
ation of fermionic molecu l es in an atomic Bq s e-Ferm i 
mixture (jOspelkaus et all l2006al: IZirbel et all l2008al ) 



and no vel quantum phases of fermions with unequal 
masses (|Petrov et aLl . l2007^ . Feshbach resonances pro- 



vide means to tune the interactions between different 
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species in order to explore quantum phases in various 
regimes. 

Atoms with magnetic moments interact via the long- 
range magnetic dipole-dipole interaction Vss in addition 
to the van der Waals and more short-range interactions. 
For alkali-metal atoms the effect of this dipole-dipole 
interaction on collective behavior is small. In atomic 
species with much larger magnetic moments, however, 
the dipole-dipole interaction can have a si gnificant im- 
pact o n the many-body behav i or of the gas (Goral all . 
I2OOOI: ISantos and Pfaul . l2006l : ISantos et al], l2003l ). In 
an atomic species like chromium (Sec. Illl.C.lj) . mag- 
netic Feshbach resonances can be used to tune the rel- 
ative strength of the short-ranged i nteractions and the 
long-range dipole- dipole interaction ( Lahave et aLl . [2007t 
lYi and Youl . l2002D . 

Another way to create exotic many-body systems is by 
pairing different atomic species into Feshbach molecules, 
which can then be converted into deeply bound molecules 
with a large electric dipole moment. Such a moment 
gives rise to large dipolar interactions, orders of mag- 
nitude larger than possible with magnetic dipole mo- 
ments. These molecules, which are either bosonic or 
fermionic, have many applications in dipolar mo l ecular 
quantum gases ( Biichler et al\. 2007 :lMicheli et all l2006f ) 



One of th ese unusual resonanc es was investigated in some 
detail in (jBeaufils et al\ . \200^ . 



and quantum computation (|DeMillel . l2002l ). 



1. Chromium 



In 2005 (jGriesmaier et all 120051 ) announced BEG of 
^^Cr atoms. This species has a magnetic moment that is 
six times l arger than that for alkali-metal atoms. Sub- 
sequently, ( Stuhler et all [2OO5I ) showed that the expan- 
sion of a cigar-shaped chromium condensate depended 
on the relative orientation of the magnetic moment and 
the elongated condensate, which for the first time showed 
the effects of the dipole-dipole interaction in a quantum- 

de generate gas. 

( Werner et all . I2005D measured fourteen magnetic 



Feshbach resonances in ^^Gr in the energetically lowest 
magnetic sublevel. Feshbach resonances were detected by 
measuring, after a fixed hold time, the number of remain- 
ing atoms as a function of magnetic field. Atom loss could 
only have occurred by enhanced three-body recombina- 
tion near the resonance. Figure[22] shows the .s-wave scat- 
tering length as a function of magnetic field derived from 
the observed locations of the Feshbach resonances and a 
multi-channel scattering model of the collision. Most of 
the resonances are g-wave resonances with the exception 
of one of the two nearly-degenerate resonances at 29 mT 
and those at 50 mT and 59 mT. These three resonances 
have c?-wave character. Moreover, (|Weruer et ai . 20051) 
also assigned two resonances as originating from d-wave 
collisions and coupling to an s-wave closed channel, i.e. 
£ = 2 and £c ~ according to the notation discussed in 
Sec. III. B. 41 Note that these resonances do not show up 
in the s-wave scattering length as displayed in Fig. [531 



2. Mixed species 

K+Rb - The first mixed system to receive a de- 
tailed effort to locate Feshbach resonances was a mixture 
fermionic *'^K and bosonic ^^Rb atoms ( Simoni et all 
I2OO3II . Based on experimentally determined inelastic and 
elastic rate coefficients at zero magnetic field they pre- 
dicted the location of fifteen Feshbach resonances with 
an uncertainty ranging from 10 G to 100 G. 

I lnouve et all 12004 ) presented the first direct deter- 



mination of the magnetic field location of three Feshbach 
resonances. The two atomic species, each in their ener- 
getically lowest hyperfine state, were optically trapped. 
As shown in Fig. [24] the Feshbach resonances were de- 
tected by measuring, after a fixed hold time, the number 
of remaining atoms as a function of magnetic field. The 
atom loss is due to three-body recombination. 

( Ferlaino et aLl . [2006l ) confirmed the positions of these 
three resonances and found nine additional resonances. 
Theoretical modeling uniquely assigned each resonance 
and determined the scattering lengths of the X^E+ and 
a'^E'*" Born-Oppenheimer potentials to about 2%. The 
difference between the experimental Feshbach locations 
and those of the best-fit theory was less than 1 G for all 
resonances. The experimental uncertainty in the rcso- 
na nce locations was 0.2 G. 

( Ferlaino et all |2006[ ) also predicted the location of 
resonances of several isotopic combinations. For the 
bosonic '^^K and ^^Rb collision with both atoms in their 
lowest hyperfine state ( Roati et all\200i\) confirmed the 
location of one such resonance. It was found at 317.9 G 



well wit hin the uncertaintie s quoted by ([Ferlaino et all 
I2OO6II . (iKlempt et all l2007^ observed a number of new 



Feshbach resonances, constructed an accurate potential 
model, and pre dicted resonances in other isotopic KRb 
combinations. ( Simoni et all [2008[ ) presented a refined 
near-threshold model for scattering and bound-state cal- 
culations for all isotopic combinations of K and Rb. 

Feshbach resonances in mixtures of and ®^Rb have 
been appli ed to the c reation of Feshbach molecules in op- 
tical traps (IZirbel et a i. 2008a. b ) and in single sites of an 
optical lattice ( Ospelkaiis et aLl . l2006a[ ): see Sec. IVI.BT] 
Recently these fermionic mol ecules have been t ransfe rred 
to more deeply bound levels ( Ospelkaus et aLl . l2008f ): see 
Sec. IV.B.4I Furthermore, interspecies interaction tun- 
ing has been exploited to study the collective behav- 
ior of a *°K+^^Rb Bose- Fermi mixture ( Ospelkaus et all 
l2006bD and to realize a tunable d ouble species BEG 



in a "^^K-|-^^Rb Bose-Bose mixture ([Thalhammer et 
I2OO8D . 



Li+Na - ([Stan et all \2004] have observed three mag- 
netic Feshbach resonances in the interaction between a 
degenerate fermionic ^Li gas and a Bose-Einstein con- 
densate of Na. The resonances were observed by detect- 



30 



115 
110 
, 105 
100 

95i 













J 


,,,,,, 

J) 




,,,,,,,,, 


7 














r 









5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 

Magnetic field [mT| 

FIG. 23 Scatterin g length of ^^C r atom s versus magnetic field. The feature near 29 mT is a pair of nearly-degenerate Feshbach 
resonances. From ( Werner et aLl [20051 ). 
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^"K in various hypcrfinc states. Their theor etical analy- 
sis 
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FIG. 24 Observation of Feshbach resonances in ^°K+*^Rb. 
The top panel shows in-trap absorption images of ''°K atoms 
after a fixed hold time of ~ 1 s at various magnetic fields. The 
label on each image gives the magnetic field in Gauss. No 
atoms could be seen at 542.2 G. The bottom panel shows the 
number of remaining atoms after the fixed hold time as 
a function of magnetic field. Narrow fea tures are observed a t 
492 G, 512 G, and 543 G. Adapted from (|lnouve et am2004D . 



ing atom loss when sweeping the magnetic field at con- 
stant rate through the resonances. Atom loss could only 
have occurred from three-body recombination or from 
molecule formation during the sweep. The observed res- 
onances at 746.0 G, 759.6 G, and 795.6 G are s-wave 
resonances. In a recent theoretical study ba sed on this 
experimental input data, ( Gacesa et aZ.I . [2008[ ) derive pre- 
cise values of the triplet and singlet scattering lengths for 
both the ^Li-Na and the ''Li-Na combination. Moreover 
they predict a variety of additional Feshbach resonances 
within an experimentally attainable field range. 

6Li-h^°K - (jWille et all l2008t) described the observa- 
tion of thirteen Feshbach resonances in fermionic ^Li and 



, w hich relies on the model developed by (jStan et al\ . 
l2004l) and discussed in Sec. III. C. 21 indicated that the 
resonances were either s- or p-wave resonances. This iso- 
topic combination is a prime candidate for the study of 
strongly interacting Fermi-Fermi mixtures. 

Na-|-Rb - Predictions of Feshbach resonance loca- 
tions based on analysis of high-resolution Fourier spec- 
troscopy of the molecular X^I]+ and a^E+ states in 
a 600 K beam of NaRb molec ules are descr i bed i n 
(jBhattacharva et all. \2004) and (jPashov et all |2005[) . 
For example ( Pashov et all l2005l ) predict for the ultra- 



cold collision between ^'^Na and ®^Rb, both in the en- 
ergetically lowest hyperfine state, s-wave Feshbach reso- 
nances at 170 G and 430 G with an uncertainty of about 
50 G. This uncertainty is sufficiently small that the pre- 
dictions will be helpful for planning experiments which 
can accurately locate the resonances. 

Initial experiments on other combinations of mixed 
atomic species have been performed. Feshbach rcs- 
ona nces have been re ported in the ^Li-|-®^Rb system 
by (iDeh et all I2008D . in the '^Li-H^'^Rb system by 
(iMarzok e^ aLl. 120091) . and in the ^"^Rb+^^^Gs system by 
Jpilch et fflZ.l . l2009D . 



3. Isotopic mixtures 

A special case of a mixed system is that where different 
isotopes of the same element are combined. In particular, 
isotopic mixtures of Rb, K, and Li have been studied. 
Isotopic mixtures of K or Li are of particular interest as 
both fermionic and bosonic isotopes exist. 

Feshbach resonances in isot opic mixtures of rubidium 
have recently been observed. ( Papp and Wiemai] . |2006| ) 
found two s-wave Feshbach resonances in the collision 
of ^■''Rb and *^Rb when both isotopes are in their low- 
est hyperfine state. Their magnetic field locations of 
265.44(0.15) G and 372.4(1 . 3) G are consistent with the 
pr edictions of ([Burke all Il998l ) . 

(jvan Kempen et all . 2004 ) have predicted the location 
of ^Li-h'^Li Feshbach resonances. When both isotopes 
arc in the lowest hyperfine state resonances occur be- 
tween 200 G and 250 G as well as between 550 G and 
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560 G. Four of the se resonances have been observed by 
((Zhan^irill,[200l). 

In predictions for isotopic mixtures "mass scaling" 
is often used. As the interatomic potentials arc to 
good approximation independent of the isotopic com- 
position of the dimer, the only thing that changes is 
the (reduced) mass of the dimer. Corrections are due 
to the bre akdown of the Bo r n-Oppenheimer a p proxi - 
mation. (jvan Kempen et all . l2002l : ISeto et cd. 1, (20001) 
showed that for Rb the experimental data on ®^Rb*^Rb, 
^^Rb^^Rb and ^^Rb^^Rb are consistent with mass scal- 
ing . A new analysis that inc ludes the 2006 observations 
of (|Papp and Wiemanl . [20061 ) is needed. 

For the lithium system ( van Kempen et all . |2004 ) 
showed that mass scaling is insufficient to explain the 
observed data for the ho monuclear ^Li-|-^Li and ^ Li-|-^Li 
systems. Consequently, ( van Kempen et all . 120041 ) quote 
a 1 G uncertainty for the location of ^Li-|-^Li Feshbach 
resonances from the breakdown of mass scaling. 



IV. CONTROL OF ATOMIC QUANTUM GASES 

Tuning two-body interactions via Feshbach resonances 
is the experimental key to control collective phenomena 
in degenerate quantum gases. This has found numer- 
ous ap plications, both with a t omic Bose- Eins t ein co nden- 
sates (ICornell and Wiemanl . l2002t iKetterld . |2002| ) and 

1,1 



1. Attainment of BEG 

Some atomic species offer favorable collision properties 
for the attainment of BEC without any necessity of inter- 
action t uning; ^^Rb and ^'^Na are the most prominen t ex- 



with degenerate Fermi gases (jlnguscio et 



The different decay properties of Bose and Fermi gases 
near Feshbach resonances play a crucial role for the ex- 
periments. For Bose gases, resonant two-body scattering 
in general leads to rapid decay via three-body collisions, 
as we have discussed in context with loss spectroscopy 
on Feshbach resonances in Sec. IIII.A.2I Three-body de- 
cay limits the practical applicability of Feshbach tuning 
to Bose gases, restricting the experiments to the dilute 
gas regime where the scattering length is small compared 
to typical interparticle separations. In contrast, Fermi 
gases can b e remarkably stabl e near s-wave Feshbach 
resonances (IPetrov et "aI]. l2004D . 

For atomic Bose- Einstein condensates in the dilute gas 
regime (Sec. lIV.Ap . the collective behavior can be de- 
scribed in a mean-field approach. In strongly interacting 
Fermi gases (Sec. IIV.B[) . the role of the scattering length 
is much more complex. Here Feshbach tuning can be 
used to control the nature of fermionic pairing in differ- 
ent supcrfluid regimes. 



A. Bose-Einstein condensates 

In experiments on atomic BEC, the role of Feshbach 
tuning can be divided into two parts. Firstly, the control 
of collision properties can be essential for the attainment 
of BEC. Secondly, the possibility to control the mean- 
field interaction opens up a variety of interesting appli- 
cations. 



amples (|Cornell and Wiemanl . I200I iKetterlj . |2002[) . In 
other cases, however, Feshbach tuning is essential either 
to produce large condensates (^Li), or to achieve BEC at 
ah (85Rb, i33cs, 39K). 

Let us first consider the general question, what is a 
"good" scattering length a for making a BEC. First, 
a should be positive, because condensates undergo col- 
lapse at negative scattering length when the number 
of con densed atoms excee ds a relatively small critical 
value ( Bradley et al\ . [19971 ) . Moreover, a should not be 
too small, because a sufficiently large elastic collision 
rate is required for evaporative cooling towards BEC; 
the cross section for elastic collisions between identical 
bosons is Sttc^. Finally, the scattering length should not 
be to o large to avo i d rapid decay by three-b ody colli- 
sions (|Roberts et aLl . l2000Hv\^ber etal. '2003b') as three- 
body decay scales as ( Braaten and Hammer. . ,2006t 
iFedichev et all Il996bl ) . In practice these conditions re- 
sult in typical values for a "good" value of a between 
a few ten and a few hundred times the Bohr radius ag. 
In detail, the optimum value for a depends on the con- 
finement properties of the trap and behavior of inelastic 
decay 

For ^Li, early magnetic trapping experiments showed 
BEC in the internal state / = 2,m = 2, where the 

scattering length a = — 27ao is small and negative 

(iBradlev et aLl . Il997l [1991 : here the number of conden- 
sate atoms was limited throug h collapse to only a few 
hundred. Later exp eriments ( Khavkovich et all . l2002t 
IStrecker et all |2002| ) on optically trapped ''Li in the in- 
ternal state / = l,m = 1 exploited the 737-G Fesh- 
bach resonance to tunc the scattering length from its 
very small background value (abg ~ -1-5 oq) to sufficiently 
large positive values, typically in a range between -1-40 oq 
and -1-200 flQ. Evaporative cooling then resulted in con- 
densates with up to 3 x 10^ atoms. More recently, 
( Gross and Khavkovichl [20081 ) exploited Feshbach tuning 
in the state / = 1, 771 = for the all-optical production of 
a BEC. They obtained favorable conditions for efficient 
evaporative cooling at 866 G, where the scattering length 
is about -1-300 ao. 

Bose-Einstein condensation of ®^Rb w as achieved by 
evapo rative cooling in a magnetic trap ( Cornish et all 
l2000f ) exploiting the broad resonance at 155 G 
(Sec. IIII.B.4|) in the state f = 2,m = -2 to tune 
the scattering length to positive values. The large 
negative background scattering length abg = —443 ao 
would limit the number of condensate atoms to less than 
one hundred. Two stages of cooling were performed. 
The first stage used a magnetic field of 250 G, where 
the scattering length a is close to its background value; 
the second was close to the resonance at 162.3 G, where 
a = -1-170 ao. This procedure optimized the ratio of 
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clastic to inelastic collision rates ( Roberts et I2OOOI) 
for the temperatures occurring during these stages. 
Fcshbach tuning played a cr ucial role for the at- 



tainment of BEG in ^^^Cs (|Kraemer et all l2004t 
IWeber et all l2003aD . The condensate was produced in 
an optical trap in the state / = 3, m = 3. In this lowest 
internal state, two-body decay is energetically forbidden, 
and the scattering length a{B) shows a large variation 
at low magnetic fields (Fig.[22|), which is due to a broad 
Feshbach resonance at — 12G. In a first evaporative cool- 
ing stage a shallow large-volume optical trap was em- 
ployed, and a large scattering length of a = +1200 ao at 
B = 73 G provided a sufficiently large elastic collisions 
rate at rather low atomic densities. The second cool- 
ing stage employed a much tighter trap. Here, at much 
higher densities, an optimum magnetic field of 21 G was 
found, where a = +210 cq. Highly efficient evaporation 
led to the attainment of BEG. Later experiments revealed 
a mi nimum of three-body d ecay in this magnetic field re- 
gion ( Kraemer et al ]. l200d) : see also Sec.lVI.G.2l 

For the attainment of BEG in ^^K, (|Roati et aLl . l2007D 
employed a mixture with ®^Rb atoms with both species 
being in their internal ground states. First, an inter- 
species Fcshbach resonance near 318 G was used to op- 
timize sympathetic cooling; the interspecies scattering 
length was tuned to +150 ao by choosing a magnetic field 
of 316 G. Then, after the removal of the Rb atoms, final 
evaporative cooling towards BEG was performed near the 
402-G resonance of ^^K fSec. lIII.B.3)) with the scattering 
length tuned to a positive value of +180 oq. 



2. Condensate mean field 

Trapped atomic B EGs in the d i lute-g as regime are 
commonly described ( Dalfovo et aLl . [l999l ) by the Gross- 
Pitaevskii equation for the condensate wave function <&, 



2y72 



2 m 



(57) 



where T4xt is the external trapping potential. Interac- 
tions are taken into account by the mean-field potential 



Klf = : 



(58) 



where the atomic number density n is related to $ by 
n = This mean field enters the Gross-Pitaevskii 

equation as a nonlinearity and leads to many interesting 
phenomena. 

In the Thomas- Fermi regime of large condensates with 
a > one can neglect the kinetic energy term and obtain 
the equilibrium density distribution of a BEG 



(/i " Vo^t) 



(59) 



which applies for /i > Voxt! otherwise n — Q. For a given 
particle number iV, the chemical potential /i follows from 
the normalization condition N = J ndPr. 
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FIG. 25 In-situ images and density profiles showing the vari- 
ation of the equilibrium size of a magnetically trapped *^Rb 
BEG close to a Feshbach resonance. With decreasing mag- 
netic field (a-e), the scattering length increases from a very 
small positive value at 165.2 G (zero crossing at 165.8 G) to a 
very large value of a « +6000 ap at 156.4 G (resonance posi- 
tion at 155.0 G). From (|Gornish et aZ.l.l2000l ). 



In many cases of practical relevance, stable conden- 
sates with positive a are confined in harmonic traps and 
are in the Thomas-Fermi regime. The condensate is then 
characterized by the Thomas- Fermi radius rxF, given as 
the radius at which the external trapping potential equals 
the chemical potential, and the peak number density 7io- 
These two quantities follow the scaling laws 
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Figure [25l illustrates how the size of a trapped ®^Rb con- 
densate increases when the Feshbach resonance at 155 G 
is approached. The in-situ mea surements were use d 
to experimentally determine a{B) ( Gornish et flZ.I . [2000f ) . 
The results are in good agreement with later measure- 
ments of the molecular binding energy, which allowed for 
a more precise determination of the scattering prop erties 
near the Feshbach resonance (jGlaussen et aLU2003l ). 

The mean-field approach is valid for scattering lengths 
which are small compared to the typical interparticle 
distance. The prospect to observe beyond- mean-field 
effects in BEGs (|Dalfovo et all Il999l ) has been an im- 
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portant motivation for experiments near Feshbach res- 
onances at large a. In atomic Bose gases, however, the 
fast inelastic d ecay makes i t very d ifficult to observe such 
phenomena. ( Papp et all . l2008l ) finally demonstrated 
beyond-mean-field behavior by Bragg spectroscopy on a 
^^Rb BEC. In molecular BECs created in atomic Fermi 
gases (Sec. lIV.B.ip . the collisional stability facilitates 
the observation of beyo nd-mean-field behavior by simpler 
means. For example, ( Altmever et all . |2007( ) observed 
beyond-mean-field behavior by studying collective oscil- 
lations of a ^Lio molecular BEC. 



3. Controlled collapse and bright solitons 

For negative scattering lengths the condensate mean 
field is attractive. The resulting nonlincarity can then 
lead to a condensate collapse and to the formation of 
bright matter-wave solitons. To study phenomena of this 
kind by Feshbach tuning, the general experimental strat- 
egy is to first produce a stable BEC at positive a. Then, 
the attractive interaction is introduced by changing a 
to negative values. E xperiments of this c l ass have been 



performe d with ^^Rb (ICornish e t ~aP.. '2006': 'Ponlc v et al. 



Weber et all . l2003aD 



bj 



2001: Roberts e.t allyMm . '^Li (.Khavkovich et al. 



Strecker et all.l2002l) . and "^Cs (jRvchtarik et al 



2002 : 



20041: 



Exploiting the 15 5-G Feshbach resonance in ^^Rb, 
( Roberts et aZ.I . [200l[ ) investigated the stability of a BEC 
with attractive interactions. They first produced the 
condensate in a magnetic trap at a moderate positive 
scattering length. They then slowly changed the atom- 
atom interaction from repulsion to attraction by ramping 
the magnetic field into the region of negative scattering 
length. With increasing attractive interaction they ob- 
served an abrupt transition in which atoms were ejected 
from the condensate. These measurements of the onset 
of condensate collapse provided a quantitative test of the 
stability criterion for a BEC with attractive interactions. 

The controlled collapse in ®^Rb following a sudden 
change of a led to the spectacular observation of a 
"Bosenova," a condensate i mplosion wit h fasc inating 
and unexpected properties ( Donley et all . l200lf ). An 
anisotropic burst of atoms was observed that exploded 
from the condensate during the early stage of collapse 
(Fig. [26|). leaving behind a highly excited long-lived rem- 
nant condensate. Strikingly, the number of atoms in 
the remnant BEC was significantly larger than the crit- 
ical number for a collapse. The surprising fact that 
the remnant BEC did not undergo further collapse was 
later explained by its fragmentation into bright solitons 
(ICornish et a^.l ■ l2006^ . 

Condensate collapse experiments were also used to de- 
tect the presenc e of a small BEC of Cs atoms in an opti- 
cal surface trap ( Rvchtarik et aZ.l . l2004[ ). While a thermal 
gas did not show loss when the scattering length was sud- 
denly switched to negative values, the condensate showed 
up in the sudden onset of collapse-induced loss. 



o 9 




'AWL 




FIG. 26 Striking phenomena have been observed in the 
controlled collapse of a ^'^Rb BEC. The images show the 
formation of 'jets', where streams of atoms with highly 
anisotropic velocities are ejected by the collapsing conden- 
sate. Reprin ted by permission from Macmillan Publishers 
Ltd: Nature (|Donlev et aZ.I . l200ll '). copyright 2001. 




FIG. 27 A train of matter- wave s olitons created from an op- 
tically trapped BEC of '^Li atoms (jStrecker et a/.l . l2002l \ The 
individual solitons contain up to about 5000 atoms. Figure 
courtesy of Randall Hulet. 



Bright solitons were observed in experiments on 
^Li atoms near the bro ad 736-G Feshbach resonance. 
(|Khavkovich et all . l2002f ) produced a BEC by evapora- 
tive cooling at a w 40 ao in an optical trap. They then 
released the BEC into a one-dimensional optical waveg- 
uide and studied the propagation of the resulting matter- 
wave packet for an ideal gas (a = 0) and a gas with a 
small attractive mean-field interaction (a = — 4ao). In 
the latter case, they observed the dispersion-free prop- 
agation that is characteristic for a soliton. In a similar 
experiment ( Strecker et aZ.! . [2002h created a train of soli- 
tons, see Fig. [27l from an optically trapped ^Li BEC by 
abruptly switching the scattering length from 200 ao to 
— 3 oq. They also observed the propagation of the solitons 
in the trap and their mutual repulsion. These spectacu- 
lar experiments on bright solitons highlight the analogy 
between bright matter-wave solitons and optical solitons 
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FIG. 28 Density profiles of a Cs BEC after 100 ms of expan- 
sion at (a) a = 210 ao and (b) a = 0. The expansion energy of 
the non-interacting co ndensate is as low as ks x 50 pK. From 
(jKraemer et all\2004 ). 



in fibers and thus the intimate connection between atom 
optics with BECs and light optics. 



4. Non-interacting condensates 

The zero crossing of the scattering length near a Fesh- 
bach resonance can be used to realize non-interacting 
ideal-gas condensates. BECs of '''Li, ^^K, ^^Rb, and ^^^Cs 
are good candidates to reduce |a| to very small values on 
the order of oq or smaller. 

To explore non-interacting condensate s with ^^^'^Cs 
( Kraemer et "aIl. l2004lWeber aU[2003aD exploited the 
zero crossing near 17 G and studied expansion after re- 
lease from the trap. The scattering length was abruptly 
changed to zero when the condensate was released from 
the trap. An extremely slow expansion was observed with 
a release energy as low as x 50 pK; see Fig. [551 The 
surprising observation that the release energy is a factor 
of five below the kinetic energy associated with the mo- 
tional ground state of the trap is explained by the fact 
that the initial size of the expanding matter-wave packet 
is determined by the repulsive condensate self-interaction 
before release, which is larger than the bare ground state 
of the trap. The momentum spread is thus significantly 
smaller. In contrast, a slow change of a to zero before re- 
lease would have ideally resulted in a wavepacket with po- 
sition and momentum spread corresponding to the bare 
ground state. 

Besides the small equilibrium size of the condensate, 



a vanishing scattering length has profound consequences 
for the collective behavior of a BEC. The sound veloc- 
ity (cx na^/^) is vanishingly small so that all excitations 
will become particle-like, and not phonon-like. Moreover, 
the healing length (cx na~^^^) becomes very large, which 
may be applied in exp eriments on rotating condensates 
( Madison et to increase the core size of vor- 



tices. 

Non-interacting condensates are promising for the ob- 
servation that are masked by interacting effects, as e.g. 
phenomena on a lower energy scale. In atom interfcrom- 
etry the mean-field interaction of a condensate is asub- 
stantial systematic error source, as ( Gupta et all . 1200^ 
observed in t he con text of photon recoil measurements. 
( Roati et al 

I l2004^ have studied Bloch oscillations in 
an optical lattice under the influence of gravity. They 
showed that, for interacting bosons, the oscillations lost 
contrast much faster than for identical fcrmions without 
s-wave interaction. A non-interacting BEC combines the 
advantages of an ultralow momentum spread with very 
long observation times. Two recent experiments have 
report ed on long-lived B l och o scillations with BECs of 
i33Cs (iGustavsson et all . l2008l ) and ^^K (|Fattori et all . 
l2008aD . which is an important advance towards high- 
precision atom interferometry. This could, for example, 
open up new possibi lities for precisiq i i mea surements of 
gravitational effects ( Carusotto et aZ.I . l2005f ). 

Another intriguing application of the zer o crossing of a 
Feshb ach resonance was demonstrated by ( Lahave et all . 
l2007ll with a BEC of ^^Cr atoms. This species exhibits a 
very large magnetic dipole-dipolc interaction because of 
its magnetic moment of 6 /is . When the isotropic con- 
tact interaction is reduced by tuning the scattering length 
close to zero, the magnetic dipole interaction dominates. 
In ^^Cr this was achieved near the 589-G Feshbach res- 
onance (Fig. I23p . The resulting dipolar quantum gas 
represents a model system for a "quantum ferrofluid" , 
the anisotropic propertie s of which have been attract- 
ing considerable interest ( Menotti et all [20081 ). In fur- 
ther wo rk on ^^Cr BEC nea r the zero of the scattering 
length, ( Lahave et all . [2008h investigated the controlled 
collapse of th e system and dem onstrated its complex dy- 
namics and ( Koch et all |2008| ) studied the stability of 
the dipolar condensate depending on the trap geome- 
try. The effect of the magnetic dipole interaction has 
also been observed in non-interacting condensates made 
of alkali atoms wi th magnetic moment s of the order of 
ItiR. for bv (iFattori et all l2008bh and for ^Li by 
( Pollack et all |2009| ) . The recent observation of Ander- 
son localizati on of matter wave s in a disordered optical 
potential by ( Roati et all l2008l ) represents a further ex- 
iting application of a non-interacting condensate. 



B. Degenerate Fermi gases 

In experim e nts on ultracold Fermi 
( Inguscio et all l2008l ). Feshbach resonances 



gases 
serve 
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as a key to explore many-body physic s in the strongly 
interacting regime (jBloch et This regime is 



realized when the scattering length exceeds the interpar- 
ticle spacing and connects the field of ultracold atoms 
to fundamental questions in various fields of physics, 
like high Tc-superconductors, nuclear matter, neutron 
stars, and the quark-gluon plasma. The f irst Feshbach 
reson a nce in a Fermi gas wa s observed by ( Loftus et all . 
I2OO2I ). (jO'Hara et aLl . l2002a[) produced the first strongly 
interacting Fermi gas. Since then the research field 
has undergon e rapid developments with many exciting 
achievements ( Inguscio et aZ.I . [20081 ) . 

The decay properties of ultracold Fermi gases 
are strongly influ er iced by Pauli 's exclusion p ri nciple 
(lEsrv et all I2OOII : iPetrovi l2003l : iPetrov et all l2004t 



ISuno et alll2003 ). Three-body recombination processes 
in one- and two-component Fermi gases necessarily in- 
volve identical particles. This generally leads to a sup- 
pression of loss as compared to Bose gases or systems 
with three nonidentical particl es. The majori t y of r ecent 
experiments on Fermi gases ( Inguscio et all T2OO8I ) has 
focussed on two-component spin mixtures of ^Li or 
with resonant s-wave interactions, realized near broad, 
entrance-channel dominated resonances. In such sys- 
tems, it is possible to realize a resonant s-wave inter- 
action (a — *■ ±00) practically without any decay. Never- 
theless, t hese resonances are acco mpanied by subtle loss 
features ( Dieckmann et aLl . l2002l ). which do not appear 
at the resonance center, but at side where the scattering 
length is positive. In contrast to the remarkable stabil- 
ity near this special s-wave scenario, three-body colli- 
sions near a p- wave Feshbach re sonanc e usually lead to 
significant loss dReeal et aLl . l2003 b: Sch unck et aZ.I . [2005t 
IZhang et a^.l ■ l2004^ . 

Here, as prominent examples for the application of 
Feshbach tuning, we review the attainment of EEC of 
molecules (Sec. lIV.B.ip and studies of the BEC-BCS 
crossover and the observation of fcrmion superfluidity 
rSec. lIV.B.2l . 



1. BEC of molecules 



Bose-Einstein condensation of molecules were cre- 



ated in atomic Fermi gases of ^Li (iBourdel et all [200I 



l^ermi ga 

IJochim et all. l2003bl : IZwierlein et all l2003l )~Mid '^°K 
( Greiner et oZ] , |2003[ ). The molecules are very weakly 
bound dimers at the side of an entrance-channel dom- 
inated s-wave Feshbach resonance where the scattering 
length is positive and very large. These dimers are 
formed in a halo state (Sec. IV.B.2|) . which is stable 
against inel astic decay in ato m-dimer and dimcr-dimcr 
collisions ((Pctrov et all [2003 l . This stability originates 
from basically the same Pauli suppression effect that also 
affects three-body decay in an atomic Fermi gas. Such 
weakly bound molecules can be detected by converting 
them back to atoms or by direct absorption imaging; see 
Sec. lV.B.ll 
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FIG. 29 Emergence of a molecular BEC in an ultracold Fermi 
gas of atoms, observed in time-of-flight absorption im- 
ages. The density distribution on the left-hand side (upper 
graph, 2D surface plot; lower graph, ID cross section) was 
taken for a weakly interacting Fermi gas which was cooled 
down to 19% of the Fermi temperature. After ramping across 
the Feshbach resonance no BEC was observed as the sample 
was too hot. The density distribution on the right-hand side 
was observed for a colder sample at 6% of the Fermi tempera- 
ture. Here the ramp across the Feshbach resonance resulted in 
a bimodal distribution, revealing the presence of a molecular 
BEC with a condensate fraction of 12%. Repri nted by permis- 
sion from Macmillan Publishers Ltd; Nature (jCreiner et all 
I2OO3I ). copyright 2003. 



In a spin mixture of ^Li in the lowest two inter- 
nal stat es, the route to mol ecular BEC is particularly 
simple ( Jochim et all l2003lJ l. Evaporative cooling to- 
wards BEC can be performed in an optical dipolc trap 
at a constant magnetic field of about 764 G near the 
broad 834-G Feshbach resonance; here a = -|-450Gao 
and Ef) = ks X 1.5 ^K. In the initial stage of evapora- 
tive cooling the gas is purely atomic and a is the rele- 
vant scattering length for elastic collisions between the 
atoms in different spin states. With decreasing tempera- 
ture the atom- molecule equilibrium (Sec. IV.A.3P favors 
the formation of molecules and, in the final evapora- 
tion stage, a purely molecular sample is cooled down to 
BEC. The large atom-dimer and dimer-dimer scattering 
lengths of 1.2 a and 0.6 a along with strongly suppressed 
loss (Sec. IV.B.3I) facilitate an efficient evaporation pro- 
cess. In this way, molecular BECs are achieved with a 
condensate fraction exceeding 90%. 

The experiments in '^'^K followed a different approach 
to achieve molecular BEC (|Greiner et aZ.l . l2003[ ). For "Or 
the weakly bound dimers are less stable because of less 
favorable short-range three-body interaction properties. 
Therefore the sample is first cooled above the 202-G Fesh- 
bach resonance, where a is large and negative, to achieve 
a deeply degenerate atomic Fermi gas. A sweep across 
the Feshbach resonance then converts the sample into a 
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partially condensed cloud of molecules. Figure[29ldemon- 
strates the emergence of the molecular BEC in *°K. 

The molecular BEC can be described in the mean- 
field approach outlined in Sec. IIV.A.2I by simply replac- 
ing the atomic with the molecular mass (m — s- 2m), 
and the atomic with the molecular scattering length 
(a 0.6a). The mean field of the molecular conden- 



sate was experimentally stu died in (jBartenstein et all . 
l2004bl : iBourdel et al\ . l2004f ). However, because of the 
large scattering length, molecu lar BECs show consider- 
able beyond-mean-field effects ( Altmever et a/.l , [20071 ) . 



2. BEC-BCS crossover and fermion superfluidity 

At a Fcshbach resonance in a two-component Fermi 
gas, different regimes of fermion pairing and superfluid- 
ity can be experimentally realized. Pairing on the side 
with a > can be understood in terms of molecule for- 
mation, and superfluidity results from molecular Bose- 
Einstein condensation. On the other side of the resonance 
(a < 0), pairing is a many-body effect and the ground 
state of the system at zero temperature is a fermionic su- 
perfluid. In the limit of weak interactions, this regime can 
be understood in the framework of the well-established 
Bardeen-Cooper-Schrieffer (BCS) theory, developed in 
the 1950s to describe superconductivity. Both limits, 
BEC and BCS, are smoothly connected by a crossover 
through a regime where the gas is strongly interacting. 
This BEC-BCS crossover has attracted considerable at- 
tention in many-body quantum physics for more than two 
decades, and has recently reviewed in (IChen et all . l2005t 
iGiorgini et "d] . l2008l : llnguscio et "d] . [20081). A theoreti- 
cal description of this challenging problem is very diffi- 
cult and various approaches have been developed. With 
tunable Fermi gases, a unique testing ground has become 
available to quantitatively investigate the crossover prob- 
lem. 

The interaction regime can be characterized by a di- 
mensionless parameter l/{kpa), where kp is the Fermi 
wave number of a non-interacting gas, related to the 
Fermi energy hyEp = h'^kj,/{2m). For l/{kpa) > 1, the 
molecular BEC regime is realized. For l/{kpa) ^ —1, 
the system is in the BCS regime. In between (l/l/cpal < 
1), the Fermi gas is strongly interacting. In the experi- 
ments, Ep/kB = 1 /iK gives a typical value for the Fermi 
energy and 1/kp « 4000 ao sets the typical length scale. 
The realization of a strongly interacting gas thus requires 
|a| ^ 4000 ao, which for the particularly broad ^Li reso- 
nance (Fig. fTU]) is obtained over a more than 100-G wide 
magnetic field range. 

A particularly interesting situation is the exact reso- 
nance case, where l/{kpa) ~ 0. Here, a is no longer 
a relevant quantity to describe the problem and scat- 
tering is only limited by unitarity (Sec. III.A.2[) . Conse- 
quently, kp^ and Ep define the relevant scales for length 
and en ergy, and the Fermi gas acquires un i versal prop- 
erties (|Giorgini et all l2008l : llnguscio et all l2008l ). For 



example, the size and shape of a harmonically trapped 
"universal Fermi gas" can be obtained just as a rescaled 
version of a non-interacting Fermi gas. 

Experimentally, various properties of strongly inter- 
acting Fer mi gases have been ex plored in the BEC-BCS 
crossover ( Inguscio et all [20081 ) . All these experiments 
have been performed on two-component spin mixtures of 
^Li near the 834-G resonance or of ^"K near the 202-G 
resonance. In both cases, the resonances have entrance- 
channel dominated character, where the two-body inter- 
action can be modeled in terms of a single scattering 
channel and universality applies fSec. III.B.2p . This con- 
dition is particularly well fulfilled for the ^Li resonance, 
where the resonance is exceptionally strong. 

I Regal et all [2004[ ) introduced fast magnetic-field 
sweeps to observe the condensed fraction of pairs in the 
crossover. Starting with an ultracold '"'K Fermi gas 
in the strongly interacting regime, they performed fast 
Fcshbach ramps into the BEC regime. The ramps were 
fast compared to the time scale of establishing a thermal 
atom-molecule equilibrium by collisions; see Sec. IV. A. 31 
However, the Fcshbach ramps were slow enough to adi- 
abatically convert fermion pairs formed in the strongly 
interacting regime into molecules. After the ramp, the 
observed molecular condensate reflected the fermion con- 
densate b efore the ramp. The f ast-ramp method was ap- 
plied by ([Zwierlein et all [2004[ ) to observe fermion con- 
densates in ^Li. 

For ^Li ([Bartenstein et all [2004bf) showed that slow 
Fcshbach ramps allow conversion of the gas in a reversible 
way from the molecular BEC to the BCS regime. Here 
the gas adiabatically follows and stays in thermal equilib- 
rium. They also observed in-situ profiles of the trapped, 
strongly interacting gas and measured its changing size 
for variable interaction strength. 

Collective modes in the BEC-BCS cro ssover were stud- 
ied in ^Li gases. ( Kinast et all [2004[ ) reported on ul- 
traslow damping in a universal Fermi gas with rcso- 
nant interactions, provid ing evidence for superfluidity. 
( Bartenstein et al\ . [2004a[ ) measured how the frequencies 
of collective modes in the crossover changed with vari- 
able interaction parameter \/{kpa). They also observed 
a breakdown of hydrodynamic behavior on the BCS side 
of the resonance, which marks a transition from the su- 
pcrfluid to the normal phase. Precision measurements of 
collective modes also revealed beyond-mean-field eff ects 
in the molecular BEC regime ( Altmever et a/.l . [2007t ). 

( Chin et "aZ[ . [2004a[ ) performed spectroscopy on 
fermion pairs by using a radio-frequency technique. 
They measured the binding energy of the pairs in the 
crossover. They showed how an effective pairing gap 
continuously evolved from the molecular regime, where 
it simply reflects the dimer binding energy, into a 
ma ny-body regime of pa iring; see also more recent work 
by ([Schunck eV aL[ [20081) 

( Partridge efaZI . [2005t) used an optical molecular- 
probe technique based on photoassociation to measure 
the closed-channel contribution of pairs in the crossover. 
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FIG. 30 Measurements of the closed-channel fraction Z for 
pairs in the BEC-BCS crossover. The experiment uses the 
834-G resonance in ®Li and a photoassociation probe that 
only couples to closed-channel singlet component. The verti- 
cal dashed lines indicate the boundaries of the strongly inter- 
acting regime, kF\a\ > 1. The dotted line shows the result of a 
coupled channels calculation for molecules; see Sec. lII.C~5l and 
Fig. 1161 Comparison with the experimental data shows that 
two-body physics describes the situation well up to close to 
the 834-G resonance. For higher fields, Z shows strong many- 
body effects. Above the resonance, where two-body pairs can- 
not exist, the many-body system shows th e closed-channel 
admi xture of the many-body pairs. From (jPartridee et all 
120051 '). 



Their results, displayed in Fig.[30l show that this fraction 
is very small in the strongly interacting regime. These 
observations strongly support single-channel descriptions 
for the crossover along with the concept of universal- 
ity. The results also demonstrate that fermionic pairing 
reaches from the strongly interacting regime well into the 
weakly interacting BCS regime. 

Superfluidity of a ^Li Fermi gas in the BE C-BCS 
crossover was observed by ( Zwierlein et al\ . [200^ . They 
produced a rotating Fermi gas and observed the forma- 
tion of vortex arrays. Here Feshbach tuning was applied 
not only to explore different regimes in the crossover, but 
also to increase the vortex cores in the expanding Fermi 
gas after release from the trap; the latter was essential 
to observe the vortices by optical imaging. 

Currently, there is considerable interest in ex- 
ploring novel regime s of pairing and superflu idity. 
(jZwierlein eFaZI . l2006h and (jPartridge et all Hooi) per- 
formed experiments with unbalanced mixtures of two 
spin states, i.e. polarized Fermi gases. This led to 



deeper insight into p h enomena like phas e separation 
(jPartridge et all l2006l ; 111 Shin et all l2008h . Ultracold 
Fermi-Fermi mixtures of different speci es, like ^Li and 
"^°K, have recentl y bec ome availabl e (iTaglieber et all . 
120081 : IVoigt et fflZ.I .' i2008l : IWiUe et all 12008D, adding the 
mass ratio and independent control of the trapping po- 
tentials for both components to our tool box to explore 



the broad physics of strongly interacting fermions. 



V. ULTRACOLD FESHBACH MOLECULES 

Cold molecules are at the center of a rapidly develop- 
ing re search field ( Dovle et~al\ . |2004 iHutson and Solddiil . 
l2006l) . offering many new opportunities for cold chem- 
istry, precision measurements, many-body physics, and 
quantum information. The coldest attainable molecules, 
at temperatures in the nanokelvin range, are created by 
Feshbach association in ultracold atomic gases. Here 
Feshbach resonances serve as the experimental key to 
couple pairs of colliding atoms into molecules. 

In 2002, (|Donlev et all l2002[ ) observed coherent oscil- 
lations between atom pairs and Feshbach molecules in a 
BEC of ^^Rb atoms. The oscillation frequency reflected 
the small binding energy of the dimer in a weakly bound 
state and provided indirect evidence for the creation of 
molecules. In 2003, several groups reported on more di- 
rect observations of Fesh bach molecvdes i n Fermi gases of 
■^°K (iReeal et fl/l.l2003aD and ^Li (.C ubiz oUes et aLl . 120031: 
IJochim et aLl. l2003at IStrecker et all. l2003h and BE Cs of 



Cs (iHerbig et all\20o£ . ^^Rb (iPgrr et a]] ■ l2004a^ . and 
Na ( Xu et aU 12003 ). This rapid development culminated 
in the attainment of molecular BEC at the end of 2003 
(Sec. lIV.B.ip . and has paved the way for numerous ap- 
plications. 

A comprehensive review of Feshbach molecules 
and their theoretic al background has been given by 
(|K6hieZeIai],[200i). We do not consider optical meth- 



ods of making molecul es, since this was recently reviewed 
by ( Jones et all l2006l ). Here we discuss various forma- 
tion methods based on magnetic Feshbach resonances 
fSec. lV.A]) and the main properties of Feshbach molecules 



A. Formation 

Various schemes to create ultracold molecules near 
Feshbach resonances have been developed in the last few 
years, most of them relying on the application of time- 
varying magnetic fields. Section IV. A. II describes the use 
of magnetic field ramps while Sec. IV. A. 2] discusses the ap- 
plication of oscillatory fields. These schemes have been 
applied to a variety of bosonic and fermionic atom gases. 
Section lV.A.3l dcscribes the formation of coUisionally sta- 
ble molecules of fermionic atoms through three-body re- 
combination and thermalization. 

Here we restrict our discussion to ultracold gases con- 
fined in macroscopic traps. The microscopic trapping 
sites of an optical lattice, where atom pairs can be tightly 
confined, constitute a special environment for molecule 
formation. This will be reviewed separately in Sec. lVI.B] 
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magnetic field 

FIG. 31 Illustration of experimental schemes to create ul- 
tracold molecules. The solid line marks the weakly bound 
molecular state |m), which dissociates into the continuum (in- 
dicated by the dotted line) at resonance B = Bq. In (a), the 
magnetic field is ramped across the resonance, which adia- 
batically converts two interacting atoms into one molecule; in 
(b), an oscillatory magnetic field drives the transition from 
the scattering state to the molecular state; in (c), three-body 
recombination results in molecule formation. 



1. Feshbach ramps 

Ramping an external magnetic field across a Fesh- 
bach resonance is the most commonly adopted scheme 
to form Feshbach molecules. This scheme, usu- 
ally referred to as a Feshb a ch ramp, was pr o posed 
by (Ivan Abeelen and Verhaail . Il999bt iMies et all l2000t 



iTimmermans et al 



In a simplified picture, il- 
lustrated in Fig. [3lT a). the resonant coupling between 
the scattering state and the molecular state opens up a 
way to adiabatically convert interacting atom pairs into 
molecules. The atomic gas is prepared at a field B away 
from resonance where the two atoms do not have a weakly 
bound state. In Fig. ISlT a). this corresponds to i? > Bq. 
The field is then ramped to a final B < Bq to make a 
Feshbach molecule. 

( Regal et'd\ . l2003al) created ultracold molecules in a 
degenerate Fermi gas of ""^K, exploiting the resonance at 
224 G in the ac channel. In the 10-G ramp across the 
resonance with a ramp speed of 25G/ms, about 50% of 
the atoms were converted to molecules. The experimen- 
tal signatures, shown in Fig. I32r a). were a disappearance 
of atoms when the field was ramped below 224 G and a 
recovery of the atoms when the field was ramped back. 
Since the Feshbach bound state exists below 224 G, their 
observation strongly suggests formation and dissociation 
of Feshbach molecules below and above 224 G, respec- 
tively. For the 202-G resonance in the ab channel of ^"K, 



( Hodbv et aZ.I . [20051 ) reported conversion efficiencies of up 
to 80%. This resonance was also used for the formation 
of a molecular BEG; see Sec. IIV.B.II 

Feshbach ramps were also applied to degenerate 
Fermi gases of ^Li in the lowest two spin states {ab 
channel), using both the narrow resonance at 543 G 
Strecker et all 20031) a nd the broad resonance at 834 G 



Cubizolles et a/.l . l2003h . T hese two resonances in Li are 



extensively discussed in Sec. lII.BTsl Both experiments re- 
vealed a remarkable coUisional stability of the molecules, 
see Sec. IV.B.3I 

In bosonic gases, an efficient atom-molecule conversion 
by a Feshbach ramp has to overcome inelastic collision 
loss during the ramp process; see Sec. IV.B.3I In general, 
ramps applied to bosonic atom samples do not only cre- 
ate weakly bound Feshbach molecules, but they also lead 
to atom loss that cannot be recovered by a reverse field 
ramp. These atoms are lost presumably to deeply bound 
molecular states as a result of atom-molecule inelastic 
collisions. To optimize the Feshbach molecule fraction, 
which corresponds to the recoverable fraction, it is cru- 
cial to optimize the ramp speed. Furthermore, a fast 
separation of the molecules from the remaining atoms is 
essential. The latter can, for example, be achieved by 
Stern-Gerlach or optical methods. 

( Herbig et al 1 120031 ) created a pure sample of about 
3000 g-wave molecules from a BEG of ^"^^Gs atoms by a 
Feshbach ramp across the narrow resonance at 20 G. The 
presence of an inhomogeneous magnetic field during the 
ramp facilitated an immediate Stern-Gerlach separation 
of the molecules from the atomic cloud. The Feshbach 
ramp removed about 60% of the atoms from the BEG, 
but only 20% were found in the weakly bound molecu- 
lar state and thus could be recovered by a reverse field 
ramp. Expansion measurements on the Cs2 cloud showed 
temperatures below 5 nK and suggested phase-space den- 
sities c lose to or exceeding unity. In later experiments on 
"^Cs (iGhin et aLl . l200llKnoop et aLl . l2008HMark et al\ . 



l2007bl ) produced Feshbach molecules in a variety of dif 
ferent internal states and partial waves including s- to 
Z-wave {i = 8) states. 

( Diirr et aZl . l2004al ) apphed a similar scheme to a ^''Rb 
BEG and formed 7000 s-wave molecules near the 1007-G 
resonance with a 7% conversion eSiciency. A magnetic 
gradient was applied right after the Feshbach ramp to 
separate atoms and molecules. The experiments showed 
remarkable oscillations of the molecular cloud in the gra- 
dient field (Fig. which appear as a result of chang- 
ing magnetic moment near an avoided crossing between 
two molecular states. The creation of ^^Rb Feshbach 
molecules has led to fascinating applications in optical 
lattices; see Sec. IVI.B.ll 

( Xu et all [2003) demonstrated molecule formation 
with a 4% conversion efficiency in a large BEG of ^^Na 
atoms by a Feshbach ramp across the 907 G resonance . 
A large number of more than 10^ Feshbach molecules 
was observed. Expansion measurements on the Na2 
cloud showed that it was in the quantum-degenerate 
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FIG. 32 Molecule formation in a "^''K gas near the 224-G Fesh- 
bach resonance, (a) The number of atoms measured after a 
ramp across the resonance as a function of the final mag- 
netic field shows the disappearance of atoms having formed 
molecules. For comparison (b) shows the magnetic-field de- 
pendence of the s-wave scattering length in units of Bohr 
radius, as measured by radio-frequency spectroscopy. From 
(|Regal et aLl. 12003^ 1. 



regime. By diffraction off an optical standing wave, 
(lAbo-Shaeer et aZ.l . l2005D demonstrated the matter-wave 
co herence of Na 2 mole cules. 

( Hodbv et all l2005t) presented a model of the atom- 
molecule conversion efficiency that is valid for both near- 
degenerate fcrmionic and bosonic atoms. Figurc lMl shows 
a comparison of this model with data from experiments 



on "K and ^°Rb. The atom-molecule conversion effi- 
ciency follows a Landau-Zener like behavior 
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where Pmax is the maximum conversion efficiency, solely 
determined by the phase-space density of the atomic 
cloud. This expression reveals a simple dependence on 
the resonance parameters A and Obg, the atomic num- 
ber density n, and the ram p speed B; the dim ensionless 
prefactor a is discussed in (jKohler et 

Improved conversion techniques have been reported for 
atomic BECs. For the case of the na rrow 20-G g-wave 
resonance of i33Cs, iMMEeLall, [2Q03) successfully con- 
verted 30% of the atoms into Feshbach molecules. Their 
scheme relied on fast switching of the magnetic field from 
21 G right on the resonance, followed by a hold time of 
~10ms and a further switch to 18 G. This scheme was 
found to be superior to using an optimized linear ramp. 

Feshbach ramps have also been applied to create 
p-wave molecules in ultracold Fermi gases. Forma- 
tion of p-wave mol ecules in a ^Li gas was reported by 
(jZhang et al\ . \2Q04 ) based on the 185-G resonance in the 




FIG. 33 Motion of *''Rb2 molecules in a magnetic field gra- 
dient. Right after the Feshbach ramp, the molecules are sep- 
arated from the atoms because of the different magnetic mo- 
ments. While the atom cloud leaves the observation region 
(see first four images) the molecules undergo an oscillatory 
motion, which is due to a changing magnetic moment caused 
by an avoiding level crossing in the molecular states. The im- 
ages correspond to steps of 1 m s, and the field of v iew of each 
image is 0.24 x 1.7 mm. From llDiirr et a/.l. l2004al '). 



ab channel. ( Fuchs et Ml, 120081 ) measured the bind- 
ing energies of such p-wave Feshbach molecules in three 
channels (aa, ab, b b) using an o scillat ing magnetic field; 
see Sec. IIII.A.5I (jlnada et all . l2008l ) studied the colli- 



sional properti e s of th ese molecules in all three channels. 
( Gaebler et al ]. 1200 it ) created p-wave molecules in a K 
gas by fast switching of the magnetic field to the p-wave 
Feshbach resonances at 198.4 and 198.8 G (both in the 
bb channel) and studied the lifetimes of the molecules. 

Similarly, Feshbach ramps have also been applied 
to ultracold atomic mixtures to create h e teron uclear 
^OK^'^Rb (iNi et all l2008^ 



molecules, 



such 



as 



and 

^Li^^K (jVoigt et a/.l . l2008f) . An isotopic rubidium mix- 
ture was used to assoc iate ®^Rb*^Rb molecules by 
( Papp and Wiemanl . l2006l ) . A variety of other heteronu- 
clear molecule systems, are currently under investigation 
in different laboratories. Feshbach ramps have become 
a standard approach to create ultracold molecules and 
serve as a starting point to investigate the dynamics and 
the interaction properties of Feshbach molecules. 



2. Oscillatory fields 

Another powerful method to produce ultracold Fesh- 
bach molecules is based on a modulation of the mag- 
netic fi eld strength ( Hanna et aLl . [20071 : [Thompson et all 
l2005bl ). The oscillating field induces a stimulated transi- 
tion of two colliding atoms into a bound molecular state; 
see Fig. IHTTb). Heating and atom loss are reduced since 
association occurs at a bias field B away from the reso- 
nance p osition Bp. For a ^^RbBEC near the 155-G reso- 
nance, ( Thompson et aLl . l2005bt ) report on high conver- 
sion efficiencies for molecules with binding energies on the 
order of 10 kHz; the molecule formation was inferred from 
the o bservation of a resonant loss signal. ( Lange et all 
I2009D used the same method to explore weakly bound 
molecular states of Cs atoms in the aa channel in an 
energy range of up to 300 kHz. 

In a "^^K spin mixture, p-wave molecules were pro- 



40 




40 






20 




40^ 


0.0 0.5 1 


0.4 


0.6 


C 



Peak phase space density 



FIG. 34 Illustration of the dependence of the atom-molecule 
conversion efficiency on the atomic phase-space density. The 
main figure shows the results obtained in a near-degenerate 
bosonic gas of *^Rb atoms, while the inset shows correspond- 
ing results on a fermionic gas of ^''K atoms. The lines refer 
to a theoretical model that is based on the pha se-space over- 
lap o f atom pairs in the trapped gases. From (|Hodbv et aZ.I . 
l2005l l. 



duccd with an oscillating magnetic field near the 198- 
G re sonance doublet in the hh channel i Gaebler et al\ . 
I2007D . The method was also applied to produce het- 
eronucle ar s-wave molecules in mi xtures of the two Rb 
isotopes] PapD and Wiemanl . [200^ and of ^^K and *^Rb 



atoms (jWeber eiaLl . l2008f) . 



In contrast to the magnetic-field modulation method, 
radio-frequency transitions in the range of tens of MHz 
that involve a change of spin channel can also be used 
to associate two atoms to make a Feshbach molecule. 
This is the inverse of the dissociat ion process de s cribed 
in Sec. IIII A.5I In this way (jKlempt et al\ . l2008t 
IZirbel et aLl . l2008bf) achieved association of ^^Rb and 



in a dipole trap. Similar association experiments in an 
optical lattice are described in Sec. IVI.B.ll 
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FIG. 35 An ultracold ®Li gas approaches a chemical atom- 
molecule equilibrium on the molecular side of the 834-G 
Feshbach resonance. The experiment starts with a non- 
degenerate, purely atomic gas at a temperature of 2.5 pX. 
and a peak atomic phase-space density of 0.04. The mag- 
netic field is set to 690 G, where the molecular binding energy 
corresponds to 15 /xK. A^at and N^o\ denote the number of 
unbound atoms and the number of molecules, respectively. 
The total number of unbound and bound atoms 2^,1101 + A*'at 
slowly decreases because ine lastic loss is not fully suppressed. 
From (jJochim et ali l2003al ;i . 



balance of these processes one can intuitively understand 
that molecule formation is favored at low temperatures 
and high number densities, i.e. at high phase-space densi- 
ties. Indeed, for a non-degenerate gas, the atom-molecule 
equih brium follows a simple relation ( Chin and Grimml . 
l2004f ) 



»mol 



^at exp 



(62) 



3. Atom-molecule thermalization 

A particular situation for molecule formation arises in 
a spin mixture of ^Li near the 834-G Feshbach resonance. 
At the low-field side of the resonance there is a broad field 
range, where the s-wave scattering is large and positive. 
Here a weakly bound state exists with a pronounced halo 
character. The molecular state shows an extraordinary 
stability against inelastic decay, which opened the way 
to efficiently create molecular BECs by straightforward 
evapor ative cooling at a constant magnetic field near 
764 G ijjochim et ai.l . l2003bl : Iz^ierlein et aL] . l2003h . 

The formation of molecules in this region can be un- 
derst ood in terms of a chemical atom-molecule e quilib - 
rium (jChin and Grimml . l2004t iKokkelmans et ali l2004h , 
where exoergic three-body recombination events compete 
with endoergic two-body dissociation processes. From a 



where ^moi and </iat denote the molecular and atomic 
phase-space densities, respectively. The Boltzmann fac- 
tor, determined by the ratio of the molecular binding 
energy Ei, and the thermal energy /sbT, enhances the 
fraction of molecules and can partially compensate for a 
low atomic phase-space density. 

The thermal atom-mo lecule equilibrium w as experi- 
mentally investigated by ( Jochim et a /.'. '2003a') in a non- 
degenerate gas of ^Li atoms. Figure [55] illustrates how 
the initially pure atomic gas tends to an atom-molecule 
equilibrium. The observation that more than 50% of the 
atoms form molecules at a phase-space density of 0.04 
highlights the role of the Boltzmann factor. The chem- 
ical atom-molecule equili brium also played an e ssential 
role in the experiment by ( GubizoUes et aZ.I . |2003| ) . where 
a slow Feshbach ramp, which kept the sample in thermal 
equilibrium, led to a conversion efficiency of 85%. 
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B. Properties 

1. Dissociation and detection 

A general way to detect Feshbach molecules is their 
controlled dissociation through reverse magnetic-field 
sweeps, followed by imaging of the resulting cloud of 
atoms. If the image is taken immediately after the forced 
dissociation, it just reflects the spatial distribution of the 
molecules before the onset of dissociation. However, if 
the image is taken after a certain time of flight, it will 
be strongly affected by a release of kinetic energy. The 
image then contains additional information on the disso- 
ciation process. 

A reverse Feshbach ramp brings the molecule into a 
quasi-bound state above the dissociation threshold, from 
which it decays into two atoms in the continuum. The de- 
cay rate T{E)/h depends on the energy E above thresh- 
old and can be calculated from Eqs. ([T7|), and 



T{E) = 2fcabgd>A = 2kaEs, 



(63) 



( Mukaivama et all . |2004[ ) gave the energy spectrum of 
hot atoms by using this Fermi's golden rule expression for 
a linear ramp and showed that it agreed w ell with mea- 
surem ents with ■^•^Na2 Feshbach molecules. ( Goral et al\ . 
|2004[ ) verified the golden rule theory with a full quantum 
dynamics calculation of Feshbach molecule dissociation. 

The mean kinetic energy released in a reverse Fesh- 
bach ramp corresponds to the typical energy that the 
molecules can reach in the quasi-bound level before the 
dissociative decay t akes place. In experiments on ^''Rb2, 
( Diirr all . l2004bl ) studied the dependence of the energy 
release on the ramp rate and on the resonance width. 
They demonstrated kinetic energy measurements after 
reverse Feshbach ramps as a powerful indirect tool to de- 
termine the widths of weak resonances with A in the mG 
range, where direct methods are impractical. They also 
demonstrated the production of a monoenergetic spheri- 
cal wave of atoms by rapidly switching the magnetic field 
instead of ramping it. 

The dissociation properties of Feshbach molecules 
can provi de add itional spectroscopic information. 
(IVolz et all l2005l ) observed interesting dissociation 
patterns of *^Rb2, when the molecular state was brought 
high above the threshold with fast jumps of the magnetic 
field. The patterns, shown in Fig. 1361 reveal a d-wave 
shape resonance. The dissociation of ^^'^Cs2 molecules 



in Z-wave states was observed by ( Mark et al\ . l2007bt ). 
The dissociation pattern showed a strikingly different 
behavior from molecules in a g-wave state and allowed to 
clearly distinguish between these two types of molecules. 

We note that direct imaging of Feshbach molecules is 
not feasible in most situations because of the absence of 
cycling optical transitions. An exception, however, is the 
direct i maging of atoms in halo states as demonstrate d 
for ^Li (jBartenstein et aLl . [2004bl : IZwierlein et aZ.l . l2003t) . 
In this special case, the extremely weakly bound dimer 
absorbs resonant light essentially like free atoms. 
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FIG. 36 Dissociation patterns of Rb2 molecules, showing 
the interference of s- and d- partial waves. At small magnetic 
field offsets B — Bq (values given in the upper left corners), 
the s-wave pattern dominates; at large offsets, d-waves are 
strongly enhance d due to a d-wave shape resonance. From 
(|Volz fflm2005l ). 



2. Halo dimers 

Broad, entrance-channel dominated s-wave reso- 
nances feature a considerable region where Feshbach 
molecules acquire u niversal properties; see Sec. |ll] and 
(|Kohler et aLl . l2006l ). "Universality" means that details 



of the interaction become irrelevant and that all proper- 
ties of the dimer are characterized by a single parameter, 
the s-wave scattering length a or alternatively the bind- 
ing energy Eb = h^/{ma'^). The reason for this great 
simplification is the fact that the wavefunction extends 
far out of the classical interaction range of the poten- 
tial. States of this kind have been coined "quantum ha- 
los" . They have attracted considerable attention in nu- 
clear physics and, more recently, in molecular physics an d 
have been extensively reviewed in ( Jensen et all [200^ . 
An early example is the deuteron, where the neutron and 
proton are likely to be found outside of the classically al- 
lowed region. Before the advent of Feshbach molecules, 
the most extended halo system experimentally accessible 
was the helium dimer (''^He2), which is about ten times 
larger than typical diatomic molecules. An extreme ex- 
ample is given by Bosc-condensed ^Li2 Feshbach dimers 
with a size of a/2 ~ 2000 aq, which exceeds the van der 
Waals length i?vdw — 30 qq by almost two orders of mag- 
nitude. 

First experiments on halo dimers have been con- 
ducted with bosonic *^Rb (jPonlev et aZ.l . l2002D and "^Cs 
(|Mark et all l2007a[) . and fermionic and ^Li; see 
Sec. IIV.B.II For ^^Rb, the dimers are not formed from 
atoms in the lowest internal states and thus have open 
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decay channels. This leads to spontaneous dissocia- 
tion without the presenc e of other atoms or molecules. 



(|Thompson et all l2005ah aave observed an a scaline; 



of the dissociation rate, which can be understood as a 
direct consequence of universality thro ugh wavefunction 
overlap arguments ( Kohler et all [2003) . For halo dimers 



of ^'^^Cs2 created from atoms in their lowest internal 
states there are no open dissociation channels. These 
molecules cannot decay spontaneously bu t decay through 
collis i ons with other atom s or molecules ( Ferlaino et al\ . 
120081 : iKnooD et aLl . 120091 ). 

A future promising direction with Feshbach molecules 
in halo states is the experimental investigation of univer- 
sal few-body physics; see Sec. IVI.CI 



3. Collision properties 

Fast collisional loss is usually observed in trapped 
samples of ultracold molecules. T his has been seen 
in ex periment s with bosonic ^^Na ^ (iMukaivama al. 



87t 



Rb9 (ISvassen et all [200 61: IWvnar et all 12000 ) 



and (IChin et all I2005D . dZkbel et all \2(i(im 



measured large rate coefficients for fermionic ^"K^^Rb 
Feshbach molecules due to collisions with ®^Rb or '^''K. 
Both atom-dimer and dimer-dimer collisions are gener- 
ally found to cause strong inelastic loss, as demonstrated 
by the example in Fig. 1371 Vibrational relaxation is the 
dominant mechanism, which leads to large loss rate co- 
efficients of the order of 10"^" cm'^/s. Rate coefficients 
of such magnitude result in molecular lifetimes on the 
order of a few ms or less for densities characteristic of ul- 
tracold gases. Rearrangement reactions, such as trimer 
formation, may also play a significant role in limiting 
molecular lifetimes. 

Inelastic collision rates of Feshbach molecules in the 
very highest bound state, when they are not univer- 
sal "halo" states, are not significantly different from 
rates for more deeply bound sta tes, for which f a st in 



elastic loss has been observed ( Staanum et 



Zahzam et all 20061) a nd predicted ( Cvitas et all 



Quemon cr et aL . 120051 ). ( Hudson et al. . 2008 ) measured 



2006; 



2007 



large inelastic collision rate coefficients for vibrationally 
excited triplet ^^Rb^'^^Cs molecules colliding with ^^'^Cs 
or *^Rb atoms. They also used a simple model to help 
understand why such large rate constants are typical for 
atom-molecule vibrational relaxation. Assuming a prob- 
ability near unity for inelastic loss when the collision 
partners approach one another in the short-range region 
of chemical bonding, the overall collision rate coefhcient 
is then determined from the threshold scattering of the 
long range van der Waals potential. If effect, the rate con- 
stant is given by Eq. (jlOp . where the length b turns out 
to be similar in magnitude to the van der Waals length 
Rvdw Such a simple model gives the typical order of 
magnitude of 10"^" cm'^/s for the vibrational relaxation 
rate constant, nearly independent of the vibrational level, 
found for the ^^Rb^^^^Cs system. 



Halo molecules comprised of two unlike fermions bound 
in an s-wave state offer an exception to the rule of 
fast inelastic dimer-dimer and atom-dimer collisions. 
This has allowed stable molecular samples and even 
molecular Bose-Ein stein condensation; see Sec. IIV.BI 
( Petrov et al showed that a combination of two 

effects explains this stability. The first effect is a small 
wavefunction overlap of a halo dimer with more deeply 
bound dimer states, and the second one is Pauli suppres- 
sion in the few-body process. For inelastic dimer-dimer 
collisions, ( Petrov et all |200"3 ) predicted the rate coeffi- 
cient for inelastic loss to scale as (a/i?vdw)~^'^^ whereas, 
for the elastic part, they obtained a dimer-dimer scatter- 
ing length of 0.6 a. For the atom-dimer interaction, the 
predicted scaling of inelastic loss is {a/ Rvdw)~^'^^ and 
the scattering length is 1.2 a. 

An interesting case is the observation of stable ^Li2 
molecules created ne ar the closed- c hanne l dominated 
resonance at 543 G ( Strecker et al I I2003D . These are 
not halo molecules and would be expected to have a 
large collisional loss rate coefficient similar to molecules 
comprised of bosons. The collision properties of these 
molecules still await detailed investigation. 

A possible way to overcome harmful inelastic col- 
lisional loss is the appl ication of an optical lattice 
( Thalhammer et al 1 1200611 . Here a pair of atoms or a 
single molecule can be trapped in an individual lattice 
site, which offers shielding from collisions with other 
molecules or atoms. Many experiments on ultracold 
molecules are now being performed with Feshbach reso- 
nances and molecules in an optical lattice; see discussion 
in Sec. IVI.BI Another way to prevent inelastic collisions 
of ultracold molecules is to transfer them to their lowest 
energy ground state, where they do not undergo vibra- 
tional, rotational, or spin relaxation. However, reactive 
collisions may still be possible. 



4. Internal state transfer 

A Feshbach resonance can serve as an "entrance gate" 
into the rich variety of molecular states below threshold, 
allowing preservation of the ultralow temperature of the 
atomic gas that is used as a starting point. The magnetic 
association technique (Sec. IV.A.ip produces a molecule 
in a specific weakly bound molecular state, i.e. the partic- 
ular molecular state that represents the closed scattering 
channel of the resonance. This leads to the question how 
a Feshbach molecule can be transferred to other states 
with specific properties of interest or, ultimately, to the 
absolute ro-vibrational ground state. Various methods 
have been developed for a controlled internal state trans- 
fer based on magnetic field ramps, radio- frequency or mi- 
crowave radiation, or optical Raman excitations. 

When a Feshbach ramp after initially associating the 
molecules is continued over a wider magnetic field range, 
the molecule will perform a passage through many level 
crossings; see e.g. Fig. [TH The ramp speed controls 
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FIG. 37 Fast coUisional decay of ^''Na2 Feshbach molecules. 
The molec ules are trapped alone (a ) or together with atoms 
(b). From (jMukaivama ef aLL[2003) . 



whether crossings are traversed diabatically (fast ramp) 
or whether they are foUowed adiabatically (slow ramp). 
(iMark et all l2007aD demonstrated the controlled trans- 
fer of ^'^^Cs2 molecules into different states by elaborate 
magnetic field ramps. In this way, they could populate 
various states from s- up to /-waves with binding ener- 
gies of up to ~10MHz. In practice, finite ramp speeds 
limit this method to rather weak cr ossings with energy 
splittings of up to typically 200kHz. ( Lang et aLl . l2008al ) 
showed how this problem can be overcome with the help 
of radio-frequency excitation. They demonstrated the 
transfer of ^^Rb2 molecules over nine level crossings when 
the magnetic field was ramped down from the 1007-G res- 
onance to zero field in 100 ms. This produced molecules 
having a binding energy Ei, = h x 3.6 GHz with a total 
transfer efficiency of about 50%. 

More deeply bound states can be reached by two- 
photon Raman transitions, as implemented in a very 
efficient wa y by stimulated Rama n adiabatic passage 
(STIRAP) (jBergmann et all Il998f) . Fig. [38] illustrates 
STIRAP between the two highest vibrational levels in 
^^Rb2 with binding energies corresponding to 24 MHz 
and 637 MHz, a s demonstrated in a proof-of-principle 
experiment by ( Winkler et all l2007l ). The experi- 
ment also highlighted the great potential of the ap- 
proach to combine Feshbach association with stimu- 
late d Raman optical transitio ns, as originally suggested 
by ( Kokkelmans et oil . l200lt ) to create deeply bound 
molecules. 

In 2008, enormous experimental progress was made 
in applications of STIRAP to transfer both homo- and 
hetcronuclear Feshbach molecules into deeply bound 
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FIG. 38 STIRAP with Feshbach molecules, (a) After cre- 
ation of Feshbach molecules in a weakly bound state \a), a 
two-photon transition is used to transfer the molecules into 
a more deeply bound state \g). In the counterintuitive pulse 
sequence of STIRAP, laser field 2 is switched on before laser 
field 1, and the molecule is always kept in a dark superposi- 
tion state of I a) and \g) during the transfer sequence. This 
suppresses the unwanted electronic excitation of state \b) and 
leads to a very high transfer efficiency. The measurements 
shown in (b) refer to an experiment in *^Rb2, where \ a) and \g) 
correspond to the last and second-to-last bound vibrational 
level. The data show the measured population in ja) during a 
0.2-ms transfer sequence and a reverse transfer sequence after 
a 5-s hold time in state \g). The observed efficiency of 75% 
for a double passage correspond s to single-passage t ransfer 
efficiency of 87%. Adapted from ijWinkler ei a/.Ll2007l ). 



states. ( Danzl et al I [looi) explored "^Csj molecules 
and demonstr ated large b i nding e nergies corresponding 
to 31.8 THz. (jLang et all . l2008bl ) reached the rovibra- 
tional ground state in the triplet potential of *^Rb2, the 
binding energy of which corresponds to 7.0 THz. The het- 
cronuclear case was suc cessfully explored with "^^K^^Rb. 
Initial experiments by ( Ospelkaus et all . l2008l ) demon- 
strated the transfer to states with a binding energy corre- 
spondi ng to 10.5 GHz. Only shortly afterwards, the same 
group ( Ni et al\ . l2008f ) demonstrated polar molecules in 
both the triplet and the singlet rovibrational ground 
state, where the binding energies correspond to 7.2 THz 
and 125 THz, respectively. These experiments opened up 
a promising new research field related to the exciting in- 
teraction properties of ground-state molecular quantum 
gases. 

All the above methods for controlled state transfer 
rely on coherent processes. Therefore they can also be 
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applied to produce coherent superpositions of molecu- 
lar states. This can, for example, be used for pre- 
cise in terferometric measu reme nts of the molecular struc- 
ture. ijMark et aZ.l . l2007bl ) and (|Lang et aZ.l . l2008al ) inves- 
tigated moleculaxjevelcrossings in this way. Using STI- 
RAP, ( Winkler et aZ.I . [20071 ) created quantum superposi- 
tions between neighboring vibrational states and tested 
their coherence interferometrically. 



VI. RELATED TOPICS 

A. Optical Feshbach resonances 

Magnetic fields have proven to be a powerful tool to 
change the interaction strength or scattering length be- 
tween ultracold atoms. As discussed at length in this 
Review this has been made possible by the presence of a 
molecular bound state that is resonantly coupled to the 
colliding atom pair. The width of the resonance (A in 
Eq. ([1])), however, is governed by the interatomic forces 
between the two atoms. Optical Feshbach resonances 
promise control of both the resonance location and its 
width. 



1. Analogies 

Figure [39] shows a schematic diagram of an op- 
tical Feshbach reson ance. As first proposed by 
( Fedichev et ^,[19963) laser light nearly resonant with a 
transition from a colliding atom pair and a ro- vibrational 
level of an excited electronic state induces a Feshbach 
resonance and modifies the scattering length of the 
two atoms. Excited electronic states dissociate to one 
ground- and one clcctronically-excitcd atom for large in- 
teratomic separations. For many kinds of atoms the pho- 
ton needed to reach such states is in the visible or opti- 
cal domain and, hence, the term "optical Feshbach reso- 
nance" has been adopted. 

The location and strength of an optical Feshbach res- 
onance is determined by the laser frequency v and in- 
tensity /, respectively. Both can be controlled experi- 
mentally. There is, however, a crucial difference between 
magnetic and optical Feshbach resonances. For optical 
Feshbach resonances the resonant state has a finite en- 
ergy width 7 and thus lifetime h/^ due to spontaneous 
emission. Hence the scattering length becomes a complex 
number. 

Note that by changing the laser frequency and si- 
multaneously detecting the population in the excited 
electronic potentials the ro-vibrational level structure of 
these potentials can be studied. This is called ultracold 
photoassociative sp ectroscopy and has been reviewed in 
(| Jones a /],l20G6f). 

( Bohn a nd JulieEEi, Il996l . I1999D obtained expressions 
for the complex scattering length in Eqs. ([M)) - (|25|) . Reso- 
nances are characterized by a width T{E) and shift hSvc, 
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R (arb. units) 



FIG. 39 Schematic of an optical Feshbach resonance. Panel 
(a) refers to two colliding ground-state atoms, \g} that are cou- 
pled through a laser at frequency v to a. vibrational level \tpc) 
of an excited electronic state that dissociates to one excited- 
state |e) and one ground-state \g) atom. The energy hfc is 
the energy of the vibrational level relative to two colliding 
atoms at zero collision energy. The excited bound state \ipc} 
irreversibly decays with an (energy) width 7. Panel (b) shows 
the dressed state picture corresponding to panel (a). The ini- 
tial state is an atom pair with n photons and the closed chan- 
nel state |C) represents \tpc) with n — 1 photons. The atom 
pair has a relative kinetic energy E and Ec = hvc — hv. In 
the dressed state picture the ground and excited interatomic 
potentials cross at a point called the Condon point. 



where = hvc — hv — hSvc- The width is 

T{E) = 2TT\{C\d-S\E)\\ (64) 

where £ is the electric field of the laser and d is the molec- 
ular electronic transition dipole moment. Both the width 
and the shift of the resonance are proportional to /. As 
in magnetic resonances \E) is the scattering wavefunction 
at collision energy E in the entrance channel. 

Figure[n]in Section llI.A.3l shows the real and imaginary 
part of the scattering length a ~ ib for an optical Fesh- 
bach resonance. The numbers are based on an analysis of 
the strength and lifetime of an expe rimentally-o bserved 
optical Feshbach resonance in 87Rb (jTheis am 20041. 
For the intensity used in the figure the optical length a^cs 
defined by Eq. ^ is 5.47 nm and To/h = 21 MHz. That 
is, arcs ~ a-hg and Fo ~ 7. Since fcobg ^ 1, the width 
r(iJ) <C 7 so there is negligible power broadening. Un- 
like for a magnetic Feshbach resonance, the real part of 
the scattering length is now finite for any detuning with 
a peak to peak variation of 2ai.cs- Tbc length b peaks 



at zero detuning. The maximum value is 2ai-cs and the 
full- width half-maximum is 7. 

In order for an optical Feshbach resonance to be prac- 
tical it is necessary that the change in the real part of 
the scattering length a — abg is large compared to b. This 
requires that the detuning hv — hvc is large compared to 
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7- On the other hand we would also like to ensure that 
a — Obg is at least on the order of Obg at such a detuning. 
In order to satisfy these requirements simultaneously we 
need a^-cs ^ Obg or equivalently Fq ^ 7. For the pa- 
rameters used in Fig. [6] this is not true. Section IVI.A.3I 
discusses how for alkaline-earth atoms it is possible to 
satisfy Oros 3> abg. 

By adding extra laser fields at different frequencies the 
scattering length can be further manipulated. Of par- 
ticular interest is the case where a second field is nearly 
resonant with the excited bound state \ipc) and a sec- 
ond molecular bound state. If this second bound state 
is in the electronic ground state, then the situation cor- 
responds to a Raman transiti on. The analytic expres- 
sion for the scatter i ng len gth ( Bohn and Julienn^ . Il999t 
iThalhammer et all |2005[) is 



a — ih = flbtr 



1 



nE) 



2k hv - hvQ + f^2/A2 + 1(7/2) 



(65) 



where ttbg, F, and 7 arc defined as before, and hv^ = 
hvc + hSvc- The two-photon detuning A2 is zero when 
the absolute value of the frequency difference of the two 
lasers equals the absolute value of binding energy of the 
ground bound state relative to two free atoms at rest. 
It is positive when the absolute value of the frequency 
difference is larger than the absolute value of the binding 
energy. The quantity is the coupling matrix clement 
between the bound levels in the ground and excited states 
and is proportional to the square root of the intensity of 
the second laser. 



2. Observations in alkali systems 

In a magneto -optical trap f i lled w ith cold (< 1 mK) 
atomic sodi um (iFatemi et ai . 20001) confirmed the pre- 



dictions of ( Fedichev et ai . 1996a^ . They observed the 



changing scattering length by detecting the correspond- 
ing change in the scattering wavefunction. Weak detec- 
tion lasers, with frequencies that differ from those used 
for the optical Feshbach resonance, induce a molecular 
ion signal that probes the change in the scattering wave 
function. For t he m olecular level used in the experiment 
( Fatemi et were able to deduce the strength of 

the resonance as well as the light shift Svc- They found 
an Ores of around 2 nm and Fq/Zi of around 20 MHz for 
the maximum reported laser intensity of 7=100 W/cm^. 
Fo r Na the backgroun d scattering length is Obg = 2.8 nm. 
( Theis et aZ.I . |200^ tuned the scattering length by op- 



tical means in a Bose-Einstein condensate. In a Rb 
condensate they were able to change the scattering length 
over one order of magnitude from 0.5 nm to 10 nm. The 
parameters of the optical resonance are given in the cap- 
tion of Fig. [HI The scattering length was measured using 
Bragg spectroscopy, where a change in condensate mean 
field energy proportional to the scattering length is mea- 
sured by a change in the frequency that determines the 
Bragg condition. 
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FIG. 40 Optical Feshbach resonance using a two-color Raman 
transition in an *^Rb BEG. Panel (a) shows the measured 
atom number after a 100 us Raman pulse as a function of 
the two- photon detuning A2. Panel (b) shows the measured 
Bragg resonance frequency and Panel (c) shows the scattering 
length, as determined from Panels (a) and (b). The open 
and filled circles correspond to the data. As indicated by 
the error bar in Panel (b) the frequency uncertainty in the 
Bragg spectroscopy is smaller than 100 Hz. The solid lines 
have been obtained from a fit to Eq. (|65|l adapted to include 
the finite linewidth of the ground-bound state. The dashed 
line in Panel (b) shows the expected signal if there was only 
loss in atom number but no change in scattering length. The 
vertical line indicates the location of maximal loss in Panel 
(a) and helps to compare the relative positions of the three 
curves. From (|Thalhammer et aZ.l . l200 



In ( Thalhammer et al I I2OO5D the scattering length 
was modified by a two-color Raman transition. As in 
(iTheis et all \2004 the change in scattering length was 
observed in an ®^Rb condensate and detected by Bragg 
spectroscopy. In fitt i ng to Eq. ([65)) a complication arose. 
([Thalhammer et all |2005[ ) could only explain their ob- 
servations if they assumed that the target ground state 
had a finite linewidth. It turned out that, even though 
this state cannot be lost by spontaneous emission, it can 
absorb a photon from the (strong) laser that photoasso- 
ciates the scattering atom pair. This process gave rise to 
a linewidth of 2 MHz. The results are shown in Fig. HOI 
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3. Prospects in alkaline-earth systems 

So far, the scattering length has been experimentaUy 
modified by optical means in ultracold alkali-metal atom 
collisions, fn these experiments the optical length arcs 
was of the same order of magnitude as the background 
scattering length Cbg so that changes in sc attering length 
were accompanied by large atom losses. ( Ciurvlo et al\ . 
I2OO5I ) showed that optical Feshbach resonances in ultra- 
cold alkaline-earth atom collisions can have Oi-cs 3> abg. 
The presence of intercombination lines in alkaline-earth 
systems make this possible. Atomic intercombination 
lines are transitions between the ground ^Sq state and the 
excited '^Pi state. The transition is only weakly allowed 
by virtue of relativistic mixing with the ^Pi state. For 
example, for Sr isotopes 7//1 = 7.5 kHz is much smaller 
than for alkali atoms. 

When 7 is very small, it is possible to use excited 
bound levels that arc very close to the excited state disso- 
ciation threshold, while simultaneously maintaining the 
large detunings that are necessary to suppress losses. Us- 
ing levels close to threshold allow a very large value of the 
ratio ro/7, and conseque ntly Orps can be orde rs of mag- 
nitude larger than Obg. ( Ciurvlo et al\ . l2005t ) illustrate 



this by model calculations of optical lengths of ultracold 
calcium, for which ^/h = 0.7 kHz. Using a model that 
assumes a level with a binding energy on the order of 100 
MHz, they predict that a^cs could be as large as 100 nm 
at the relatively low i ntens ity of / = 1 W/cm^. 

( Zelevinskv et al\ . l2006t ) obtained photoassociation 
spectra near the intercombination line of **Sr and mea- 
sured the strength of various transitions. They found 
that the last bound state of the excited potential had an 
optical length a^^s = 24 /im at / = 1 W/cm^. The very 
large value of Ores implies that practical changes in the 
scattering length should be feasible in this species. Sim- 
ilar photoassociation spectra ha ve been observed for two 
different isotopes of Ytterbium ( Toio "all . l2006[ ). which 
has electronic structure like that of the alkaline earth 
atoms. Optical control of both bosonic and fcrmionic 
isotopic specie s may become possible with alkaline earth 
or Yb atoms. ( Enomoto et aZ.l . l2008( ) have demonstrated 
optically induced changes in scattering length for ^^^Yb 
and i^^Yb. 



B. Feshbach resonances in optical lattices 

Ultracold atoms in optical lattices arc of great interest 
because of the exciting prospects to simulate a variety of 
condensed matter phenomena , to realize la, rge scale quan- 
tum information processing (|Blochl . [2001 . and to form 
ultracold molecules in individual lattice sites to avoid 
detrimental collision instability. In all these research di- 
rections, Feshbach resonances will provide excellent tools 
to control the interaction of the constituent atoms and to 
explore the transition between different quantum regime 
and quantum phases. In the following sections we review 



atom-atom scattering in optical lattices and describe the 
role of Feshbach resonances therein. 

Optical lattices are realized by standing-wave laser 
fields, which result in spatially-periodic potentials for 
the atoms. The atoms are confined in the individual po- 
tential minima or sites of the lattice potential. One-, 
two-, and three-dimensional lattices can be created in 
this way. In experiments with three-dimensional config- 
urations the optical lattice can be filled with only one or 
two atoms per site. 

Section fVl . B . 1 1 discusses that when two atoms are held 
in a single lattice site a Feshbach resonance can be used 
to very efficiently produce stable molecules. One advan- 
tage of lattice confinement is that such molecules are pro- 
tected from harmful coUisional losses with a third body. 
Section rVl.B. 21 discusses the possibility that confinement 
in one spatial direction can induce resonant behavior in 
scattering along the remaining directions. Finally, Sec- 
tion IVI.B.3I describes uses of Feshbach resonances in op- 
tical lattices where tunneling between lattices sites is im- 
portant. 



1. Atom pairs and molecules 

When an atom pair is trapped in a single site of a three- 
dimensional optical lattice the motion is fully quantized. 
For deep optical lattices the confining potential is har- 
monic, and the center-of-mass and relative motion of the 
atom pair separate. In other words the six-dimensional 
wavefunction of the two atoms becomes a product of a 
center-of-mass and relative wavefunction. The center-of- 
mass motion is harmonic and is solved trivially. The rel- 
ative motion is determined by a potential that is the sum 
of the atom-atom interaction potentials and a harmonic 
potential. 

The atom-atom interactions between alkali-metal 
atoms are independent of the relative orientation of the 
atoms when very weak spin-dependent interactions, T4s, 
are ignored. Consequently, for a spherically symmet- 
ric harmonic trapping potential the three-dimensional 
relative motion can be further simplified. The angu- 
lar motion can be solved analytically and only a radial 
Schrodinger equation for the atom-atom interaction po- 
tential plus iiuP'B? jl needs to be solved. Here, /i is the 
reduced mass and uj is the oscillation frequency in the 
trap. 

Figure WS\ shows cigenenergies for two Na atoms with 
zero relative orbital angular momentum = 0) in a 
spherically-sym metric harmonic tra p as a function of 
magnetic field ( Tiesinga tt al\ . I2OOOI ) . The atoms are in 
their lowest hyperfine state and the energies are obtained 
from coupled-channels calculations. For these atomic 
states there is a Feshbach resonance near 910 G (91 mT) . 
The zero of the vertical axis corresponds to zero rela- 
tive kinetic energy in the absence of a trapping potential. 
Hence, positive energies correspond to atoms in the trap 
and negative energies correspond to molecules bound in 
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FIG. 41 The energy in the relative motion of two trapped 
interacting Na atoms in their lowest hyperfine state \a), as a 
function of magnetic field. The trapping frequency f = uj/2n 
is 1 MHz. The full lines correspond to energies obtained 
from exact numerical calculations. The dotted lines corre- 
spond to eigenenergies for trapped Na atoms interacting via 
a regularized delta-function potential with a magnetic field 
dependent scattering length given by the inset. This inset 
shows the exact scattering length for two freely scattering \a) 
states near a Feshbach resonance. In the theoretical model 
-Bo = 90.985 mT, whereas experimentally Bo = 90.7 mT; see 
Section Fill. B. 21 The long-dashed lines correspond to energies 
of the l=Q, '^~9; ^? ^^'^ ^ harmonic oscillator states. From 
(jTiesinga 



the atom-atom interaction potential. Avoided crossings 
between the closed channel Feshbach state and the trap 
levels are clearly visible. 

An analytic solution for two intera cting atoms 



trapp ed in a harmonic trap was found by (jBusch et al\ . 
Il998l ). The atom-atom interaction was modeled 
by a so-called regularized delta-function potential, 
ATTn'^/{2n)a6{R){d/dR)R, where a is the scattering 
length. This approximation of the actual interaction po- 
tential is valid in the Wigner threshold regime. In partic- 
ular, they showed that the eigenenergies for the relative 
motion are given by 



a _ r(-^/2 + 1/4) 
cr ^ 2r(-£;/2 + 3/4) 



(66) 



where F is the Gamma function, the energy E is in units 
of hoj , and a = ^JhJJJujj) is the harmonic oscillator length 
for the relative motion. 

Near a magnetic Feshbach resonance the scattering 
length is a rapidly-changing function of i3, as for ex- 
ample shown in the inset of Fig. 21] The dotted lines 
in Fig. [JT] are the eigenenergies found from combining 
Eq. ([M)) and the a{B) for the resonance. The exact 



and model calculations do not agree where the scattering 
length is large. A s shown by ( Blume and Green3 . l2002t 
iBolda et aLl . |2003^ this is due to the breakdown of the 
Wigner threshold regime at the finite zero-point energy o f 
the atoms near the resonance. ( Blume and Gree"ni . 120021 ) 
introduced an energy-depende nt scattering l ength , based 
on the effective range theory. (|Bolda et all l2002h found 
that an energy-dependent regularized delta-function po- 
tential reproduces well the exact results in Fig. 141] 

(|Dickerscheid et all l2005t iGubbels et di |2006[ ) de- 
veloped an analyti c approach extending the theory of 
( Busch et gLl . llQQSh to the states of two trapped atoms in 
a single lattice site that interact strongly through a Fesh- 
bach resonance. They applied their two-body theory to 
the broad ^Li resonance near 834 G. They also showed 
how to incorporate this theory into a many-body Hub- 
bard model that treats the tunneling of atoms between 
lattice sites. 

( Mies e^"aZI . l2000l ) have theoretically shown that by 



varying the magnetic field in time two atoms in a single 
lattice site can be converted into a molecule with near 
100% efficiency. The idea is to prepare the atoms in the 
lowest trap level at a magnetic field where the Feshbach 
state has a higher energy. In Fig. [41] this corresponds 
to the nominally n = state at for example -6=912 G. 
The magnetic field is then varied. If the ramp is suf- 
ficiently slow the atom pair will be adiab atically con- 
verted into a molecule. ( Mies et "oil. I2000t) also showed 
that this process can be described by a Landau-Zcner 
cur ve crossing model. A more recent derivation is given 
in (| Julienne et"aZI, |2004[ ). 

( Widera etaR - hm^) used a magnetic Feshbach reso- 
nance to entangle two ^^Rb atoms in a site of an optical 
lattice. In a Ramsey-type interferometer a sequence of 
microwave pulses manipulate and control the superposi- 
tion between the |0) = \F = l^mp = 1) and |1) = |F = 
2,mp ^ —1) hyperfine states of each ^'^Rb atom. Initially 
the atoms are in state |0). The population in the two hy- 
perfine states after the pulses will depend on the atom- 
atom interaction, which entangles the two atoms. By 
controling the scattering length by t he magnetic field and 
the hold time between the pulses, ( Widera et al\ . 12004] ) 
were able to create maximally-entangled Bell states. 

( Stoferle et al\ . 12006] ) spectroscopically mapped the 
avoided crossing between the Feshbach and the low- 
est harmonic oscillator state as a function a magnetic 
field. They prepared fermionic ^"K atoms in a three- 
dimensional optical lattice configured such that the bot- 
tom of the lattice sites are spherically symmetric. The 
two atoms in each lattice site are in different hyperfine 
states. They confirm that the model of trapped atoms 
interacting via an energy-dependent delta- function po- 
tential agrees with the experimental observations. 

( Thalhammer et all 12006] ) showed in an experiment 
with ^^Rb in its lowest hyperfine state that the Landau- 
Zener model for a time-dependent sweep of the magnetic 
field through t h e Feshbach resonance, as developed by 
(|julienne et all \2004 iMies et all l2000f) . is vahd. The 
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FIG. 42 Energy of heteronuclear ''"K - ^'^Rb molecules versus 
B in an optical lattice for two different lattice depths C/iat in 
units of the *^Rb recoil energy Er — hP'k'^ /2m-iih (Er/h « 
3.7 kHz). The center of the Feshbach resonance is located 
at Bo = 546.8(1) G. The zero of energy corresponds to two 
atoms in the lattice cell with B far from Bo . Molecules which 
are stable in free space are observed for B < Bo where a > 
0. Confinement-induced molecules are observed for B > Bo 
where a < and no threshold bound state exists in free space. 
In addition , repulsively interacting pairs are observed for B < 
Bo. From (|Ospelkaus et a/.ll2006al '). 



Feshbach resonance near 1007 G was used. In the exper- 
iment the ramp speed dB /dt was varied over four orders 
of magnitude and for the slowest ramp speed of 2 x 10'^ 
G/s a 95% conversion efficiency was observed. The exper- 
iment demonstrated a dramatic increase in the lifetime of 
the trapped molecules, where the lattice protected them 
from harmful collisions. Molecular lifetimes up to 700 ms 
were observed. 



( Ospelkaus et all [2006al ) were able to make heteronu- 
clear molecules by associating atoms of two different 
species, ^"K and '^^Rb, trapped on the same lattice site. 
They used a rf association technique both to form the 
molecule and to measure its binding energy. Figure [321 
shows the energy of the near-threshold states of the atom 
pair as a function of B and illustrat es an avoided crossin; 
simil ar to that shown in Fig. [4TJ ( Deuretzbacher et al. 
I2008D developed a theoretical model to account for an 
harmonic corrections, which couple center-of-mass and 
relative motion of the atoms in the trap. 



2. Reduced dimensional scattering 

Optical lattices that confine atoms in only one or two 
directions in combination with magnetic Feshbach reso- 
nances lead to controll able quasi-2D o r quas i-lD scat- 
tering, respectively. ( Yurovskv et al\ . l2008f ) have re- 
cently reviewed such reduced dimensional scattering. By 
integrating out the confined spatial direction, effective 
one- and two-dimensional atom-atom potentials can be 
derived. Their strength is related to the magnetic-field 
de pendent scattering length for free scatt ering. 

shaniil . Il998t ) derived the ef- 



( Bergeman et 



m g leng th : 
H. 120031 lOl 



fective atom-ato m potential for qua si one-dimensional 
scattering. As in ( Busch et al\ . [ToOSl ) the starting point 
is a regularized three-dimensional delta-function poten- 
tial for the atom-atom interaction potential. The trap- 
ping potential along the two confined dimensions is the 
same and harmonic with frequency lo^. They find that 
for an atom pair in the lowest harmonic oscillator state 
of the confined directions the atoms interact via a one- 
dimensional delta-function potential giY)5{z)^ with cou- 
pling constant 



fi^ a 



9lB 



/i cr^ 1 — Ca/ai_ ' 



where C = 1.4602 . . . and (t± = ^fij {^u!±). The cou- 
pling constant is singular when a — a\_IC and ap- 
proaches the negative finite v alue g^o ~ —2h^/{^C(7i^ 
for a ±oo. ( Olshanii Il998^ has called the singularity 
a confinement-induced resonance. In practice, the res- 
onance condition can be fulfilled by changing a with a 
magnetic Feshbach resonance. For fermionic atoms in 
quasi one-dimensional confinemen t the effective atom- 
atom potential has been derived by ( Granger and Blumd . 


(jPetrov et all l2000l : iPetrov and Shlvapnikovl . l200ll ) 
derived a similar coupling constant for one-dimensional 
confinement or two-dimensional scattering. In this 
case the resonance location not only depends on a but 
also logarithmically on the relative wavenumber be- 
tween the a toms along the two free spatial directions. 
(jNaidon et al.. . .2007i) describe how these reduced dimen- 
sional treatments can be extended to much tighter con- 
finements than previously thought and made more ac- 
cura te by using the ene r gy-de pendent scattering length 
of l|Blume and Greend . I2OO2I ) with the effective range 
ex pansion for the scat tering phase. 

( Moritz et al\ . [20051 ) presented experimental evidence 
for a confinement-induced bound state in a one- 
dimensional system. They confirmed the existence of the 
bound state of the one-dimensional Hamiltonian _ffiD = 
—T? I I dz^ -I- g\\y8{z) by changing both a and u^. 
The experiment was performed by employing an ar- 
ray of ID tubes, each containing about 300 atoms, 
equally divided between the |/ = 9/2, m = —9/2) and 
1/ = 9/2, TO ~ —7/2) hyperfine state, were held in a two- 
dimensional harmonic trap with frequency w^/(27r) = 69 
kHz. The scattering length for the collision between these 
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FIG. 43 Measured molecular bound state binding energy ver- 
sus magnetic field B near the 202. 1-G resonance between the 
two lowest spin states of *°K. The upper curve labeled 3D is 
for a free space gas, whereas the lower curve labeled ID is 
for a quasi-lD gas with tight confinement in two dimensions. 
The points show measured bin ding energies and th e lines show 
theoretical predictions. From (|Moritz et aZ.I . [20051 ) . 



two hypcrfinc state was varied using the magnetic Fesh- 
bach resonance at Bq = 202 G. Figure B51 shows that the 
measured energy of the tightly confined atom pair varies 
as predicted by theory. Confinement-induced molecules 
exist in reduced dimensio n for B > Bp, where they do 
not exist in free space. (jPickerscheid and Stooi l2005l ) 
develop an analytical nonperturbative two-channel the- 
ory of the binding energy that is in excellent agreement 
with the data. 

In ijGiinter et a/.l . l2005t) fermionic K atoms were pre- 
pared in a single hyperfine state and held in either a 
one- or two-dimensional optical lattice. By virtue of 
Fermi statistics the atoms can only collide via odd par- 
tial waves, which for ultracold collision energies have 
very small cross sections. One might expect the atomic 
gas to be a nonint eracting Fermi gas. Nevertheless, 
( Giinter et ^,[200i) could observe a Feshbach resonance 
in the p-wave collision by the losses it induced. The losses 
were even sensitive to the orientation of the magnetic 
field relative to the principal axis of the trap. Moreover, 
they showed that the magnetic-field location of the reso- 
nance is modified by the confinement. 

( Nvgaard et all . l2008al ) investigate the effect of tun- 
ing a narrow Feshbach resonance across the Bloch band 
of a one-dimensional optical lattice for the case when 
the resonance width is small compared to the width of 
the band, such as the 414 G ^^Rb resonance studied 
by ( Svassen et all [20071 ) . They investigate the changes 
in scattering and bound states due to the band struc- 
ture in the periodic structure and characterize the time- 
depend ent dynamics o f sweepi ng the resonance across the 
band. ( Nvgaard et aLl . r2008bt ) extend this work, develop 
the concept of a generalized scattering length at the band 
edges, and show the existence of a "universal" bound 



state near the top and bottom band edges at the field 
strength where the resonance emerges from the band. 



3. Scattering in shallow lattices 



The previous two Sections I VI. B. II and I VI . B ."2] discussed 
deep optical lattices where the tunneling between lat- 
tices sites could safely be neglected. For weaker optical 
lattices, atoms tunnel from site to site and then inter- 
act with all other atoms. This leads to many-body sys- 
tems that can be described by either a mean-field Gross- 
Pitaevs kii equ a tion o r a Bose- or Fermi-Hubbard Hamil- 
tonian (iBlochl . I2005D . The presence of Feshbach reso- 
nances has added and continues to add new twists to 
these kinds of Hamiltonians. 

We will discuss the more simple situa tion where only 
two atoms scatter in an optical lattice. ( Fedichev et all 



I2004D studied the case of two atoms scattering in a weak 
three-dimensional optical lattice of cubic symmetry and 
interacting with the regularized delta-function potential. 
They predict the presence of a geometrical resonance in 
the 3D lattice based on a derivation of an effective atom- 
atom interaction between the atoms, which by virtue of 
the periodic potential have an effective mass m* that is 
much larger than their atomic mass m. Figure shows 
the results of their calculation. The resonance occurs at a 
scattering length a = I* = ~{Tr/{2ln2)){m/m*){a/d)'^a, 
where a is the atomic harmonic oscillator length for 
motion in a single lattice site and d is the latt i ce pe- 
riod. Note that, in practice, I* ^ a. ( Orso e^"all . |2005[ ) 
presented a similar analysis for a one-dimensional opti- 
cal lattice. ( Grupp et all l200l1 ) studied the effect of a 
very-narrow Feshbach resonance in scattering in an one- 
dimensional lattice. 



C. Efimov states and universal few-body physics 

Feshbach resonances provide experimental access to 
systems with very large values of the scattering length. 
Such systems are governed by "universal" physics, i.e. 
their low-energy ob servables are independent of details 
of the interaction ( Braaten and Hammed . [2006[) . Uni- 
versality appears as a consequence of the quantum-halo 
character of the wave function carrying its dominant part 
far out of the classically allowed region. In this case, de- 
tails of the interaction potential become irrelevant and 
the system can be described by a few global parame- 
ters. Halo dimers (Sec. IV.B.2P are the most simple ex- 
ample. For addressing universal physics with ultracold 
gases, Feshbach resonances that are strongly entrance- 
channel dominated (sres 2> 1) are of particular interest, 
as they allow a description in terms of a single-channel 
model with a large range of universal behavior; see dis- 
cussion in Sections III. B. 21 and III.C.5] 
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FIG. 44 The dependence of the effective scattering length 
acff on the absolute value of the microscopic scattering length 
\a\. The dashed line corresponds to a repulsive interaction 
potential (a > 0). The solid line corresponds to an attractive 
potent ial (a < 0) wi t h a geometric resonance at \a\ = h 
From (|Fedichev et all 12004 



1. Efimov's scenario 

Efim ov quantum states i n a system of three identical 
bosons ( Efimovl . Il970l Il97ll ) are a paradigm for universal 
few-body physics. These states have attracted consider- 
able interest, fueled by their bizarre and counter-intuitive 
properties and by the fact that they had been elusive 
to experimenta l ists f or more than 35 years. In 2006, 
( Kraemer et aLl . l2006l ) reported on experimental evidence 
for Efimov states in an ultracold gas of cesium atoms. By 
Feshbach tuning they could identify a pronounced three- 
body resonance, which occurs as a fingerprint of an Efi- 
mov state a t the three-body sc attering threshold. Two 
years later, ( Knoop et aLl . l2009l ) presented additional ev- 
idence for Efimov-like trimer states, reporting on the ob- 
servation of a decay resonance in atom-dimer scattering. 

Efimov's scenario is illustrated in Fig. [45l showing the 
energy spectrum of the three-body system as a func- 
tion of the inverse scattering length 1/a. For a < 0, 
the natural zero of energy is the three-body dissociation 
threshold for three atoms at rest. States below are trimer 
states and states above are continuum states of three free 
atoms. For a > 0, the dissociation threshold is given by 
—Eh ~ —h^/{nio?) where iJf, is the universal binding en- 
ergy of the weakly bound halo dimer; at this threshold a 
trimer dissociates into a dimer and an atom. All states 
below threshold are necessarily three-body bound states. 
Efimov predicted that in the limit a — > ±oo there would 
be an infinite sequence of weakly bound trimer states 
with a universal scaling behavior. Each successive Efi- 
mov state is larger in size by a universal scaling factor 




Inverse scattering lengtli 1/a 

FIG. 45 Efimov's scenario; Appearance of an infinite series 
of weakly bound Efimov trimer states (ET) for resonant two- 
body interaction. The binding energy is plotted as a function 
of the inverse two-body scattering length 1/a. The shaded 
region indicates the scattering continuum for three atoms 
(a < 0) and for an atom and a dimer (a > 0). The ar- 
row marks the intersection of the first Efimov trimer with 
the three-atom threshold. To illustrate the series of Efimov 
states, the universal scaling factor is artificially reduced from 
22.7 to 2. For comparison, the dashed line indicates a more 
tightly bound non-Efimov trimer (T) which does not cross the 
scattering continuum. Adap ted by permission fro m Macmil- 
lan Publishers Ltd: Nature (| Kraemer copyright 
2006. 



energy by a factor of (22.7)^ « 515. 

Efimov states exist on both sides of a resonance, and 
Fig. |45] shows the adiabatic connection between both 
sides. For a > 0, an Efimov state near the atom- 
dimer dissociation threshold can be regarded as a weakly 
bound state of an atom and a dimer with a size set 
not by a but by the even larger ato m-dimer scattering 
length (jBraaten and Hammer . I2006D. For g < 0, Efi - 
mov states are "Borromean" states (Jensen et ^,|200l, 
which means that a weakly bound three-body state ex- 
ists in the absence of a weakly bound two-body state. 
This property that three quantum objects stay together 
without pairwise binding is part of the bizarre nature of 
Efimov states. 



Resonant scattering phe nomena arise a s a natural con- 
sequence of this scenario ( Efimovl . Il979l ). and they are 
closely related to the basic idea of a Feshbach resonance. 
When an Efimov state intersects with the continuum 
threshold for a < 0, three free atoms resonantly couple to 
a trimer. This results in a "triatomic Efimov resonance" . 
When an Efimov state intersects with the atom-dimer 
threshold for a > 0, the result is an "atom-dimer Efimov 



22.7 (sq = 1.00624) and has a weaker binding resonance" (INielsen et al liool). 
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FIG. 46 Observation of an Efimov resonance in three-body 
decay of an ultracold gas of cesium atoms. The data 
are presented in terms of a reco mbination length ps = 
[2m/{V3h) LaY^'^ ijEsrv et adll999t l. The general a*-scaling 
of Z/3 corresponds to a linear behavior in p3{a) (straight lines). 
The filled circles represent measurements taken at tempera- 
tures around 10 nK, whereas the filled triangles and open dia- 
monds refer to measurements in the range of 200 nK - 250 nK. 
The solid line is a fi t to the low-temperature dat a based on 
effective-field theory (|Braaten and Hamniei1 . l200d ). The inset 
shows an expanded view of the region of positive scattering 
lengths up to 600 oq. Repri nted by permiss i on fro m Macmil- 
lan Publishers Ltd: Nature (|Kraemer et aU |2006| ). copyright 
2006. 



2. Observations in ultracold cesium 

In an ultracold atomic gas with resonant interactions, 
Efimo v physics manifests itself in three-body decay prop- 
erties (iBedaaue et al . 2000l: Braa ten an d Hammeil . l200ll 
120061: lEsrv et all . Il999l : iNielsen a nd Mace k|, 119991) . The 
three-body loss coefficient L3 fScc. lIILA.2[) can be conve- 
niently expressed in the form L3 ~ 'iC{a)'ha'^ jm^ which 
separates an overall a^-scaling from an additional de- 
pendence C(a). Efimov physics is reflected in a loga- 
rithmically periodic behavior C(22.7a) = C(a), corre- 
sponding to the scaling of the infinite series of weakly 
bound trimer states. A tria tomic Efimov resonance leads 
to giant recombin ation loss ( Braaten and Hammeil . I2OOII : 
lEsrv et all . ll999t ). as the resonant coupling of three atoms 
to an Efimov state opens up fast decay channels into 
de eply bound dim e r stat es plus a free atom. 

(|Kraemer et all I2006D observed a triatomic Efimov 
resonance in an ultracold, thermal gas of Cs atoms. They 
made use of the strong variation in the low-field re- 
gion (Fig. W2\ . This tunability results from a strongly 
entrance-channel dominated resonance at — 12G with 
Sics = 566 (Table HV|) . which provides a broad range 



of universal behavior. By applying magneti c fields be- 
tween G and 150 G, ( Kraemer et al\ . [20061 ) varied the 
s-wave scattering length a between —2500 ag and 1600 oq, 
large enough to study the universal regime, which re- 



quires \a\ i?vdw — 100 oq. The occurrence of one tri- 
atomic Efimov resonance could be expected in the ac- 
cessible negative-a region. The position, however, could 
not be predicted from knowledge of the scattering length 
alone as, for a three-body process, a second parame- 
ter is required to characteri ze the universal properties 
( Braaten and Hammed . l2006t ). 

Figure HSl shows the results of ( Kraemer et a/.l . [2006( ). 
The three-body loss resonance was found at a magnetic 
field of 7.5 G; corresponding to a scattering length of 
— 850ao. The behavior of loss at temperatures around 
10 nK close l y rese mbles the theoretical predictions of 
(|Esrv et ad il999f) . who numerically solved the three- 
body Schrodinger equation for a generic two-body model 
potential. The observed behavior is also well fit with 
a universal analytic expression obtained in the frame- 
work of effective-field theory ( Braaten and Hammed . 
|2006[ ). Experimental data taken at higher temperatures 
demonstrated the un itarity limitation of three-body loss 
( DTncao et Ml, 120041 ) and showed how the Efimov res- 
onance evolved into a triatomic continuum resonance 
(|Bringas et al\ . \2004 . 

For positive scattering lengths, theory predicts a vari- 
ation of C (a) between very small va l ues and a maxi- 
mum of about 70 (.Bcdaque et all I2OOOI : lEsrv et all Il999t 
iNielsen and Macekl . Il999[ ). The resuhs in Fig. [35] are 
consistent with the upper loss limit, represented by the 
straight line for a > 0. For a below 600 oq, the exper- 
imentally determined recombination length significantly 
drops below this limit, as seen in the inset of t he figure. 
Further measurements of ( Kraemer et ^,[200i) revealed 
the existence of a loss minimum at B = 21 G, where 
a = -|-210ao. It is interestin g to note that earlier ex- 
periments by the same group ( Kraemer et aLl . [20041 ) had 
identified 21 G as an optimum magnetic field for evap- 
orative cooling of cesium and attainment of BEG; see 
Sec. IIV.A.II The nature of the minimum may be inter- 
preted in the framework of universal physics, following 
theoretical predictions of an interfere nce effect between 
two different recombinatio n pathways ( Esrv et al\ . Il99^ 
[Nielsen and Macel^ . [1999[ ). However, as the minimum 
occurs at a scattering length which is only a factor of 
two larger than i?vdw — 100 ao (Table [l| the application 
of univ ersal theory to describe this feature is question- 
able. ( Massignan and Stooj . [2008t) presented an alter- 
native theoretical approach, which reproduced both this 
minimum and the maximum observed for negative a on 
the basis of the two-body physics of the particular Fesh- 
bach resonance. 

In a pure sample of trapped atoms, as discussed 
so far, three-body recombination is the only probe for 
Efimov physics. Mixtures of atoms and dimers can 
provide complementary information on Efimov states 
through resonances in inelastic atom-dimer co llisions 
( Braaten and Hammed. [2007t [Nielsen et aZ.I . [2002[). In a 
recent experiment, ( Knoop et al\ . 2009t ) prepared an op- 
tically trapped mixture of Gs atoms and Gs2 halo dimers. 
Their measurements revealed an atom-dimer scattering 
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resonance, which is centered at a large value of the two- 
body scattering length, a ~ -1-390 at a magnetic field 
of 25 G. This observation provides strong evidence for a 
trimer state approaching the atom-dimer threshold. The 
situation is close to the atom-dimer resonance in Efimov's 
scenario, but it probably remains a semantic question 
whether, at a ~ 4_Rvdw, the underlying trimer state may 
be called an Efimov state. 

In the Cs experiments described in this Section, the 
Efimov state that causes the observed triatomic reso- 
nance at 7.5 G does not connect to the state that causes 
the atom-dimer resonance at 25 G when the magnetic 
field is varied. This is because these two cases are 
separated by a zero crossing in the scattering length 
(Fig. [22)) and not by the pole as in Efimov's sce- 
nario in Fig. US] A universal relation between thes e 
regimes m a y nev ertheless exist ( Kraemer et al\ . |2006[ ). 
( Lee et gTI . 2007t ) provided a further interpretation of 
these observations in terms of the underlying CS3 states 
and point out the analogies to trimer states of helium 
(ISchoUkopf and Toenniesl . I1994D . 



3. Prospects in few-body physics 

Ultracold gases with resonantly tuned interactions of- 
fer many opportunities to study universal, Efimov-related 
few-body physics. Cesium alone has much more to offer 
than the experiments could explore so far. Moreover, 
several other systems with broad Feshbach resonances 
promise new insight into this field. 

In c esium, a predicte d broad Feshbach resonance near 
800 G (iLee et all [2OO7I) in the aa channel offers similar 
properties as the low-field region explored in previous ex- 
periments, but overcomes the disadvantage that only the 
tail of the resonance is accessible at low fields. The broad 
155-G resonance in ^^Rb fSec. IIII.B.4P might be another 
interesting candidate, but experiments may suffer from 
strong two-body decay which is absent for the discussed 
cesium resonances. A further interesting candidate is the 
402-G resonance in (Sec. IIII.B.3|) : (jZaccanti et all 
I2OO8I ) have studied three-body decay near this resonance 
and found features strongly indicative of Efimov physics. 

Many more opportunities for studying Efimov-related 
physics in ultracold gases with resonant interactions are 
offered by mixtures of different spin states or differ- 
ent species. In ^Li, all three combinations of the low- 
est three spin states (chann els ah, ac, and be ) have 
broad Feshbach resonances (IBartenstein et al I I2005D 
that overlap in a magnetic-field range between 650 G 
and 850 G. For such a three-com ponent fermionic spin 
mixture I Luu and SchweiiS . '2007*1 predicted a novel Bor- 
romean three-b o dy st ate. (Ottcnstc in_ej al . . 2008j) and 
( Huckans et all . l2009l ) experimentally investigated the 
stability of a ®Li three-component spin mixture and 
found evidence for a three-body resonance at 130 G. 
Its interpretation in terms of coupling to a three- 
body bound state is supported by several theoretical 



studies (iBraaten et all l2008l : iNaidon and Uedal . l2008t 
ISchmidt et al. , 20081 ). A particularly interesting situa- 
tion arises in mixtures of atoms with different masses. 
With increasing mass ratio the Efimov factor substan- 
tially decreases from its value of 22.7 at equal masses to 
values as l ow as 4.9 for the mass ra tio of ^"^^Cs combined 
with ^Li. ( D'Incao and Esrvl . l2006t ) pointed out that this 
will substantially enhance the observability of the Efimov 
effect in terms of the logarithmically periodic variation of 
the three-body loss coefficient with increasing two-body 
scattering length. 

Four-body processes at large values of the s-wave scat- 
tering length a represent a logical next step in under- 
sta nding universal few-body physics. Theoretical stud- 
ies (iHammer and Platteil . l2007t Ivon Stecher et aZ.l . l2009t 
IWang and Esrvl . |2009[ ) predict the existence of univer- 
sal four-body states, and consider the process of atomic 
four-body recombination. A first experimenta l step 
into this field was made by (jFerlaino et al I I2008D . who 
studied collisions of Cs2 halo dimers at large positive 
a. They observed a loss minimum in the same re- 
gion where atom-dimer scattering shows a maximum 
(a ~ 500 qq at 30 G) , which may be related to a uni- 
versal connection between four- and three-body phy sics 
(|Hammer and Platted . l2007t Ivon Stecher et aZ.I . l2009t) . 

Optical lattices have proven a powerful tool for 
the manipulation of ultracold Feshbach molecules, see 
Sec. IVI.B.ll and may also open up new possibilities for 
the creation of Efimov trimers and, more generally, for 
the controlled production of few-body quantum states 
(iLuu and Schwenkl . l2007t IStoU and KohleJ . I2005D . 



D. Molecular resonances and cold chemistry 

While we have reviewed the formation of cold 
molecules from ultracold atoms, parallel progress has 
been made by other techniques for preparing samples of 
cold molecules that extend the range far beyond alkali- 
metal species. These advances have been made possi- 
ble by Stark decele ration of m olecules such as ND3 , OH, 
and formaldehyde (jvan de Me cra kker et at. 2006) or by 



buffer gas cooling with liquid helium (jPeCarvalho et all 
I1999II . In contrast to the association of cold atoms to 
make Feshbach molecules, which have a high level of 
vibrational excitation, these methods can produce cold 
mol ecules in the rotational and vibrational ground state . 
See (iDovle et aLl . 120041; [Hutson and Soldanl . l2007tlKremsl . 
l2005r for overviews of the issues involved in trapping, 
cooling, and colliding such molecules. 

Resonances will play a very prominent role in atom- 
molecule and molecule-molecule collisions. The com- 
plexity of these systems will increase the number of 
closed channels and lead to numerous resonances with di- 
verse properties. Their presence will make it possible to 
change molecular scattering properties as well as to cre- 
ate more complex molecules. Both static magnetic and 
electric fields can be used to tune the molecule-molecule 
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resonance and provide control over collisions, as many 
molecules not only have a magnetic moment but an elec- 
tric dipole moment as well. Here wc will briefly review 
some of this work. 

( Forrev et al\ . Il998h pointed out that Fcshbach res- 
onances occur in ultracold atom-diatom scattering, 
giving an example from collisions of H2 with He. 
The effect of resonance s tates in chemical reactions 
has been studied by (IBalakrishnan and Dalgarno . 
20011) for F-fH2^FH-fH and (|Weck and Balakrishnan . 
2005,) for Li-l-HF^H-l-LiF reactions. Recent cou 



1.0 



pled channel s models of Rb-|-OH (iLara et al. ^ 20071 ) 



and He-|-NH I Gonzalez-Martinez and Hutsonl . 2007 ) col- 



lisions have been developed to give a realistic assessment 
of atom-molecule scattering. The latter study demon- 
strated the effect of magnetic tuning of a decaying res- 
onance across a threshold, using the resonance length 
formalism described in Sec. III. A. 31 

In 2002 (jBohn et all l2002f ) realized that, unhke for 
atomic systems where Feshbach resonances originate 
from the hyperfine structure of the atoms, for molecules 
the resonances can also be due to rotational states. For 
many molecules the rotational spacing for low-lying ro- 
tational levels is of the same order of magnitude as hy- 
perfine interactions in atoms. Based on the rotational 



splitt ings and a potential energy surface (jBohn et al\ . 
1200211 estimated the mean spacing and widths of the reso- 
nances and found for collisions between oxygen molecules 
as many as 30 resonances for collision energies below 
E/ kB = 1 K. 

(IChin et all |2005^ observed magnetic Feshbach-like 
resonances between two weakly bound ^^'^Cs2 molecules 
that temporarily form a tetramer during a collision. 
Their data are shown in Fig. [47l The ^^'^Cs2 molecules 
in this experiment are bound by no more than E/h=5 
MHz and have a temperature and peak density of 250 
nK and 5 x 10^° cm~^, respectively. As the magnetic 
field is varied near i? =13 G the lifetime of the molecules 
rapidly changes, indicating two resonances. 

Heteronuclear molecules can be manipulated by static 
electric fields in addition to magnetic fields. The electric 
field Stark shifts the rotational levels of the molecule. 
These level shift s can then give rise to ele ctric Fesh- 
bach resonances ( Avdeenkov and Bohnl . [2002? ) . In atoms 
the levels can also be sufficiently Stark shifted to in- 
duce co Uisional resonances but rat her large fields are re- 
quired ( Marinescu and Youl . Il998l ). 

Figure [35] shows the results of a calculation on re- 
action rate coefficients of formald ehyde H2CO reactin g 
with OH to yield HCO and H2O (iHudson all [200i). 
Both H2CO and OH are in their lowest vibrational state 
of their ground electronic configuration. Multiple res- 
onances o ccur for electric fields up to 2 kV/cm. More 
recently, ( Tscherbul and Kremsl . |2008[ ) have studied re- 
action rates of LiF with H in the presence of electric 
fields. 

We can expect atom-molecule and molecule-molecule 
collisions to exhibit a rich variety of resonance phenom- 
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FIG. 47 Observation of Feshbach-like resonances in the col- 
lisions of ^^^Cs2 molecules. The figure shows the number of 
remaining ^'^''Cs molecules after fixed storage time in an op- 
tical trap as a function of external magnetic field. The two 
features near 13 G shown magnified in t he inset have bee n 
interpreted as tetramer resonances. From 
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FIG. 48 Calculated chemical reaction cross section of 
HaCO+OH^ HCOH-HaO at a collision energy of E/kg =1 
mK as a function of applied electric field. Strong-/weak-field 
seekers correspond to states of H2CO and OH that can be 
held in an electric trap where the center of the trap has the 
larger/smallest elect ric field strength, r espectively. Figure 
and model are from (|Hudson et al. 1 120061 ). 



ena in elastic, inelastic, and reactive collisions. Such phe- 
nomena are likely to become progressively more impor- 
tant to understand as sources of trap loss and for coherent 
control of molecular ensembles. 
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APPENDIX: Tables of selected resonances 

Tabic IIVI lists positions and properties of resonances 
for various species. The data is a combination of experi- 
mentally determined as well as theoretically derived val- 
ues. Most of the magnetic field locations are experimen- 
tally determined. Most of the widths A and background 
scattering lengths Obg are determined from theoretical 
calculations. Where unavailable from the literature we 
calculated values based on the best Born-Oppenheimer 
potentials obtained from the literature. For a complete 
list as well as error bars on the resonance locations the 
reader is referred to the literature. The notation defined 
in Sees. III. A. 31 and riI.B.4l has been used in the table. 

The table shows a richness in the kinds of resonances 
available for magnetic field values that are relatively eas- 
ily created in laboratories. Some of the resonances are 
very narrow with A on the order of a niG. Others are 
very broad with A larger than 100 G. The background 
scattering length can be either negative or positive, its 
absolute value ranging from a few tens to several thou- 
sands Bohr radii. The magnetic moment of the reso- 
nance state is always on the order of the Bohr magneton, 
which reflects the form of the Zeeman interaction. The 
partial wave of the resonance states ranges from zero to 
four {£c — 0...4). Finally, the resonances are character- 
ized in terms of their background scattering lengths Obg, 
their strengths Srcs (Sec. III.B.2[) . and the parameter C 
fSec. lII.C.5|) . 

For atomic cesium, a resonance location is given with a 
negative magnetic field value. This is not an experimen- 
tal value. Here Bq is determined from a fit of Eq. ^ to 
the slowly varying scattering length as shown in Fig. 
(Vogcls et a/.,'l99J) give the physical interpretation of a 
negative Bq, namely, taking B < corresponds to the 
case for B > with the spin projections of each atom 
reversed in sign. For the case of cesium, a negative mag- 
netic field in the aa channel corresponds to a positive 
field in the gg channel. 
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FIG. 49 Overview of selected Feshbach resonances in terms 
of widths A and strengths Sres . 



Figure SH] illustrates the rich variety of Feshbach res- 
onances in terms of their widths A and strengths Sres- 
Both parameters change over six orders of magnitude. 
Resonances with A > 1 G tend to be entrance-channel 
dominated (sios > !)• A notable exception is the ^Li 
737-G resonance mentioned in Sec. III.BTsI 
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TABLE IV Properties of selected Feshbach resonances. The first column describes the atomic species and isotope. The next 
three columns characterize the scattering and resonance state, which includes the incoming scattering channel (ch.), partial 
wave i", and the angular momentum of the resonance state 1^. This is followed by the resonance location Bo, the width A, 
the background scattering length Obg, the differential magnetic moment J/x, the dimensionless resonance strength Sros, the 



background scattering length in van der Waals units rbg = Obg/o, and the bound state parameter from Eq. (|52|l . Here ao is 
the Bohr radius and /is is the Bohr magneton. Definitions are given in Sec. [Ill The last column gives the source. A string 
"na" indicates that the corresponding property is not defined. For example Obg is not defined for p-wave scattering. 
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